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ABSTRACT 


Elastic  protein-based  polymers  comprised  of  repeating  pentapeptide  sequences,  GXGXP)n, 
exhibit  mechanical  resonances  that  have  been  observed  to  date  with  frequency  maxima  near  5 
MHz  and  3  kHz.  Because  the  3  kHz  resonance  is  in  the  middle  of  the  acoustic  frequency  range, 
the  purpose  here  is  to  substantiate  the  relevance  of  the  3  kHz  resonance  to  acoustic  absorption 
and  to  demonstrate  means  of  improving  mechanical  properties  for  the  sound  absorption 
application.  Previously  reported  loss  factor  data  in  the  100  Hz  to  10  kHz  range  is  substantiated 
by  relevant  but  distinctly  different  measurements  of  loss  shear  modulus  and  loss  permittivity. 
Furthermore  cross-linking  approaches  are  reported  that  result  in  increase  elastic  moduli  by  an 
order  of  magnitude  to  4  x  106  Pa  at  20%  strain  and  that  increase  break  stress  by  two  orders  of 
magnitude  to  1.3  x  107  Pa  while  exhibiting  break  stain  values  of  several  hundred  per  cent. 
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SUMMARY 

The  problem  addressed  in  this  effort  is  the  development  of  elastomeric  polypeptides  for  acoustic 
absorption  applications.  The  basic  problem  facing  this  objective  is  the  continued  fixation  by 
authorities  in  the  area  on  the  idea  that  elastic  materials  must  be  (must  have  the  properties  of) 
random  chain  networks.  The  staid  authorities  simply  will  not  entertain  the  possibility  that 
elastic  polymers  can  have  mechanical  resonances  and  that  such  resonances  can  exhibit  intense 
absorptions  (loss  factors)  in  the  acoustic  frequency  ranges.  A  quotation  due  to  Maurice 
Maeterlinck,  seems  particularly  appropriate,  “At  every  crossway  on  the  road  that  leads  to  the 
future,  tradition,  has  placed  against  each  of  us,  10,000  men  to  guard  the  past.”  This  report, 
buttressed  by  three  papers  in  the  APPENDIX,  is  believed  to  demonstrate  the  unique  opportunity 
that  elastic  protein-based  polymers  present  to  satisfy  sound  absorption  needs  of  the  Navy. 

Specifically,  elastic  protein-based  polymers  comprised  of  repeating  pentapeptide  sequences, 
(GXGXP)n,  exhibit  mechanical  resonances  that  have  been  observed  to  date  with  frequency 
maxima  near  5  MHz  and  3  kHz.  As  the  3  kHz  resonance  is  in  the  middle  of  the  acoustic 
frequency  range,  the  principal  result  of  this  report  substantiates  the  relevance  of  the  3  kHz 
resonance  to  acoustic  absorption.  Additional  results  demonstrate  means  of  improving  mechanical 
properties  for  the  sound  absorption  application. 

In  particular,  previously  reported  loss  factor  data  in  the  100  Hz  to  10  kHz  range  was  declared 
to  be  artifact  by  one  particularly  influential  man  of  the  10,000  men  that  “guard  the  past.”  The 
principal  result  of  this  report  is  that  loss  shear  modulus  and  loss  permittivity  data  substantiate 
the  intense  absorption  capacity  of  y-irradiation  cross-linked  (GVGIP)32o  centered  near  3  kHz 
with  moduli  in  the  range  of  5  x  106  Pa  and  with  the  potential  to  reach  107  Pa.  Furthermore, 
cross-linking  approaches  are  reported  that  result  in  materials  exhibiting,  in  the  most  favorable 
examples  so  far,  elastic  moduli  and  break  stresses  as  high  as  2  x  107  Pa  while  exhibiting  break 
strain  values  of  several  hundred  per  cent. 
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We  recommend  that  the  NAVY  give  serious  consideration  to  further  development  of  these 
unique  sound  absorbing  protein-based  materials,  and  that  those  who  believe  that  they  have  the 
instrumentation  with  which  to  measure  artifact  free  acoustic  absorption  of  elastomeric  materials 
find  it  in  their  wisdom  to  reconsider  doing  so. 


INTRODUCTION 


When  seeking  materials  suitable  for  low  frequency  sound  absorption,  compliant  or 
elastomeric  materials  become  a  natural  consideration.  Classical  elastomers,  or  rubbers,  are 
reasonable  and  commonly  used  materials,  but  they  have  a  fundamental  limitation.  They  are 


composed  of  random  chain 
networks  that  necessarily  exhibit 
low  intensity,  very  broad  and 
essentially  featureless  frequency 
dependences  of  relaxation  (or 
absorption).  Because  of  this,  the 
possibilities  of  tuning  or  designing 
classical  elastomers  to  maximize 
sound  absorption  over  a  desired 
frequency  range  are  limited. 

On  the  other  hand  certain  elastic 
protein-based  polymers  can  exhibit 
good  elastic  moduli  and  yet  contain 
dynamic  structural  regularities  that 
exhibit  mechanical  modes.  In  fact 
elastic  protein-based  polymers  can 
be  designed  to  exhibit  mechanical 
resonances  with  high  loss  factors 
over  interesting  frequency  ranges. 

Preliminary  acoustic  absorption 
studies  on  elastic  protein-based 
polymeric  materials: 

The  elastic  protein-based 
polymer,  (glycyl-valyl-glycyl- 
isoleucyl-prolyl)n  abbreviated  as 
(GVGIP)n,  emerged  as  a  particularly 
interesting  composition;  it  exhibited 
a  large  hysteresis  in  macroscopic 
stress-strain  studies  and  in  atomic 
force  microscopy(AFM)-based 
single-chain  force-extension 
studies . 1  Therefore,  (G  V  GIP) 


Figure  1:  A.  Acoustic  absorption  of  20  Mrad  y-irradiation 
cross-linked  (GVGIP)260 
B.  5MHz  dielectric  relaxation  of  poly(GVGIP) 

(imaginary  part,  6”) 
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might  be  expected  to  exhibit  substantial  loss  factors  as  required  for  good  acoustic  absorption. 

Accordingly,  two  y-irradiation  (20  Mrad)  cross-linked  cylinders  of  X20-(GVGIP)260 
were  prepared  with  diameters  of  2.5  cm  and  lengths  of  0.5  cm  and  1.0  cm.  Lev  Sheiba  and  Jack 
Dea  examined  the  absorption  properties  of  these  cylinders  at  the  SPAWAR  (Spatial  Warfare) 
Systems  Center  in  San  Diego.  The  results  are  given  in  Figure  1A  where  the  elastic  protein-based 

polymer  is  compared  to  natural  rubber  and  an  absorption  enhanced  polyurethane.  As  to  be 
expected  for  the  random  chain  networks  of  natural  rubber  and  polyurethane,  both  exhibited  broad 
featureless  dependences  on  frequency  and  relatively  low  intensity  loss  factors  (sound 
absorption/unit  volume).  On  the  other  hand,  the  elastic  protein-based  polymer  exhibited  a 
localized  relaxation  that  grew  in  intensity  as  the  temperature  increased,  just  as  had  been  observed 
in  the  dielectric  relaxation  spectra,  specifically  in  the  imaginary  part  of  the  dielectric  permittivity, 

of  (GVGIP)n  at  5  MHz  (See  Figure  IB).2 


Molecular  structure  of  (GXGXP)-type  elastic  protein-based  polymers: 

(GXGXP)-type  protein-based  elastomers  form  (3-spiral  structures  described  as  a  series  of  (3- 
tums  arranged  in  helical  array  with  dynamic  suspended  segments  between  the  (3-tums.  For 
discussion  of  the  single  secondary  structural  feature  of  the  (3- turn,  the  repeating  conformational 
feature  is  considered  in  terms  of  the  sequence  permutation,  X1P2G3X4G5.  With  this  residue 
numbering,  the  (3-tum  is  a  ten-atom  hydrogen-bonded  ring  involving  the  C-0  of  the  X  residue 
preceding  the  proline  residue  and  the  NH  of  the  X4  residue.  The  (3-spiral  structure  in  band 
representation  of  the  backbone  is  shown  on  the  left  side  of  Figure  2,  and  cross-eye  stereo  views  of 
the  detailed  molecular  structure  are  seen  in  side  and  spiral  axis  views  in  the  center  structures  of 


Figure  2:  Molecular  structure  and  acoustic  absorption  cartoon  for  elastic  protein- 
based  polymers  of  the  GXGXP-type  repeating  pentapeptide  sequence. 
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Figure  2.  At  the  molecular  level,  there  occur  suspended  segments  connecting  the  (3-tums  that 
contain  peptide  moieties  capable  of  large  amplitude  rocking  motions.  The  suspended  segments 
exhibit  large  amplitude  torsional  oscillations  that  become  damped  on  extension.  Thus,  the 
suspended  segments  that  link  the  |3-tums  function  as  entropic  molecular  springs. 

The  P-tums  themselves  provide  dynamically  active  surface  elements,  effectively  forming  a  series 
of  plates  on  the  surface  of  the  wide  dynamic  helix  that  is  the  P-spiral.  A  macroscopic-classical 
analog  would  be  equivalent  to  springs  connecting  the  (5-tum  plates.  The  P-tum  surface  elements 
may  be  thought  of  as  molecular  acoustic  veins  attached  to  springs  that  are  tunable  to  resonate  at 
desired  frequencies.  A  cartoon  demonstrating  these  structural  features  is  shown  as  the  right-most 
illustration  of  Figure  2.  The  four  spiral  structures  are  separated,  such  that  when  the  two  central 
detailed  structures  are  viewed  in  stereo,  the  backbone  band  representation  and  the  cartoon  also 
overlap  with  detailed  molecular  structures. 

Mechanism  of  elasticity  and  mechanical  resonances  in  repeating  peptide  sequences: 

By  means  of  AFM  (atomic  force  microscopy),  single-chain  force-extension  curves  have  been 
obtained  on  Cys-(GVGIP)nx260-Cys  and  Cys-(GVGVP)nx25l-Cys,  where  Cys  is  the  cysteinyl 
residue  with  a  sulfhydryl  for  attachment  to  surfaces.1  In  the  experiment,  one  end  of  the  single  chain 
is  attached  to  the  cantilever  tip  and  the  other  end  to  a  surface.  Force  is  measured  as  a  function  of 
chain  length,  and  elastic  moduli  of  several  pN/chain  have  been  obtained.  In  this  experiment,  a  single 
chain  held  fixed  at  both  ends  gives  rise  to  entropic  elastic  force.  Since  total  molecule  rotations  and 
translations  are  not  possible,  the  decrease  in  entropy,  that  gives  rise  to  an  entropic  elastic  force  on 
extension  to  a  particular  fixed  length,  must  come  from  the  damping  of  internal  chain  dynamics. 
Therefore,  it  is  established  that  entropic  elasticity  can  occur  without  random  chain  networks  and 
without  a  random,  or  Gaussian,  distribution  of  end-to-end  chain  lengths. 

Importantly,  the  intensity  of  the  5  MHz  dielectric  relaxation  of  (GVGVP)n  can  be  calculated 
from  the  molecular  dynamics  of  the  structure  of  the  repeating  unit  in  Figure  2  using  Onsagers’s 
equation  for  polar  liquids.4’5  Quoting  from  the  McGraw-Hill  Concise  Encyclopedia  of  Science 
and  Technology,  2nd  Edition,  1989,  page  1599,  “When  a  mechanical  or  acoustical  system  is  acted 
upon  by  an  external  periodic  driving  force  whose  frequency  equals  a  natural  free  oscillation 
frequency  of  the  system,  the  amplitude  of  the  oscillation  becomes  large  and  the  system  is  said  to  be 
in  a  state  of  resonance.”  The  repeating  structural  units  of  Figure  2  effectively  resonate,  and  the 
absorption  (loss  factor)  peaks  in  the  relaxation  spectra  may  be  called  mechanical  resonances. 
Furthermore,  the  occurrence  of  mechanical  resonances  in  elastic-protein-based  polymers  makes 
the  classical  theory  of  rubber  elasticity  only  marginally  relevant  to  these  polymers  and  the  use  of 
related  considerations  in  treatments  of  acoustic  absorption  by  such  polymers  would  also  be 
suspected  of  similar  limitations.  Finally,  the  occurrence  of  mechanical  resonances  in  these  elastic 
protein-based  polymers  presents  unique  opportunities  for  designing  materials  with  acoustic 
absorption  properties  that  can  be  improved  over  current  elastomeric  materials. 

Acoustic  absorption  data  claimed  to  be  artifact: 

In  spite  of  the  occurrence  of  the  same  temperature  dependent  phenomena  near  5  MHz 
(Figure  IB)  and  near  3  kHz  (Figure  1A)  using  the  same  elastic  protein-based  polymer 
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composition,  the  acoustic  absorption  data  of  Figure  1 A  on  the  elastic  protein-based  polymer  has 
been  declared  to  be  an  artifact  by  academic  professionals  concerned  with  acoustic  absorption 
applications  for  the  Navy.  This  is  in  spite  of  the  fact  that  the  similarly  determined  data  for 
natural  rubber  and  polyurethane  did  not  appear  to  be  in  question.  At  the  same  time  that  the  data 
from  SPA  WAR  has  been  declared  to  be  artifact,  the  party  making  the  claim  refuses  to  measure,  in 
their  “uniquely  artifact-free”  instrumentation,  the  base  polymer  composition  as  well  as  additional 
analogues  designed  to  extend  further  the  acoustic  absorption  properties  of  this  family  of  elastic 
protein-based  polymers.  Because  of  this,  independent  approaches  have  been  sought  and  they 
have  been  found  to  evaluate  the  significance  of  the  data  in  Figure  1  A. 

Content  of  this  report: 

This  report  presents  a  family  of  ten  elastic  protein-based  polymers  that  have  been  designed 
to  achieve  even  more  potent  acoustic  absorption  materials.  They  have  been  prepared  using 
recombinant  DNA  technology  and  transformation  of  E.  coli  and  have  been  expressed  in  large 
quantities  through  E.  coli  fermentation.  Improved  mechanical  properties  -  such  as  essential  work 
of  fracture,  elastic  modulus,  break  stress  and  break  strain  -  have  been  achieved  by  obtaining 
higher  molecular  weight,  monodisperse  polymers  and  by  improved  cross-linking  approaches. 

The  design  for  improved  chemical  cross-linking  utilized  the  concept  of  apolar-polar  repulsive  free 
energies  of  hydration.  Importantly  for  developing  improved  acoustic  absorption  materials,  using 
the  base  composition  of  the  family,  X20-(GVGIP)32o,  a  combination  of  loss  shear  modulus,  G”, 
measurements  and  loss  permittivity,  8”,  data  have  been  obtained  that  substantiate  the  existence 
of  an  intense  temperature  dependent  acoustic  absorption,  that  is  a  mechanical  resonance,  centered 
near  3kHz. 

METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

Production  of  polymers  by  recombinant  DNA  technology: 

The  list  of  elastic  protein-based  polymers  of  Table  1  were  prepared  by  recombinant  DNA 
technology.  This  involves  the  construction  of  a  monomer  gene  from  two  chemically  synthesized, 
overlapping  DNA  sequences  that  were  made  into  the  complete  double-stranded  DNA  by  the 
polymerase  chain  reaction  (PCR).  Trace  amounts  of  the  monomer  genes  so  prepared  were 
inserted  into  an  appropriate  plasmid  and  a  strain  of  E.  coli  was  transformed  that  produced  a  large 
number  of  copies  of  the  monomer  gene.  The  isolated,  relatively  large  quantity  of  monomer  gene 
was  polymerized  (in  terms  molecular  biology,  concatemerized  or  concatenated)  and  ultimately 
used  to  transform  a  strain  of  E.  coli  that  was  suitable  for  expression  of  the  elastic  protein-based 
polymer.6  Table  I  gives  the  size  of  the  gene  for  the  resulting  elastic  protein-based  polymer  in 
terms  of  the  number  of  repeating  peptide  units  and  the  amount  of  polymer  prepared  for 
subsequent  characterization. 

Preparation  of  Polymers  I’  through  VI’  was  for  chemical  cross-linking  efforts,  whereas 
Polymers  VH’  through  XI’  and  XV’  and  XVII’  were  prepared  for  systematic  change  in 
hydrophobicity  and  for  inclusion  of  a  sulfate  to  introduce  the  effect  on  acoustic  absorption  of 
MgSC>4  association/dissociation  in  polymers  of  different  hydrophobicities.  Finally,  Polymer 
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XVII’  was  specifically  prepared  for  its  higher  molecular  weight  and  for  a  more  monodisperse 
molecular  weight  than  had  been  prepared  previously. 


Table  1: 

prepared 

pentar 

grams 

perpo 

Polymer  I' 

[(GVGIP  GEGIP  GVGIP)3]n 

17 

240 

34 

108 

Polymer  IF 

[(GVGIP  GVGIP  GEGIP  GVGIP  GVGIP  GVGIP)]n 

58 

276 

Polymer  III' 

[(GEGIP  GVGIP  GEGIP  GVGIP  GVGIP  GVGIP)]„ 

31 

180 

Polymer  IV' 

[(GVGIP  GKGIP  GVGIP)3]„ 

0 

Polymer  V' 

[(GVGIP  GVGIP  GKGIP  GVGIP  GVGIP  GVGIP)]n 

20 

252 

20 

144 

Polymer  VT 

[(GKGIP  GVGIP  GKGIP  GVGIP  GVGIP  GVGIP)]„ 

0 

Polymer  VII' 

[(GVGIP)21(GYGIP)]n, 

60 

242 

Polymer  IX' 

[(GVGIP),  ,(GYGIP)]n, 

40 

216 

Polymer  XI' 

[(GV  GIP)g(G  Y  GIP  )]„, 

47 

288 

Polymer  XIII' 

[(GVGIP)5(GYGIP)]„, 

30 

240 

Polymer  XV' 

{[(GVGIP)2(GYGIP)]3}„, 

trace 

135 

Polymer  XVII' 

[(GVGIP), 0]n 

247 

320 

Polymer  VIII' 

[(GVGIP)2,(GY  {S04=}GIP)]n 

19 

242 

Polymer  X' 

[(GVGIP),  ,(GY{S04=}GIP)]n 

19 

216 

Cross-linking  of  polymers: 

Preparation  of  y-irradiation  cross-linked  matrices:  In  preparation  for  y-irradiation  cross- 
linking,  the  polymer  is  dissolved  in  water  at  low  temperature.  On  raising  the  temperature  above 
that  of  the  inverse  temperature  transition  for  hydrophobic  association,  phase  separation  occurs. 
The  phase-separated  state  is  then  exposed  to  20  Mrad  of  y-irradiation  from  a  cobalt-60  source. 

Preparation  of  chemically  cross-linked  fibers:  Chemical  cross-linking  was  achieved  by 
extrusion  of  solutions  of  E-  and  K-  pairs  of  polymers  at  equal  concentrations  of  functional 
groups  into  a  saturated  solution  at  50°C  of  water  soluble  carbodiimide  (EDC,  l-(3- 
dimethylaminopropyl)-3-ethylcarbodiimide)  to  form  amide  bonds  between  the  carboxyl  of 
glutamic  acid  (E)  residues  and  the  e-amino  groups  of  lysine  (K)  residues.7  When  in  an  adequately 
hydrophobic  elastic  protein-based  polymer,  the  charged  carboxylate  and  amino  functions 
experience  a  driving  force  for  ion-pairing.  The  force  arises  due  to  an  apolar-polar  repulsive  free 
energy  of  hydration,  AGap,  that  is  measurable  from  the  hydrophobic-induced  pKa  shifts. 
Cross-linking  occurred  on  extrusion  to  form  fibers  and  also  in  a  manner  to  result  in  disks.  When 
using  microbially  prepared  elastic  protein-based  polymers  with  functional  groups  occurring  every 
thirty  residues  and  with  the  favorable  AGap,  materials  of  exceptional  mechanical  properties  are 
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obtained,  for  example  fibers  exhibited  elastic  moduli  and  break  stresses  that  are  one  to  two  orders 
of  magnitude  greater  than  obtained  with  y-irradiation  cross-linking. 

Determination  of  shear  moduli: 

The  instrument  for  determination  of  shear  moduli  was  a  Rheometric  Scientific  dynamic 
mechanical  thermal  analyzer,  model  DMTA  V.  Round  shear  sandwich  geometry  was  used.  The 
instrument  was  inverted  so  that  the  sandwich  fixtures  and  sample  were  in  water.  A  water- 
jacketed  1000  mL  Pyrex  cylinder  supplied  by  Rheometric  Scientific  allowed  control  of 
temperature  with  a  circulating  temperature  bath.  A  sinusoidal  linear  shear  was  applied  by 
moving  a  flat  plate  between  two  identical  disk-shaped  samples  over  a  specified  range  of 
frequencies.  The  two  identical  disk-shaped  samples  were  sandwiched  between  the  moving  plate 
and  two  12  mm  diameter  plates  (called  studs)  fastened  to  a  frame.  The  dynamic  frequency 
sweeps  to  obtain  the  loss  shear  modulus,  G”,  from  0.16  Hz  to  318  Hz  were  reported  in  log 
scale.9  For  our  purposes  Initial  static  force  was  0.05  N.  Sample  size  was  12  mm  diameter  and 
0.7  mm  thickness.  The  sample  was  equilibrated  at  40°C  for  12  hours  prior  to  starting  the  series 
of  measurements.  Shear  moduli  were  measured  from  high  to  low  temperature  with  an 
equilibration  time  of  2  hours  at  each  temperature.  The  sequence  of  measurements  was  40°C, 
30°C,  25°C,  and  15°C. 

Determination  of  low  frequency  dielectric  relaxation: 

Low  frequency  dielectric  relaxation  data  was  obtained  with  a  Rheometric  Scientific  dielectric 
thermal  analyzer  (DETA)  system,  model  MK-II.  The  unit  consists  of  a  temperature 
programmer,  dielectric  thermal  analyzer,  and  water-jacketed  inverted  can-shaped  enclosure,  called 
a  furnace.  Within  the  inverted  can-shaped  enclosure  is  a  pair  of  disk-shaped  parallel  plates 
between  which  the  sample  is  placed.  The  sample  was  sealed  between  0.01mm  Teflon  sheets 
(Goodfellow,  Cambridge,  England)  and  the  sealed  packet  was  placed  between  the  two  15mm- 
diameter  parallel  plates.  An  alternating  electric  field  is  applied  to  the  sample  by  means  of  a 
sinusoidal  voltage  and  an  out  of  phase  current  is  measured  to  obtain  the  loss  permittivity,  8“. 

An  external  temperature  bath  was  used  for  temperature  control.  Sample  size  for  X  -(GVGIP)320 
was  15mm  diameter  and  0.6  mm  thickness.  The  system  was  equilibrated  overnight  at  40°C. 
Temperature  was  decreased  in  5-degree  increments  with  2  hours  of  equilibration  between 
temperature  changes  and  3  sets  of  measurements  were  made  from  20  Hz  to  100  kHz  at  each 
temperature. 

Basis  for  correlation  between  loss  shear  modulus  and  loss  permittivity: 

When  a  polymer  exhibits  a  maximum  in  the  imaginary  part  of  the  dielectric  permittivity  (the 
loss  permittivity,  £")  at  frequencies  less  than  200  Hz,  it  becomes  possible  to  make  comparisons 
with  the  frequency  dependence  of  shear  moduli  and  most  specifically  with  the  loss  shear 
modulus,  G”.  This  has  been  done  for  polypropylene  diol,  also  called  poly(oxypropylene),  where 
there  is  reported  a  near  perfect  superposition  of  the  frequency  dependence  of  the  normalized  loss 
shear  modulus  with  that  of  the  normalized  loss  permittivity  as  reproduced  in  Figure  3.  ’  The 
acoustic  absorption  frequency  range  of  interest  here  is  100  Hz  to  10  kHz,  yet  present 
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Figure  3:  Superposition  of  loss 
permittivity  and  loss  shear  modulus 
for  polypropylene  diol.  Adapted 
with  permission  from  references  10 
and  1 1 . 


macroscopic  loss  shear  modulus  data  can  be 
determined  at  most  up  to  a  few  hundred  Hz. 
Nonetheless,  for  X20-(GVGIP)320  there  is  a  maximum 
in  loss  permittivity,  8“,  near  3  kHz  that  develops  on 
raising  the  temperature  through  the  temperature  range 
of  the  inverse  temperature  transition.  With  the  width 
of  the  loss  permittivity  curve  a  distinct  set  of  curves 
as  a  function  of  temperature  become  defined  as  low 
as  the  100  to  200  Hz  frequency  range.  Accordingly, 
if  these  loss  permittivity  curves  are  sufficiently 
distinctive  and  compare  favorably  with  the  exactly 
corresponding  set  of  loss  shear  modulus  curves,  the 
melding  of  the  two  sets  of  data  in  the  overlap  region 
provide  for  calibration  of  the  loss  permittivity  curves 
in  terms  of  loss  shear  modulus  of  recognized 
relevance  to  acoustic  absorption  data. 

Determination  of  essential  work  of  fracture: 

Samples  of  X20-(GVGIP)320  were  run  on  an  in 
house  custom-built  vertical  stress-strain  apparatus.  A 
water-jacketed  chamber  allowed  grips  and  sample  to 
be  submerged  at  all  times  with  temperature  control. 
Sample  size  was  4x5x0.5  mm.  Samples  were  double 
notched  with  a  fresh  razor  blade  and  stretched  at  a 

speed  of  2.2  mm/min  until  fractured.12  Several 
different  ligament  lengths  (distance  between  notches) 
were  used  and  work  of  fracture  was  plotted  against 
ligament  length. 


RESULTS  AND  DISCUSSION 

Unique  temperature  and  frequency  dependence  of  loss  shear  modulus  data: 

It  is  generally  recognized  that  the  loss  shear  modulus,  G”,  data  in  the  1  Hz  to  200  Hz 
frequency  range  is  relevant  to  acoustic  absorption  in  the  100  Hz  to  10  kHz  frequency  range. 
Commonly,  data  obtainable  at  lower  frequencies,  generally  100  Hz  or  less,  has  been  extrapolated 
into  the  acoustic  frequency  range  in  order  to  evaluate  candidates  for  good  acoustic  absorption 
properties.  Such  an  extrapolation  requires  an  assumption  about  the  frequency  dependence. 
While  this  becomes  a  reasonable  practice  for  classical  elastomers  that  occur  as  random  chain 
networks  with  the  broad,  featureless  frequency  dependence  of  natural  rubber  and  polyurethane 
seen  in  Figure  1  A,  it  is  not  a  reasonable  assumption  for  elastic  protein-based  polymers  that 
exhibit  mechanical  resonances  as  seen  in  Figure  IB  for  the  loss  permittivity  of  (GVGIP)n. 
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The  variable  temperature  dependence  of  the  steepness  of  curvature  for  the  loss  shear  modulus 
curves  reported  in  Figure  4  causes  concern  for  an  approach  that  extrapolates  to  higher  frequencies 


A. 


Figure  4:  A.  Discs  of  X20-(GVGIP)32o .  B.  Loss  shear  modulus,  G”  vs.  frequency  for  (GVGIP)32o 

without  a  better  understanding  of  the  frequency  dependence.  Because  of  the  correlation  of  loss 
shear  modulus  and  loss  permittivity  data  for  polypropylene  diol  of  Figure  3,  loss  permittivity  data 
in  the  acoustic  frequency  range  has  particular  significance  for  elastic  protein-based  polymers  that 
exhibit  a  high  temperature  dependence  for  the  loss  permittivity  as  seen  in  Figure  IB.  In  fact,  as 
shown  immediately  below,  low  frequency  dielectric  relaxation  data  provides  a  means  more  effective 
than  computational  extrapolation  with  which  to  evaluate  acoustic  absorption  properties  of  elastic 
polymers  that  exhibit  mechanical  resonances. 


Melding  of  loss  shear  modulus  and  loss  permittivity  data  to  evaluate  acoustic  absorption 
potential  of  elastic  protein-based  polymers: 

The  temperature  dependence  of  the  loss  shear  modulus  data  of  Figure  4  is  quite  distinctive. 
At  200Hz  the  temperature  dependence  is  one  of  increasing  G”  on  increasing  the  temperature  from 
15  to  25  °C,  but  then  the  next  larger  G”  occurs  at  40°C  with  the  value  of  G”  at  200  Hz  for  the 
30°C  data  being  more  than  twice  as  large.  On  the  other  hand  at  5  MHz  in  Figure  IB,  the 
increase  in  intensity  of  the  loss  permittivity  continues  to  increase  on  raising  the  temperature  from 
10  to  60°C. 

Nonetheless,  the  merging  of  loss  shear  modulus,  G”,  and  loss  permittivity,  8”,  data  for  the 
overlapping  frequencies  in  Figure  5  is  truly  remarkable.  Explanation  of  the  distinctive  40  C  curve 
now  becomes  obvious.  The  explanation  arises  due  to  the  occurrence  of  the  temperature  interval  of 
the  inverse  temperature  transition8  being  from  10  to  30°C.  During  the  inverse  temperature 
transition  of  hydrophobic 
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Figure  5:  Melding  of  G”  and  £”  data  for  (GVGIP)32o 

based  polymer  can  be  chosen  such  that  the  temperature  interval 
intense  absorption  would  exist  at  temperatures  above  5°C. 


folding  into  the  structure  of  Figure 
2,  the  frequency  maximum  of  the 
loss  permittivity  peak  decreases  as 
the  elastic  protein-based  polymer 
folds  into  a  more  regular  structure. 
Once  the  inverse  temperature 
transition  is  complete,  the  normal 
increase  in  frequency  with  increase 
in  temperature  ensues. 

A  coarse  calibration  of  the 
higher  frequency  loss  permittivity 
data  becomes  possible  by  using  the 
overlap  with  the  loss  shear 
modulus  values.  It  seems 
reasonable  to  conclude  that  the  loss 
shear  modulus  values  would  reach 
5xl06  Pa  or  more  near  3  kHz. 

The  comparison  also  suggests  that 
the  intensity  of  the  peak  in  Figure 
1 A  would  almost  double  on 
continuing  from  20  to  40°C.  It 
should  also  be  appreciated  that  the 
composition  of  the  elastic  protein- 
is  below  5°C  in  which  case  the 


Mechanical  properties  of  cross-linked  fibers: 

When  utilizing  the  apolar-polar  repulsive  free  energy  for  hydration,  AGap,  the  free  energy 
gained  on  ion-pairing  aligns  and  orders  cross-linking  groups,  E  and  K.  This  is  demonstrated  for 
the  polymer  pair. 


Polymer  i: 

Polymer  x’: 

and  for  the  pair, 
Polymer  ii: 

Polymer  xi’: 


(GVGVP  GVGVP  GEGVP  GVGVP  GVGVP  GVGVP)36(GVGVP) 
(GVGVP  GVGVP  GKGVP  GVGVP  GVGVP  GVGVP)22(GVGVP) 

(GVGVP  GVGFP  GEGFP  GVGVP  GVGVP  GVGVP)40(GVGVP) 
(GVGVP  GVGFP  GKGFP  GVGVP  GVGVP  GVGVP)22(GVGVP) 


E/OF 

K/OF, 

E/2F 

K/2F. 
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Micrographs  of  the  fibers  formed  and  represented  mechanical  properties  are  given  in  Figure  6. 
Among  the  several  experiments,  values  of  the  elastic  modulus,  E,  at  20%  of  greater  than  10  Pa  have 
been  obtained  and  break  stresses  greater  than  1.5  x  107  Pa  have  also  been  obtained  for  both  fiber 
compositions.  These  values  are  two  orders  of  magnitude  greater  than  obtained  by  means  of  y 
irradiation  cross-linking.  Generally,  the  more  hydrophobic  phenylalanine  (Phe,  F)  containing 

polymers  gave  the  higher 
values. 


Essential  work  of  fracture 

The  stress-strain  curves  given 
in  terms  of  the  work  in  J/m2 
versus  extension  to  the  point  of 

breakage  is  seen  in  Figure  7A. 
When  the  specific  work  of 
fracture  is  plotted  as  a 
function  of  ligament  length 
resulting  from  notching  each 
side  of  the  elastic  strip, 
extrapolation  to  the  intercept 
gives  the  essential  work  of 
fracture.  The  results  for  both 
X20-(GVGIP)260  and 
X20-(GVGIP)32o  are  shown  in 
Figure  7B.  While  the  gene  for 
(GVGIP)260  indicates  a  chain 
length  of  260  pentamers,  the 
expressed  protein-based 
polymer  contained  a  smear  of 
chains  of  lower  molecular 
weights.  In  making  a  new 
monomer  gene  encoding  for  10 
pentamers,  the  same  codons 
were  used  for  the  two  G 
residues  within  the  GVGIP 


Figure  6:  Mechanical  properties  of  chemically  cross-linked  fibers. 


pentamer  but  different  codons 
for  G  were  used  between  the 


pentamers.  On  multimerization  a  more  stable  higher  molecular  weight  gene  product  was 
expressed,  namely  (GVGIP)32o-  On  20  Mrad  y-irradiation  cross-linking 
an  elastic  matrix  formed  with  twice  the  essential  work  of  fracture  that  was  more  consistent  with 
chemically-synthesized  poly(GVGIP).  Thus,  a  tougher  matrix  was  obtained  that  would  be  more 
suitable  for  acoustic  absorption  applications. 
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CONCLUSIONS 

Elastic  protein-based  polymers  of  the  (GXGXP)-type  exhibit  mechanical  resonances  that 
open  a  new  vista  for  the  development  of  unique  acoustic  absorption  materials.  Previously 
determined  temperature  dependent  loss  factor  (absorption/unit  volume)  curves  for 
X20-(GVGIP)260  in  the  100  Hz  to  10  kHz  range  have  been  dismissed  as  artifact  in  spite  of  the 
fact  that  exactly  parallel  behavior  as  a  function  of  temperature  is  observed  for  a  mechanical 

resonance  of  (GVGIP)n  at  5  MHz  by  means 
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of  loss  permittivity  measurements.  Melding 
of  loss  shear  modulus  and  loss  permittivity 
data  confirmed  the  presence  of  an  intense 
mechanical  resonance  in  the  acoustic 
absorption  range  for  matrices  of 
X20-(GVGIP)32o-  The  cylindrical  sample  of 
Figure  8  stands  ready  for  acoustic 
absorption  measurements  by  instrumentation 
claimed  to  be  uniquely  free  of  artifact. 

Using  an  understanding  of  the  apolar- 


Figure  8:  Cylinder  of  20  Mrad  y- 

irradiation  cross-linked  (GVGIP)32o  prepared 
in  a  form  suitable  for  acoustic  absorption 
measurements  by  the  “uniquely  artifact  free” 
instrumentation. 


Ligament  Length  (mm) 

Figure  7:  Essential  work  of  fracture  for 
X20-(GVGIP)260  and  X20-(GVGIP)320. 


polar  repulsive  free  energy  of  hydration,  AGap,8  elastic  protein-based  polymers  can  be  prepared 
with  impressive  mechanical  properties  of  elastic  moduli  and  break  stresses  up  to  2  x  10  Pa  and 
break  strains  of  several  hundred  percent. 

RECOMMENDATIONS 

Based  on  the  data  presented  in  this  Final  Report,  we  respectfully  recommend  that  the  Office  of 
Naval  Research  and  the  Navy  consider  elastic  protein-based  materials  for  their  sound  absorption 
needs.  One  possible  area  of  application  would  be  hearing  protection.  To  the  best  of  our 
knowledge  the  property  of  mechanical  resonances  with  high  loss  factors  in  the  acoustic 
absorption  range,  as  exhibited  by  (GXGXP)-type  protein-based  elastomers,  is  simply  not 
available  in  random  chain  network  elastomers. 

ACKNOWLEDGMENTS:  The  authors  gratefully  acknowledge  the  Office  of  Navy  Research 
under  Contract  Number  N00014-00-C-0178  for  the  support  of  these  materials  developments. 
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During  the  last  half  century,  identification  of  an  ideal  (predominantly  entropic)  protein  elastomer  was 
generally  thought  to  require  that  the  ideal  protein  elastomer  be  a  random  chain  network.  Here,  we  report 
two  new  sets  of  data  and  review  previous  data.  The  first  set  of  new  data  utilizes  atomic  force  microscopy 
to  report  single-chain  force-extension  curves  for  (GVGVP)251  and  (GVGIP)260,  and  provides  evidence 
for  single-chain  ideal  elasticity.  The  second  class  of  new  data  provides  a  direct  contrast  between  low- 
frequency  sound  absorption  (0.1-10  kHz)  exhibited  by  random-chain  network  elastomers  and  by  elastin 
protein-based  polymers. 

Earlier  composition,  dielectric  relaxation  (1-1000  MHz),  thermoelasticity,  molecular  mechanics  and 
dynamics  calculations  and  thermodynamic  and  statistical  mechanical  analyses  are  presented,  that  combine 
with  the  new  data  to  contrast  with  random-chain  network  rubbers  and  to  detail  the  presence  of  regular 
non-random  structural  elements  of  the  elastin-based  systems  that  lose  entropic  elastomeric  force  upon 
thermal  denaturation. 

The  data  and  analyses  affirm  an  earlier  contrary  argument  that  components  of  elastin,  the  elastic  protein 
of  the  mammalian  elastic  fibre,  and  purified  elastin  fibre  itself  contain  dynamic,  non-random,  regularly 
repeating  structures  that  exhibit  dominantly  entropic  elasticity  by  means  of  a  damping  of  internal  chain 
dynamics  on  extension. 

Keywords:  atomic  force  microscopy;  acoustic  absorption;  dielectric  relaxation;  thermoelasticity; 

entropic  elasticity;  p-spiral 


1.  INTRODUCTION 

(a)  Definition  of  ideal  or  perfect  elasticity 
Ideal  elasticity  is  the  property  whereby  the  energy 
expended  in  deformation  of  the  elastomer  recovers  com¬ 
pletely  on  the  removal  of  the  deforming  force.  As  the 
energy  expended  in  deformation  is  given  by  the  area  under 
the  force,  /,  versus  the  increase  in  length,  AL,  curve,  a 
perfectly  reversible  force-extension  curve  means  complete 
recovery  on  relaxation  of  the  energy  expended  on  defor¬ 
mation.  Therefore,  ideal  elastomers  exhibit  perfectly 
reversible  force-extension  curves. 

Perhaps  our  earliest  perspective  of  the  mechanism 
underlying  ideal  elasticity  comes  from  a  fundamental 
observation  concerning  rubber  elasticity.  In  the  mid- 19th 
century.  Joule  and  Thomson  noted  a  quantitative  corre¬ 
lation  between  the  increase  in  temperature  of  the  elasto¬ 
mer  due  to  stretching  and  the  increase  in  force  due  to 
increasing  the  temperature  (Flory  1968).  Thermodyn¬ 
amics  provides  for  the  analysis  underlying  this  correlation, 
and  the  Boltzmann  relationship  provides  the  bridge 
between  experimental  thermodynamic  quantities  and  a 
statistical  mechanical  description  of  molecular  structures. 
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Continuing  qualitatively  with  the  Joule  and  Thomson 
correlation,  heat  produces  motion,  and  the  energy  rep¬ 
resented  by  heat  distributes  into  the  various  available 
degrees  of  freedom  in  the  chain  molecules  comprising  the 
elastomer.  Accordingly,  the  release  of  heat  on  stretching 
correlates  with  a  loss  of  motion.  By  means  of  statistical 
mechanics,  the  loss  of  motion  is  seen  as  a  decrease  in 
entropy  on  extension.  In  addition,  should  a  solvent  be 
essential  for  elasticity,  this  requires  explicit  consideration. 

However,  there  is  more  to  elastomeric  force  than 
entropy  changes  arising  from  changes  in  mobility.  There 
is  also  internal  energy,  and  a  proper  understanding  of  elas¬ 
ticity  requires  the  delineation  of  the  internal  energy  and 
entropy  components  of  the  elastomeric  force.  As  the 
internal  energy  component  of  the  elastic  force  increases, 
irreversible  processes,  such  as  chain  breakage,  become 
more  probable.  Therefore,  greater  durability  results  as 
elastomers  become  more  dominantly  entropic  and  dura¬ 
bility  becomes  a  feature  of  an  ideal,  or  more  perfect,  elas¬ 
tomer. 


(b)  Delineation  of  internal  energy  and  entropy 
components  of  elastomeric  force 
The  delineation  of  internal  energy  and  entropy  compo¬ 
nents  of  elastomeric  force  can  begin  with  the  definition  of 
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the  Helmholtz  free  energy.  A,  also  referred  to  as  the  maxi¬ 
mal  work  function, 

A  =  E  -  TSs  (1.1) 

where  E  is  the  internal  energy,  T  is  the  absolute  tempera¬ 
ture  (°K),  and  S'  is  the  entropy  for  the  system.  The  differ¬ 
ential  of  the  maximum  work  function,  d  A ,  can  be  written 
to  include  the  work  done  on  an  elastomer  by  application 
of  a  force,  /,  over  the  change  in  length,  dL,  i.e. 

d  A  =  -  PdV  -  SdT  +  /dL,  (1.2) 

where  P  is  the  pressure,  and  V  is  the  volume.  Writing  the 
partial  differential  of  A  in  equation  (1.2)  with  respect  to 
length  at  constant  V,  T  and  composition,  k,  gives 

0 AfdL)ViTtn  =  (dE/dL) v>r>n  -  TQS/dQy, r>M,  (1 .3) 

but  by  equation  (1.2),  ( dA/dL)ViTj„=fi  such  that 

/=  (dE/dL) VtT>n  -  T(dS/dL)v,Tt„ .  (1.4) 

Accordingly,  the  force  is  seen  to  comprise  two  compo¬ 
nents,  an  internal  energy  component,  and  an  entropy 
component,  fSy  i.e. 

/=/*+/*.  (1.5) 

Now  it  becomes  useful  to  find  an  expression  that  will  allow 
experimental  estimation  of  the  relative  magnitude  of  the 
internal  energy  and  entropy  components  of  the  force.  Fol¬ 
lowing  Flory  et  al.  (1960),  when  the  functions  exist  and 
are  continuous,  the  order  of  a  partial  differential  does  not 
matter,  and  it  can  be  shown  that  (dS/dL)KT  =  (dfdT)KLi 
such  that 

/=  0 E/dL)VtTi„  +  T(df/dT)VtL>n.  (1.6) 

Interestingly,  this  can  be  rewritten  as 

-  T(d\nUIT]ldT)ytL,n,  (1.7) 

which  allows  for  an  experimental  estimate  of  the  fjf  ratio 
from  the  slope  of  a  plot  of  ln(//7)  versus  temperature 
under  conditions  of  constant  V,  L  and  composition,  n. 

By  means  of  an  approximate  correction  term,  the  fjf 
ratio  is  estimated  under  the  more  usual  experimental  con¬ 
ditions  of  constant  P,  L  and  at  equilibrium  (denoted  by 
subscript  eq),  e.g.  with  surrounding  solvent,  i.e. 

M=  -  T(d]n\fT\/dT)P,^ .  -  ft, T/[a\VJV )  -  1],  (1.8) 

where  /3eq  =  (31n V/dT)P>Lieq  is  the  thermal  expansion  coef¬ 
ficient;  a  is  the  fractional  increase  in  length,  L/L b  with  the 
subscript  i  indicating  initial  length;  and  Vy  and  V  are  the 
volumes  of  the  elastomer  before  and  after  elongation, 
respectively  (Hoeve  &  Flory  1962;  Dorrington  &  McCrum 
1977).  The  correction  term  was  derived  under  the 
assumption  of  random-chain  networks  with  a  Gaussian 
distribution  of  end-to-end  chain  lengths. 

As  we  will  see,  a  random-chain  network  comprised  of 
a  Gaussian  distribution  of  end-to-end  chain  lengths 

between  cross-links  is  neither  an  accurate  structural 
description  of  the  elastin  protein  nor  of  its  simpler  models. 
In  addition,  there  are  important  solvation  changes  that  can 
occur  under  experimental  conditions.  Nonetheless,  the 
insight  into  the  magnitude  of  the  entropic  component  of 
the  elastomeric  force  provided  by  equation  (1.8),  with 
experimental  estimates  of  the  quantities  in  the  correction 


term  or  when  neglecting  the  correction  term  altogether, 
is  informative,  but  only  when  the  temperature  range  for 
determination  of  the  slope  is  chosen  carefully. 

(c)  Basic  statistical  mechanical  expression  for 
entropy 

(i)  The  Boltzmann  relationship 

The  starting  point  for  the  statistical  mechanical 
expressions  of  entropy  is  the  Boltzmann  relationship, 

5  =  Plnl^,  (1.9) 

R  (1.987  cal  deg  mole-1)  is  the  gas  constant,  R  =  Nk,  N 
being  Avogadro’s  number  (6.02  x  1023  mole-1)  and  k  the 
Boltzmann  constant  (1.38  xlO-16  erg  deg  K-1),  and  Wis 
the  number  of  a  priori  equally  probable  states  accessible 
to  the  system  (Eyring  et  al.  1964).  W  is  the  volume  in 
phase  space  occupied  by  a  particular  state  of  a  molecular 
system.  In  practice,  W  becomes  the  product  of  partition 
functions  for  each  of  the  degrees  of  freedom  of  the  mol¬ 
ecular  system.  There  are  3  n  —  k  d.f.,  where  n  is  the  num¬ 
ber  of  atoms  in  the  molecule  and  k  is  the  holonomic 
constraints  on  the  system.  In  general,  there  are  three 
translational  and  three  rotational  degrees  of  freedom,  and 
the  remainder  are  vibrational  degrees  of  freedom.  For 
large-sized  molecules,  the  vibrational  degrees  of  freedom 
also  include  motions  of  rotation  (or  torsional  oscillations) 
about  bonds.  Due  to  the  constraints  of  the  elasticity 
measurement  k  =  6;  the  ends  of  the  molecules  are  effec¬ 
tively  fixed  in  space  such  that  there  are  neither  whole  mol¬ 
ecule  translational  nor  rotational  degrees  of  freedom. 

(ii)  Fundamental  expression  for  the  change  in  entropy  on 
extension 

Fortunately,  when  calculating  entropic  contributions  to 
the  elasticity,  interest  centres  on  the  change  in  entropy  on 
extension  from  a  relaxed,  ST ,  to  an  extended,  Se,  state  of 
the  elastomer, 

A 5=  (Se  -  ST )  -  R\n(Wc/VT).  (1.10) 

As  the  change  in  entropy  is  a  ratio,  it  becomes  relatively 
straightforward  to  calculate  the  contribution  to  the  change 
in  entropy  of  a  particular  expression  of  motion  accessible 
to  the  representation  of  the  molecular  structure  under 
examination  in  its  relaxed  and  extended  states.  For  our 
purposes  the  P-spiral  structure  of  the  (Gly-Val-Gly-Val- 
Pro)„,  or  poly(GVGVP),  model  of  elastin  is  used  and  two 
different  means  of  representing  that  structure  are  con¬ 
sidered.  The  first  is  the  enumeration  of  states  in  configur¬ 
ation  (<t>~  &  torsion  angle)  space  within  a  given  energy 
cut-off,  or  weighted  by  the  Boltzmann  distribution,  for 
relaxed  and  extended  states,  and  the  second  is  to  use  the 
molecular  dynamics  approach  to  determine  the  change  in 
root  mean  square  (rms)  torsional  oscillations  that  occur 
on  extension  by  the  same  amount  as  used  in  the  enumer¬ 
ation  of  states  approach.  The  magnitudes  of  the  calculated 
entropy  changes  on  extension  obtained  by  each  approach 
can  then  be  compared. 

(d)  Historical  notes  of  proposed  mechanisms  for 
the  entropic  elasticity  of  elastin 
(i)  Classical  (random  chain  network)  theory  of  rubber 
elasticity 

During  the  last  half  of  the  last  century,  identification  of 
ideal  (dominantly  entropic)  protein  elasticity  required  (by 
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Figure  1.  Plot  of  the  percentage  water  by  weight  as  a 
function  of  temperature  for  poly(GVGVP)  showing  a  stable 
intermediate  structured  state  of  63%  water  by  weight 
between  30  °C  and  60  °C,  a  state  that  slowly  denatures 
above  60  °C  to  form  a  state  containing  32%  water  by 
weight.  Graph  plotted  from  the  data  of  Urry  et  al.  (19856). 
Amphiphilic  LCST  polymers  such  as  PNIPAM,  exhibit  only 
a  single,  sharper  transition,  to  what  is  generally  considered  a 
disordered  state,  with  ca.  30%  water  by  weight.  Curve 
plotted  from  the  data  of  Grinberg  et  al.  (1999). 


the  recognized  authorities  of  the  period)  the  conclusion 
that  the  ideal  elastomer  comprised  of  random  chain  net¬ 
works  with  a  Gaussian  distribution  of  end-to-end  chain 
lengths  between  cross-links.  For  elastin,  the  protein  of 
interest  here,  this  point  of  view  began  with  a  paper  by 
Hoeve  &  Flory  (1958).  This  perspective  became 
entrenched  in  the  minds  of  the  interested  scientific  com¬ 
munity  by  the  award,  in  the  autumn  of  1974,  of  the  Nobel 
Prize  to  Paul  Flory  following  reaffirmation  of  the  random- 
chain  network  conclusion  (Hoeve  &  Flory  1974).  The 
message  of  the  paper  is  unmistakable.  In  the  synopsis 
(p.  677)  it  is  stated,  ‘A  network  of  random  chains  within 
the  elastic  fibres,  like  that  in  a  typical  rubber,  is  clearly 
indicated.*  Furthermore,  in  fig.  1  on  p.  678  of  their  paper 
there  is  a  structural  representation  of  the  chains  between 
cross-links  with  a  statement  in  the  figure  legend  that, 
‘Configurations  of  chains  between  cross-linkages  are  much 
more  tortuous  and  irregular  than  shown.* 

(ii)  Decrease  in  solvent  entropy  (due  to  hydrophobic 
hydration)  on  extension 

The  purpose  of  the  Hoeve  &  Flory  (1974)  paper,  which 
presented  no  additional  data,  was  to  refute  an  earlier  pub¬ 
lication  by  Weis-Fogh  &  Andersen  (1970),  which  had  sug¬ 
gested  an  alternative  mechanism  for  the  elasticity  of 
elastin.  The  alternative  mechanism  proposed  that 
hydration  of  hydrophobic  side  chains  of  the  protein,  which 
become  exposed  to  solvent  on  extension,  would  be 
responsible  for  the  stretch-induced  decrease  in  entropy, 


i.e.  by  this  proposal  changes  in  solvent  entropy  became  the 
source  of  the  entropic  component  of  the  elastomeric  force. 

(iii)  Damping  of  internal  chain  dynamics  on  extension 

Just  over  a  decade  later  yet  another  perspective  was  put 
forward.  Considering  structural  studies  and  molecular 
mechanics  calculations  of  the  most  prominent  repeating 
sequence  of  bovine  elastin,  (Gly-Val-Gly-Val~Pro)n, 
entropic  elasticity  was  described  as  arising  from  a  decrease 
in  available  configuration  space  on  extension  (Urry  et  al. 
1982).  Equivalently,  treating  the  experimentally  and  com¬ 
putationally  derived  regularly  repeating  structure  from  the 
perspective  of  molecular  dynamics,  the  damping  of 
internal  chain  dynamics  on  extension  described  the  same 
decrease  in  entropy  for  the  same  degree  of  extension 
(Chang  &  Urry  1989). 

In  §  4,  these  three  mechanisms  are  considered  in  more 
detail,  following  consideration  of  new  experimental  data 
on  models  of  elastin  together  with  previously  published 
relevant  data. 

(e)  Inverse  temperature  transition  behaviour  of 
elastin  and  its  models 
(i)  Increase  in  order  with  increase  in  temperature 
Filament  formation  on  raising  the  temperature 

When  the  temperature  of  an  aqueous  solution  of  the 
precursor  protein,  tropoelastin,  is  increased,  the  protein 
aggregates  to  form  a  denser  viscoelastic  phase,  called  a 
coacervate.  When  a  droplet  of  the  aqueous  suspension  of 
incipient  aggregates  is  placed  on  a  carbon-coated  grid, 
negatively  stained  with  uranyl  acetate  and  oxalic  acid  at 
the  appropriate  pH,  and  examined  under  TEM,  fibrils 
comprised  of  parallel-aligned  filaments  are  observed  with 
a  5  nm  periodicity  (Cox  et  al.  1974).  Similar  results  are 
obtained  for  a-elastin,  a  70  000  Da  fragment  of  the  elastic 
fibre  (Cox  et  al  1973),  for  high  polymers,  poly(GVGVP) 
and  poly(GVGVAP)  of  repeating  sequences  of  elastin 
(Volpin  et  al  1976).  Fibrils  of  similar  dimensions  have 
been  reported  for  fibrous  elastin  itself  using  similar  tech¬ 
niques  (Gotte  et  al.  1974). 

Crystallization  of  cyclic  analogues  on  raising  the  temperature 

Most  strikingly,  cyclic  analogues  of  repeating  sequences 
of  elastin  crystallize  on  raising  the  temperature  (Urry  et  al. 
1978;  Cook  et  al.  1980)  and  redissolve  on  lowering  the 
temperature.  Without  ambiguity,  these  polymers  increase 
order  with  increase  in  temperature  through  a  transition 
temperature  range. 

We  refer  to  this  phase  transition  as  an  inverse  temperature 
transition.  While  the  peptide  component  of  this  water- 
peptide  system  increases  in  order  on  increasing  the  tem¬ 
perature,  ordered  water  molecules  surrounding  hydro- 
phobic  residues  become  less-ordered  bulk  water  as 
hydrophobic  groups  separate  from  solution.  Accordingly, 
the  overall  change  effected  by  the  phase  transition  is 
towards  less  order,  in  keeping  with  the  second  law  of  ther¬ 
modynamics. 

Indeed,  the  phase  diagram  for  HMW  polymers  of 
repeating  sequences  of  elastin  is  inverted  with  the  soluble 
phase  below,  at  lower  temperature,  and  the  insoluble 
phase  above  the  binodal  or  coexistence  line  and  with  an 
inverted  coexistence  line,  i.e.  with  the  curvature  of  the 
coexistence  line  convex  to  the  volume  fraction  axis 
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(Sciortino  et  al.  1993 ;  Manno  et  al.  2001).  Conversely,  the 
usual  circumstance  for  polymers  is  for  the  coexistence  line 
to  be  concave  to  the  volume  fraction  axis  and  for  the  poly¬ 
mer  to  be  insoluble  below  and  soluble  above  the  coexist¬ 
ence  line  (Flory  1953). 

(f)  Composition  of  poly  (GVGVP)  in  water  as  a 
function  of  temperature 

The  temperature  dependence  of  composition  of  the 
poly(GVGVP)-water  system  provides  a  clear  visualization 
of  a  water-containing,  structured  state  at  intermediate 
temperatures  (Urry  et  al.  19856).  Figure  1  schematically 
illustrates  the  plot  of  percentage  water  by  weight  as  a  func¬ 
tion  of  temperature.  As  the  polymer  is  miscible  with  water 
in  all  proportions  below  20  °C,  the  plot  arbitrarily  starts 
at  ca.  90%  water  by  weight  near  0  °C.  A  rise  in  tempera¬ 
ture  from  20  to  30  °C  causes  a  reversible  phase  separation 
to  form  a  state  that  is  63%  water  by  weight.  The  percent¬ 
age  water  by  weight  of  this  intermediate  water-retaining 
state  remains  essentially  unchanged  until  the  temperature 
rises  above  60  °C.  Then  a  very  slow  irreversible  transition 
occurs  on  raising  the  temperature  from  60  to  80  °C  to 
form  a  state  that  is  32%  water  by  weight.  Dialysis  of  the 
polymer  against  100  000  molecular  weight  cut-off  mem¬ 
branes,  after  prolonged  heating  at  80  °C,  resulted  in  a  less 
than  1  %  loss  of  polymer.  Prolonged  heating  did  not  result 
in  chain  breakage,  and  therefore,  the  slow  irreversible 
transition  was  not  due  to  chain  breakage. 

In  figure  1,  the  results  on  poly(GVGVP)  are  compared 
to  those  of  Grinberg  et  al  (1999)  on  the  petroleum-based 
polymer,  PNIPAM.  This  polymer  shows  an  inverted 
phase  diagram  like  that  of  the  elastin-based  polymers. 
Though  with  a  composition  more  hydrophobic  than  that 
of  poly(GVGVP),  the  transition  for  this  LCST  polymer 
occurs  at  a  higher  temperature  (34  °C),  is  much  sharper, 
and  results  directly  in  the  formation  of  the  state  with  ca. 
30%  water  by  weight.  A  reasonable  description  of  this 
PNIPAM  state  of  30%  water  by  weight  is  as  a  random- 
chain  network.  It  is  clear  that  poly(GVGVP)  forms  an 
intermediate,  water-containing  state  that  holds  a  fixed 
amount  of  water  until  irreversible  denaturation  occurs  at 
high  temperature.  This  state  of  63%  water  by  weight  is  an 
elementary  property  that  simply  is  not  consistent  with  the 
recently  computed  ‘compact  amorphous  globule*  for 
(GVGVP) ,  8  attributed  to  Li  et  al.  (2000). 

(i)  A  structured  state  at  intermediate  temperatures  for 
elastin-based  systems 

Both  the  elastin  fibre  and  a-elastin  at  intermediate  tem¬ 
peratures  form  similar  water-containing  states  of  ca.  50% 
water  by  weight  (Partridge  et  al.  1955;  Partridge  &  Davis 
1955;  Partridge  1966,  1967).  It  will  also  be  demonstrated 
below  (see  figure  8)  that  both  cross-linked  poly(GVGVP) 
and  fibrous  elastin  exhibit  a  slow,  irreversible,  loss  of  elas¬ 
tic  force  due  to  prolonged  exposure  to  temperatures  above 
60  °C  when  at  fixed  extension.  These  are  further  demon¬ 
strations  of  the  slow  irreversible  denaturation  observed  in 
figure  1,  that  occurs  above  60  °C,  for  poly  (GVGVP). 
Furthermore,  detailed  dielectric  relaxation  studies  on  the 
elastin-based  systems — fibrous  elastin,  a-elastin,  and 
poly(GVGVP) — demonstrate  molecular  motions,  relax¬ 
ations,  limited  to  localized  frequency  ranges  (see  figures  6 
and  7),  and  these  could  only  be  the  result  of  non-random, 
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Figure  2.  Schematic  diagram  of  an  AFM  apparatus  adapted 
for  varying  the  ^-direction  to  obtain  force-extension  curves 
for  single  chains. 


regularly  structured,  albeit  dynamic,  conformational 
states.  These  experimental  results  are  not  consistent  with 
random-chain  networks  nor  are  they  consistent  with  com¬ 
pact  amorphous  globules. 

2.  MATERIAL  AND  METHODS 

(a)  Preparation  of  elastin  models 

In  the  native  state  of  elastin,  the  precursor,  tropoelastin,  a  pro¬ 
tein  of  about  70  000  Da,  occurs  as  the  overwhelmingly  domi¬ 
nant  component  (90+%  by  weight)  of  the  mammalian  elastic 
fibre.  In  the  vascular  wall,  for  example,  the  elastic  fibre  resides 
within  the  extracellular  matrix  as  a  structure  several  microns  in 
diameter.  No  functional  role  is  given  to  non-elastin  components 
within  the  elastic  fibre.  During  cross-linking  by  the  enzyme,  lysyl 
oxidase,  it  can  be  trapped  within  the  fibre,  and  a  fine  coating 
of  microfibres,  some  hundreds  of  nanometers  in  thickness,  can 
surround  the  fibre.  Due  to  the  cross-linking,  elastin  is  an  insol¬ 
uble  material  of  the  extracellular  matrix.  Efforts  to  isolate  it,  in 
order  to  define  its  properties,  require  Draconian  methods.  For 
example,  refluxing  in  sodium  hydroxide,  denaturants  and  pro¬ 
teolytic  enzymes  are  used  to  remove  all  of  the  surrounding 
matrix  components.  This  insoluble  product  is  called  purified 
fibrous  elastin. 

Alternatively,  for  the  study  of  elastin,  model  approaches  have 
been  used  such  as  isolating  and  characterizing  (i)  the  precursor 
protein,  tropoelastin,  which  may  be  obtained  in  ill-defined  inter¬ 
mediate  states  of  cross-linking,  (ii)  the  pure  precursor  protein 
prepared  by  microbial  biosynthesis,  (iii)  chemical  degradation 
products  of  purified  fibrous  elastin  which  have  been  prepared, 
and  (iv)  chemically  and  microbially  synthesized  repeating 
sequences  found  within  the  protein  sequence  as  well  as  interest¬ 
ing  analogues  which  have  been  prepared. 

(i)  Synthesis  of  model  systems 
Chemical  synthesis 

Initially,  several  thousand  protein-based  polymers  were  syn¬ 
thesized  in  order  to  determine  the  conformation  and  function 
of  these  model  proteins.  This  involved  the  synthesis  of  oligomers 
and  high  polymers  of  repeating  tetra-,  penta-,  hexa-,  nona-  and 
even  decapeptides  reported  in  the  sequence  of  elastin.  By  means 
of  chemical  synthesis,  many  compositions  could  be  prepared 
and  characterized,  and  axioms  were  developed  for  the  function 
of  these  elastic  model  proteins  in  energy  conversion  and  for  their 
use  as  biomaterials  for  medical  and  non-medical  applications. 
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Figure  3.  Single-chain  force-extension  curves  for  Cys- 
(GVGVP)rtx25i-Cys  at  a  temperature  below  the  onset 
temperature  for  hydrophobic  folding  and  assembly  for  this 
composition.  The  initial  trace  appears  at  the  lowest  position 
on  the  graph  and  subsequent  traces  of  the  same  chain, 
without  intervening  detachment,  are  displaced  250  pN 
vertically.  The  second  and  fifth  curves  exhibit  perfect 
reversibility.  The  curves  marked  by  an  asterisk,  and  the 
initial  curve,  were  the  result  of  a  period  of  at  least  30  s  in 
the  relaxed  state,  where  one  possibility  is  that  some 
hydrophobic  back  folding  occurs. 

Microbial  biosynthesis 

Once  a  specific  composition  was  identified  as  being  of  suf¬ 
ficient  interest,  it  was  prepared  by  means  of  recombinant  DNA 
technology.  With  this  approach,  hundreds  of  grams  could  be 
prepared  in  a  single  fermentation.  On  the  other  hand,  chemical 
synthesis  at  the  laboratory  bench  would  require  much  of  a  year 
to  prepare  such  a  quantity. 

(ii)  Preparation  of  natural  materials 
Isolation  of  the  precursor  protein ,  tropoelastin 

Initially,  the  precursor  protein  was  prepared  from  animals 
rendered  copper  deficient  and  treated  with  inhibitors  of  enzy¬ 
matic  cross-linking.  The  product  was  meagre  in  amount  and  a 
mixture  of  partially  cross-linked  and  oxidized  lysyl  side  chains.  It 
was  nonetheless  helpful  in  determining  many  of  the  fundamental 
physical  properties  relating  to  structure  and  elastogenesis. 

Preparation  of  tropoelastin  by  microbial  biosynthesis 

In  an  important  development,  Weiss  and  co-workers  utilized 
recombinant  DNA  technology  to  prepare  the  precursor  protein, 
tropoelastin,  in  Escherischia  coliy  which  enabled  the  preparation 
of  pure  protein  in  large  quantities  (Vrhovski  et  al  1997;  Wu  et 
al  1999).  Important  aspects  deduced  by  the  model  systems  and 
by  preparations  of  model  structures  derived  from  the  elastin  and 
from  sequence  knowledge  have  been  confirmed  and  extended 
by  recombinantly  prepared  tropoelastin. 

Preparation  of  a-elastin 

A  model  system  useful  in  stepping  systematically  from  poly¬ 
mers  of  model-repeating  sequences  to  the  natural  fibrous  elastin 
has  been  a  chemical  fragmentation  product  from  fibrous  elastin 
called  a-elastin,  prepared  and  chemically  characterized  by  Par¬ 
tridge  and  co-workers  (Partridge  et  al.  1955;  Partridge  &  Davis 
1955;  Partridge  1966,  1967).  It  is  a  multiply  cross-linked, 
70  000  molecular  weight  fragment  comprising  16  chains. 


(iii)  Purification  using  phase  transitional  behaviour 
In  all  cases  these  protein-based  polymers  are  soluble  in  water 
at  a  lower  temperature  and  phase  separate  on  raising  the  tem¬ 
perature  above  that  for  the  onset  of  a  hydrophobic  folding  and 
assembly  transition.  This  phase  separation  process  has  been 
termed  coacervation.  It  is  a  fundamental  property  of  these 
polymers  having,  as  they  do,  the  proper  balance  of  apolar 
(hydrophobic)  and  more  polar  residues. 

(b)  Cross-linking  of  elastin  models 

(i)  General  cross-linking  procedure 

Cobalt  60  7-irradiation  was  used  to  cross-link  the  phase-separ¬ 
ated,  coacervate  state  with  its  interpenetrating  polymer  chains. 

(ii)  Efforts  to  cross-link  a-elastin 

Even  though  a-elastin  undergoes  phase  separation,  it  does  not 
form  an  elastic  matrix  on  y-irradiation.  As  a  cluster  of  some  16 
cross-linked  chains,  with  a  total  molecular  weight  of  70  kDa, 
a-elastin  does  not  cross-link  by  y-irradiation  because  the  cluster 
of  chains  of  one  a-elastin  molecule  does  not  interpenetrate  suf¬ 
ficiently  with  surrounding  molecules.  This  is  in  contrast  with 
similar  molecular  weights  of  poly(GVGVP)  which  readily  cross¬ 
link  by  y-irradiation  when  in  the  intermediate  structured  state 
between  30  and  60  °C. 

(iii)  Efforts  to  cross-link  heat  denatured  poly(GVGVP) 

While  the  phase-separated  state  of  poly(GVGVP)  with  63% 

water  by  weight  (figure  1)  cross-links  very  effectively,  the  heat- 
denatured  state  obtained  on  prolonged  heating  above  60  °C 
does  not.  Even  though  the  density  of  chains  is  twice  as  great  for 
the  denatured  state  formed  above  60  °C  (a  state  of  32%  water 
by  weight  as  seen  in  figure  1),  the  chains  no  longer  interpen¬ 
etrate  and  associate  sufficiently  to  result  in  an  elastic  matrix  on 
y-irradiation.  This  is  simply  another  demonstration  that  the 
intermediate  state,  between  30  and  60  °C,  is  a  structured  state 
of  interpenetrating  filaments  rather  than  being  globular  struc¬ 
tures  with  limited  contact  between  units,  as  occurs  with  a-elastin 
or  as  would  be  the  case  for  amorphous  compact  globules. 

3.  SPECIALIZED  METHODOLOGIES  AND  ANALYSIS 
OF  RESULTS 

(a)  Single  molecule  force-extension  curves  by  AFM 

(i)  Preparation  of  the  sample  for  AFM 

Cleaned  glass  microscope  slides  were  coated  with  ca. 
30  nm  of  gold  in  a  home-built  evaporation  chamber.  Si3N4 
AFM  tips  were  used  (Microlevers,  Park  Scientific  Instru¬ 
ments,  Sunnyvale,  CA),  also  coated  with  gold.  For  the 
preparation  of  mixed  self-assembled  monolayers,  quan¬ 
tities  of  1  mg  of  the  polypentapep tides  and  0.5  mg  of 
methoxy-PEG-thiol  (Mw  5000)  were  dissolved  in  1  ml 
Milli-Q  water  unless  otherwise  noted.  A  20  pi  portion  of 
this  solution  was  incubated  on  the  gold-coated  slide  for 
about  30  min  at  3  °C  and  then  rinsed  with  Milli-Q  water. 

(ii)  The  AFM  instrument 

A  detailed  description  of  the  AFM  force  spectroscopic 
technique  and  details  of  the  instrumental  set-up  have  been 
supplied  elsewhere  (Oesterhelt  et  al.  1999;  Clausen- 
Schaumann  et  al.  2000).  Briefly,  the  tip  of  a  cantilever  is 
brought  into  contact  with  molecules  on  the  surface  and 
then  retracted,  its  deflection,  and  therefore  the  force,  are 
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detected  with  a  laser  by  optical  lever  detection  (Cleveland 
et  al.  1993),  while  the  ^-distance  is  controlled  by  a  strain 
gauge  (figure  2).  The  nominal  spring  constants  of  cantile¬ 
vers  used  in  the  experiments  were  10  mN  m-1.  Before  the 
first  approach  of  the  AFM  tip  to  the  surface,  the  spring 
constants  of  each  lever  were  individually  calibrated  by 
measuring  the  amplitude  of  their  thermal  oscillations 
(Butt  &  Jaschke  1995).  The  sensitivity  of  the  optical  lever 
detection  was  measured  by  indenting  the  AFM  tip  into  a 
hard  surface.  All  experiments  were  conducted  in  Milli-Q 
water  at  room  temperature  (21  °C),  unless  otherwise 
specified. 

(iii)  Obtaining  the  single-chain  force-extension  curve 

For  the  characterization  of  elastic  properties,  the  gold- 
coated  AFM  tip  was  brought  into  contact  with  the  poly- 
pentapeptides  on  the  microscope  slide  by  manual  control. 
The  AFM  tip  and  the  polymer  layer  were  kept  in  contact 
under  a  contact  force  of  several  nanoNewtons  for  ca.  30  s 
to  allow  for  an  attachment  of  polypeptide  chains  to  the  tip. 
Usually,  one  or  several  chains  adsorbed  non-covalently  at 
some  position  on  the  tip,  possibly  a  chemical  reaction 
between  the  gold-coated  cantilever  and  the  cysteine  on  the 
end  of  a  polypentapeptide  sometimes  allowed  for  very  high 
rupture  forces.  Upon  retraction  of  the  cantilever,  individ¬ 
ual  polypentapeptides  were  stretched  between  the  surface 
and  the  AFM  tip. 

Note  that  the  first  force-distance  profile  recorded  after 
tip-substrate  contact  can  be  rather  complex,  consisting  of 
contributions  from  stretching  several  polypentapeptides, 
desorption  from  the  substrate  and/or  cantilever,  bond  rup¬ 
ture  of  short  strands,  as  well  as  interchain  aggregation  and 
entanglements.  Therefore,  in  each  measurement  the  can¬ 
tilever  was  first  retracted  from  the  substrate  to  a  distance 
at  which  unspecific  adhesion  was  no  longer  observed. 
Then,  in  successive  retraction-approaching  cycles  the  dis¬ 
tance  range  was  increased  while  trying  to  avoid  contact 
between  the  tip  and  additional  polypentapeptide  strands 
at  the  substrate  surface,  until  only  one  polymer  strand 
remained  between  the  tip  and  the  substrate.  The  force- 
distance  profile  of  this  strand  was  then  measured  repeat¬ 
edly  until  rupture. 

(iv)  Analysis  of  the  results 

The  experimental  traces  were  fitted  by  an  extended 
WLC  model  (Bustamante  et  al.  1994;  Odijk  1995)  includ¬ 
ing  linear  elastic  contributions  arising  from  the  stretching 
of  bond  angles  and  covalent  bonds: 

_ —  4. _ l _ i  (2  1) 

kBT  L  Ko  4(1  -RJL+FIKo)2  4* 

In  this  expression,  Rz  is  the  measured  end-to-end-dis- 
tance  at  any  given  force,  i7,  and  L  is  the  contour  length 
of  the  stretched  chain  (polypentapeptide)  under  zero  force 
(F=0).  It  should  be  appreciated  that  the  WLC  model 
entirely  neglects  any  discrete  molecular  structure  along  the 
chain,  and  describes  the  polymer  as  a  continuous  rod  of 
constant  bending  module.  The  characteristic  length  scale 
expressing  the  polypeptide’s  bending  rigidity  is  the  persist¬ 
ence  length,  Lp,  which  is  defined  as  the  decay  length  of 
the  directional  correlation  along  the  polymer  chain. 
Finally,  the  chain’s  extensibility  upon  stretching  is 
described  by  the  segment  elasticity,  J5^,  which  is  intro¬ 


duced  into  equation  (2.1)  as  a  linear  term.  Hereby,  K 0  can 
be  understood  as  the  inverse  of  the  normalized  compliance 
of  a  Hookean  spring;  the  spring  constant  of  the  polymer 
chain  is  given  by  KJL . 

This  resulted  in  values  for  the  persistence  length  Lp  of 
(GVGVPJnxasi:  LP  =  0.4  nm  in  the  low  force  regime  and 
Lp  =  0.6  nm  when  fitted  to  higher  forces,  which  are  com¬ 
parable  to  values  measured  previously  for  other  polypep¬ 
tide  backbones,  e.g.  of  proteins  (Rief  et  al .  1997).  The 
reversible  traces  provide  evidence  for  single-chain  ideal 
elasticity,  i.e.  any  molecular  processes  related  to  the  elong¬ 
ation  and  relaxation  of  the  polypeptide  chains  must  be  fast 
on  the  time-scale  of  the  AFM  experiment.  Apart  from  the 
entropic  contributions  at  low  forces  below  ca.  20  pN, 
which  are  related  to  the  unfolding  of  the  polymer  coil 
(‘classical  conformational  entropy’),  it  has  been  previously 
shown  that  ideal  elastic  behaviour  in  the  high  force  regime 
may  also  result  from  enthalpic  bond-angle  deformations 
and  configurational  changes  along  the  polymer  backbone 
(Rief  et  al.  1997).  As  will  be  shown  below,  additional 
entropic,  and  therefore  reversible  elasticity,  in  the  high 
forces  regime  can  also  arise  from  backbone  torsional 
movements  (rocking),  wherein  the  contribution  to  the 
entropic  component  of  the  force  increases  with  a  decrease 
in  the  frequency  of  the  oscillator.  Temperature-dependent 
AFM  measurements  should  allow  for  a  proof  of  this  find¬ 
ing  at  the  single  molecule  level. 

The  small  deviations  from  perfect  reversibility  in  traces 
3,  4  and  6  of  figure  3  (marked  by  the  asterisk)  occurred 
when  the  chain  was  held  for  at  least  30  s  in  the  relaxed 
state.  These  deviations  could  have  several  origins.  There 
could  be  a  configurational  transition,  e.g.  an  opening  of 
p-tums  and/or  hydrophobic  interaction,  a  back  folding, 
between  different  parts  of  the  chain.  The  latter  is  favoured 
by  the  observed  rate  dependence  of  this  non-equilibrium 
contribution. 

The  persistence  length  LP  of  (GVGIP)mx260  is  difficult 
to  determine,  when  intra-  and/or  interchain  aggregation 
gives  rise  to  a  strong  hysteresis  even  when  located  far  away 
from  the  surface  within  two  successive  cycles  (figure  4). 
It  is  determined  to  exhibit  an  Lp-  0.7  nm,  which  might 
be  slightly  higher  than  with  (GVGVP)roc25i3  but  this  value 
lies  well  within  the  error  range  of  ca.  20%.  This  strong 
aggregation  of  (GVGIP)*^  occurs  at  the  same  concen¬ 
tration,  at  which  the  aggregation  for  (GVGVP)„x251  is 
much  smaller.  This  reflects  the  increase  in  hydrophobicity 
due  to  an  additional  CH2  group  in  GVGIP,  which  is  the 
only  difference  with  GVGVP.  Unfortunately,  it  is  not  yet 
clear  whether  the  aggregation  is  mainly  due  to  back  fold¬ 
ing  within  a  single  chain  or  due  to  association  between 
different  chains.  This  is  under  investigation. 

(v)  Comparison  of  AFM  single-chain  and  macroscopic 
elastic  moduli 

For  the  characterization  of  elastomers  at  the  macroscopic 
scale,  force  versus  elongation  is  plotted,  but  instead  of  being 
represented  as  an  actual  increase  in  length,  as  in  figures  3 
and  4,  elongation  is  stated  in  terms  of  the  ratio  of  the  initial 
length  to  the  extended  length.  Furthermore,  the  defor¬ 
mation  is  generally  uniform  so  that  the  initial  cross-sectional 
area,  As  is  used.  Therefore,  the  Young’s  elastic  modulus,  Y, 
is  the  force,  /,  times  the  initial  length,  L^  divided  by  the 
elongation,  L,  and  the  cross-section  of  the  elastomer,  that 
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is,  Y=fLJAL.  The  units  are  Nm"2  or  dynes  cm-2  where 
Nm"2  =  10  dynes  cm-2  and  N  stands  for  Newtons.  The 
elastic  modulus  of  20Mrad  ^-irradiation  cross-linked 
(GVGVP)251,  i.e.  X20  (GVGVP)251,  is  1.6  x  105Nm-2. 

In  order  to  utilize  the  data  in  figure  3  to  obtain  an  esti¬ 
mate  of  the  single-chain  elastic  modulus,  estimates  of  the 
L  JL  ratio  and  of  A  are  required.  Initially,  the  total  number 
of  repeats  that  are  being  extended  must  be  deduced.  In 
figure  3,  the  WLC  fit  of  the  force-distance  trace  shown 
gives  an  end-to-end-distance  of  L  =  733  nm,  which  we 
assume  is  the  fully  extended  chain.  This  is  too  long  for  a 
chain  of  251  pentamers,  but  appears  to  be  too  short  for  a 
chain  of  753  pentamers.  By  elimination  of  chain  lengths 
with  251  and  753  pentamers,  this  indicates  that  the  struc¬ 
ture  would  be  (GVGVP)2X25i  giving  502  pentamers,  i.e. 
2510  residues  with  a  value  of  0.29  nm  per  residue.  This 
is  shorter  than  the  usual  value  of  0.35  nm  per  residue, 
which  could  be  due  to  the  presence  of  a  proline  residue 
in  every  pentamer.  This  results  from  the  cis  Val-Pro  pep¬ 
tide  bond,  which  introduces  a  kink  in  the  extended  back¬ 
bone  for  every  pentamer.  It  should  be  noted,  however, 
that  if  all  of  the  p-tums  were  retained  at  rupture,  the 
3  x  251  chain  length  could  apply  with  ca.  1  nm  per  penta¬ 
mer,  or  just  greater  than  700  nm.  This  scenario  was  not 
considered  because  the  forces  for  detachment  can  be 
greater  than  500  pN,  and  retention  of  the  p-tum  was  con¬ 
sidered  less  likely  at  such  forces. 

It  should  also  be  noted  that,  in  general,  one  cannot 
exclude  a  priori  the  possibility  of  polypeptides  being  picked 
up  at  any  position  along  the  chain  due  to  strong  unspecific 
adsorption.  In  the  case  of  (GVGVP),^!,  however,  the 
simplicity  and  uniformity  of  the  repeating  sequence  with¬ 
out  the  presence  of  functional  groups  except,  at  the  very 
ends  of  the  chain,  are  expected  to  limit  strong  unspecific 
adsorption.  Because  the  polymer  length  estimated  from 
the  WLC  fit  is  sufficiently  close  to  what  would  be  expected 
for  2x251  proline-containing  pentamers  it  warrants,  in 
our  view,  proceeding  with  this  assumption. 

Perhaps  the  most  significant  factor  in  obtaining  a  rel¬ 
evant  value  for  the  single-chain  elastic  modulus  for  com¬ 
parison  to  the  macroscopic  case  is  the  choice  of  the 
structure  before  elongation  in  order  to  provide  an  estimate 
of  an  initial  length  and  cross-sectional  area.  In  this  regard 
two  quite  different  initial  lengths  may  be  considered,  one 
for  the  random  coil  state  below  the  hydrophobic  folding 
transition  and  another  for  the  P-spiral  structure  that  could 
only  occur  above  the  onset  temperature  for  the  hydro- 
phobic  folding  transition. 

Calculation  of  elastic  modulus  from  AFM  data  assuming  a 
random  coil  state  for  (GVGVP)nx25} 

At  temperatures  lower  than  the  onset  temperature  for 
the  hydrophobic  folding  and  assembly  transition,  the  more 
correct  representation  of  the  unstretched  single  chain 
would  be  the  random  coil  state.  Based  on  photon  corre¬ 
lation  spectroscopic  studies,  the  hydrodynamic  diameters 
of  disordered  polymeric  coils  resulting  from  502  penta¬ 
mers  (some  200  000  Da  molecular  weight)  below  the 
transition  would  be  ca.  18  nm  (San  Biagio  et  al.  1988).  By 
way  of  example,  at  a  distance  of  378  nm  in  figure  3,  the 
UJL  ratio  would  be  1  : 20.  In  this  distance  range  (i.e. 
between  300  and  500  nm),  the  slope  of  the  force  curve  is 
evaluated  as  0.0810.03  pNnm-1.  With  a  slope  of  0.08  ± 


Figure  4.  Single-chain  force-extension  curves  for  Cys- 
(GVGIP)mt260-Cys  at  a  temperature  above  the  onset 
temperature  for  hydrophobic  folding  and  assembly  for  this 
composition.  In  each  case  the  energy  expended  in 
deformation  is  greater  than  that  recovered  on  relaxation.  In 
addition,  there  appears  to  be  a  complex  relaxation  curve. 


0.03  pN  nm-1  from  figure  3,  at  a  distance  of  378  nm  from 
the  substrate,  the  increase  in  force  would  be 
28.81 10.8  pN;  i.e.  upon  elongation  from  18  to  378  nm, 
the  relative  increase  in  length,  A L/Lb  would  be  a  factor  of 
20,  so  that  a  force  of  1.441 0.54  pN  is  estimated  for  a 
100%  extension  of  the  18nm  diameter  polymeric  coil. 

Next,  this  quantity  should  be  divided  by  an  appropriate 
cross-sectional  area,  that  is,  rr  multiplied  by  the  hydrody¬ 
namic  radius,  squared,  of  the  random  coil  prior  to  exten¬ 
sion  of  the  single  chain,  i.e.  n r2  =  tt(9  nm)2  -  254  nm2. 
Given  the  foregoing  assumptions,  the  calculated  single¬ 
chain  elastic  modulus  would  be  5.7(1 2.1)  x  103  N  m"2. 
As  might  have  been  expected,  without  taking  polymer 
density  and  cross-link  density  into  consideration  (Urry  et 
al  1984),  this  value  is  a  factor  of  30  from  the  experimental 
value  at  37  °C  for  X20-(GVGVP)25l  of  1.6  x  105Nm~2, 
but  more  closely  approximates  what  is  expected  below  the 
transition  temperature  or  when  X20-(GVGVP)25t  has  been 
thermally  denatured,  as  seen  in  figure  8 b.  Concerning 
polymer  density,  it  may  be  noted  that  the  density  of  the 
low-temperature  random-coil  state  is  about  one-quarter 
that  of  die  phase-separated  state. 

Calculation  of  elastic  modulus  from  AFM  data  assuming  the 
fi-spiral  structure  for  (GVGVP)mi2si 

For  comparison,  it  is  useful  to  consider  the  structure 
representative  of  the  state  above  the  temperature  of  the 
inverse  temperature  for  hydrophobic  folding  and 
assembly.  This  addresses  the  importance  of  choosing 
structure  when  converting  the  force-extension  curve  into 
a  macroscopic  elastic  modulus.  If  the  single  (3-spiral  struc¬ 
ture  were  assumed  before  extension  with  three  pentamers 
per  turn  and  a  translation  along  the  axis  of  1  nm  per  turn, 
the  length  of  a  502  pentamer  would  be  ca.  170  nm.  This 
means  that  an  elongation  of  100%  would  occur  at  340  nm, 
and  the  change  in  force  for  extension  of  the  (3-spiral  from 
170  to  340  nm  would  be  170  nm  x  0.08  ±  0.03  pN  nm-1 
as  13.6  ±  5.1  pN  per  (3-spiral. 

The  next  issue  is  an  estimate  of  the  number  of  p-spirals 
per  m2  for  comparison  with  the  cross-linked  matrix  from 
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(a) 


(b) 


Figure  5.  (a)  Block  diagram  of  the  instrumentation  used  to  determine  the  absorption  per  unit  volume  of  elastic  samples  in  the 
acoustic  range.  (6)  Example  of  the  use  of  two  samples  that  differ  in  length  by  a  factor  of  2  as  required  to  permit  the  analysis 
of  Sheiba  (1996)  to  be  utilized. 


which  the  experimental  macroscopic  value  of 

1.6  x  105Nm"2  is  obtained.  For  a  macroscopic  state  of 
37%  water  by  weight  (see  figure  1),  an  assumed  density 
of  1.33  g  cm-2  for  the  polymer,  and  a  cross-sectional  area 
for  the  p-spiral  of  3.2  nm2,  the  number  would  be  ca . 

8.7  x  1016  p-spirals  m-2.  The  single-chain  elastic  modu¬ 
lus,  calculated  on  the  basis  of  the  p-spiral  structure,  would 
be  (13.6  ±  5.1  x  10“12Nper  p-spiral)  x  (-8.7  x  1016P-spi- 
rals  m-2)  =  1.2  x  106  N  m"2. 

Thus,  by  assuming  a  different  structure,  instead  of 
being  low  by  a  factor  of  30  as  was  the  case  above  for  a 
random-coil  structure,  which  occurs  below  the  transition 
temperature,  the  assumption  of  the  p-spiral  structure 
resulted  in  a  value  for  the  single-chain  elastic  modulus  that 
is  seven  and  a  half  times  larger  than  the  macroscopic 
value.  We  note  that  the  orientation  of  a  random  coil  is 
not  relevant  in  the  argument  relating  single-chain  data  to 
macroscopic  data.  For  the  y-irradiation  cross-linked 
macroscopic  matrix,  however,  the  P-spirals  would  be  ran¬ 
domly  orientated  with  respect  to  the  direction  of  exten¬ 
sion.  This  would  decrease  the  contribution  of  each  p-spiral 
not  aligned  in  the  direction  of  extension.  It  will  now  be 
particularly  interesting  to  obtain  the  elastic  modulus  for 
the  twisted  filament  structure  depicted  in  figure  9/  Work 
is  in  progress  for  the  preparation  of  multistranded 
twisted  filaments. 


(b)  The  acoustic  absorption  experiment 
(i)  Sample  preparation  and  experimental  set-up 

The  polymer  is  prepared  in  the  phase-separated  state, 

1. e.  the  intermediate  structured  state  with  the  interpen¬ 
etrating  chains  depicted  in  figure  1,  and  20Mrad  y- 
irradiation  cross-linked  in  the  form  of  the  desired  cylindri¬ 
cal  shape  and  with  path  lengths  that  differ  by  a  factor  of 

2,  as  shown  in  figure  5b. 

The  Naval  Research  Laboratory  at  the  Spatial  Warfare 
Systems  Center — San  Diego  (SSC — SD)  has  the  capa¬ 
bility  to  measure  accurately  the  elastic  properties  of 
materials  by  using  a  broadband  technique  developed  by 
Sheiba  (1996).  Older  techniques  are  more  cumbersome 
and  less  reliable.  The  frequency  of  measurement  ranges 
from  subsonic  to  ca.  1 5  kHz.  The  measurements  are  per¬ 
formed  in  air  inside  an  environmental  chamber.  Figure  5 
also  shows  the  block  diagram  of  the  instrumental  set-up 
for  these  measurements. 

The  network  analyser  sends  sweeping  signals  to  the 
shaker.  The  transducer  (sometimes  called  an 
accelerometer)  and  impedance  head  record  the  velocities 
at  the  base  and  at  the  top  of  the  sample.  After  recording 
data  during  this  frequency  sweep,  a  new  set  of  data  is  then 
obtained  at  a  new  temperature  setting.  Commonly  tem¬ 
peratures  can  range  from  -30°  to  70  °C.  The  velocity  rec¬ 
ord  constitutes  the  transfer  matrix.  By  mathematically 
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frequency  (MHz) 

Figure  6.  (a)  Plots  of  the  loss  factor  (absorption  per  unit 
volume)  for  a  series  of  elastomers.  Crossed  squares  represent 
the  curve  for  natural  rubber,  and  filled  diamonds  that  of  the 
structureless  curve  for  polyurethane,  which  is  recognized  as  a 
rubber  with  a  high  loss  factor  in  the  acoustic  absorption 
range.  The  upper  four  curves  are  for  20  Mrad  cross-linked 
(GVGVP)26o  (open  squares,  20  °C;  open  diamonds,  10  °C; 
open  circles,  0  °C;  open  triangles,  -5  DC).  As  the 
temperature  rises  above  the  temperature  for  the  inverse 
temperature  transition  of  hydrophobic  folding  and  assembly 
for  this  composition,  the  intensity  for  a  localized  relaxation 
grows  dramatically.  To  the  best  of  our  knowledge,  these 
model  elastic  protein-based  polymers  exhibit  the  greatest 
sound  absorption  for  the  100-1000  Hz  frequency  range. 

This  demonstrates  the  formation  of  a  regular  structure  for 
(GVGIP)„  in  contrast  to  the  lack  of  a  regular  structure  for 
the  classical  rubbers.  ( b )  The  imaginary  (absorptive)  section 
of  the  dielectric  permittivity  of  poly(GVGIP).  The  same 
phenomenon  of  increased  absorption  on  passing  through  the 
temperature  range  of  the  inverse  temperature  transition  is 
observed  in  the  1-100  MHz  range  as  exhibited  by  the  same 
composition  in  the  acoustic  absorption  range.  Reproduced, 
with  permission,  from  Buchet  et  al  (1988).  Such  low 
frequency  motions  provide  a  source  of  entropic  elasticity  and 
low  energy  for  the  structural  state  of  the  elastic  protein- 
based  polymers. 


Figure  7.  (a)  The  real  pan  of  the  dielectric  permittivity  of  a- 
elastin  showing  increasing  relaxation  intensity  at  ca.  5  MHz 
with  increasing  temperature  from  below  to  above  the 
temperature  range  for  the  inverse  temperature  transition. 

The  data  for  a-elastin  are  nearly  identical  to  that  found  for 
the  polypentapeptide  of  elastin,  poly(GVGVP)  shown  in  the 
inset  for  a  single  temperature  below  (20°)  and  one  above 
(40°)  that  of  the  inverse  temperature  transition.  The 
temperatures  for  the  a-elastin  data  on  the  left  edge  of  the 
curves  listed  as  (a)-(m)  are  68,  12.5,  20,  25,  30,  35,  40,  45, 
50,  55,  60,  65  and  70  °C,  respectively.  This  demonstrates 
that  the  same  relaxation  occurs  in  a  fragment  of  the  natural 
elastic  fibre  as  has  been  so  thoroughly  characterized  in  the 
model  elastic  protein,  poly(GVGVP).  Reproduced,  with 
permission,  from  Urry  et  al  (1985a). 

(b)  The  imaginary  (absorptive)  part  of  the  dielectric 
permittivity  of  natural  fibrous  elastin  showing  an  increasing 
relaxation  at  ca.  5  MHz  with  increasing  temperature  from 
below  (23  °C)  to  above  (58  °C)  the  temperature  range  for 
the  inverse  temperature  transition.  The  same  relaxation  is 
seen  at  58  °C  with  a  somewhat  higher  intensity  and  a  slightly 
higher  frequency  is  found  for  the  20  Mrad  cross-linked 
poly(GVGVP).  It  seems  quite  evident  that  dynamic 
structured  states  occur  in  the  natural  elastic  fibre  as 
demonstrated  for  the  model  elastic  protein-based  polymer. 
Adapted,  with  permission,  from  Luan  et  al  (1988). 
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inverting  the  transfer  matrix,  it  is  possible  to  extract  the 
complex  modulus  and  Poisson  ratio.  These  data  provide 
a  complete  description  of  a  material’s  elastic  properties. 
Further  analysis  yields  the  shear  modulus,  (i ,  the  Poisson 
ratio,  vy  and  the  loss  factor,  tj,  each  as  a  function  of  fre¬ 
quency  and  temperature.  The  loss  factor,  17,  is  the  main 
parameter  of  interest.  It  represents  the  energy  absorption 
per  unit  volume  at  a  given  frequency.  The  loss  factor  ver¬ 
sus  frequency  and  temperature  were  determined  for  sev¬ 
eral  different  elastic  protein-based  polymer  compositions. 
The  data  for  (GVGIP)„x26o  are  given  in  figure  6 a. 
Mathcad  7.0  was  used  to  perform  the  mathematical 
analysis. 

(ii)  Comparison  of  the  acoustic  absorptions  of  (GVGIP)  and 
natural  rubber 

Natural  rubber  has  been  well  characterized  as  a  ran¬ 
dom-chain  network.  This  means  that  there  are  no  regular 
repeating  conformational  features  in  the  polymer  chain. 
The  consequence  of  a  random  chain  is  a  broad,  slowly 
changing  dependence  on  frequency  of  physical  properties, 
such  as  absorption.  Simply  stated,  the  absence  of  structure 
means  that  the  barrier  to  rotation  about  each  bond  will 
tend  to  be  different.  This  is  what  is  observed  in  figure  6 a 
for  natural  rubber  and  a  synthetic  petroleum-based  poly¬ 
mer,  polyurethane,  chosen  for  its  relatively  high  loss  factor 
in  the  frequency  range  reported. 

At  frequencies  lower  than  those  of  the  IR  spectroscopic 
range,  however,  any  material  that  exhibits  a  localized 
absorption  band  as  a  function  of  frequency  must  do  so  by 
virtue  of  a  regular  repeating  structure.  Interestingly,  this 
is  seen  in  Figure  6 a  for  X20-(GVGIP)26o  in  the  range 
200  Hz  to  7  kHz,  and  the  intensity  of  the  absorption  per 
unit  volume  increases  as  the  temperature  is  raised  through 
the  temperature  interval  of  the  inverse  temperature  tran¬ 
sition.  Acoustic  absorption  increases  as  the  structure 
becomes  more  regular!  Besides  demonstrating  the  remark¬ 
able  acoustic  absorption  property  of  this  material,  the 
acoustic  absorption  data  provide  an  additional  demon¬ 
stration  of  the  existence  of  a  regular  structure  for 
poly(GVGIP)  and  for  related  polymers.  Actually,  this  is 
simply  further  confirmation  of  earlier  dielectric  relaxation 
data,  which  demonstrated  the  development  of  an  intense 
localized  relaxation  at  ca.  5  MHz  as  the  temperature  was 
raised  through  the  interval  of  the  inverse  temperature  tran¬ 
sition,  as  reviewed  below. 

(c)  Previous  dielectric  relaxation  studies  on 
elastin-related  systems 

(i)  Dielectric  relaxation  and  acoustic  absorption  data  on  the 
same  polymer 

Figure  6b  is  a  three-dimensional  plot  of  the  imaginary 
(absorption)  part  of  the  dielectric  permittivity  of 
poly(GVGIP)  as  a  function  of  frequency  and  temperature. 
Just  as  for  the  acoustic  absorption  study  in  figure  6 ay  there 
is  an  increase  in  the  absorption  (imaginary)  component  of 
the  dielectric  relaxation  as  the  temperature  rises  through 
the  temperature  interval  of  the  inverse  temperature  tran¬ 
sition  for  hydrophobic  folding  and  assembly.  Thus,  the 
data  on  the  elastin  models  can  now  be  directly  compared 
to  a-elastin  and  fibrous  elastin,  itself.  Mechanical  reson¬ 
ances  due  to  regularly  repeating  structures  have  now  been 
observed  both  in  the  vicinities  of  1  kHz  and  5  MHz. 


(ii)  Comparison  of  poly  (GVGVP)  and  a  -elastin 

Figure  la  is  a  graph  of  the  real  part  of  the  dielectric 
permittivity  as  a  function  of  temperature  for  a-elastin  and 
for  poly(GVGVP)  in  the  inset.  In  both  cases,  there  is 
development  of  an  intense  relaxation  at  ca.  5  MHz  as  the 
protein-based  elastomers  hydrophobically  fold  and 
assemble  on  increasing  the  temperature  through  the  range 
of  the  inverse  temperature  transition.  Both  elastomeric 
molecular  systems  exhibit  the  same  development  of  reg¬ 
ular  structure  on  raising  the  temperature  through  the 
range  of  the  inverse  temperature  transition. 

(iii)  Comparison  of  cross-linked  poly  (GVGVP)  to  fibrous 
elastin 

Figure  lb  is  a  graph  of  the  imaginary  (absorption)  part 
of  the  dielectric  relaxation  for  purified  fibrous  elastin 
(from  bovine  ligamentum  nuchae)  as  a  function  of  tem¬ 
perature  from  23  °C  to  58  °C.  Figure  lb  also  contains  a 
plot  of  the  same  data  at  58  °C  for  the  elastomer,  20  Mrad 
cross-linked  poly  (GVGVP).  While  there  is  a  small  increase 
in  mean  frequency  and  intensity  for  the  X20-(GVGVP)25i 
relaxation,  it  is  certain  that  there  is  hydrophobic  folding 
and  assembly  resulting  in  the  formation  of  regular  struc¬ 
tures  in  the  natural  fibrous  elastin  that  are  similar  to  those 
of  cross-linked  poly  (GVGVP).  Incidentally,  the  decrease 
in  intensity  at  1  GHz  for  fibrous  elastin  reflects  the  loss  of 
hydrophobic  hydration  as  the  hydrophobic  folding  pro¬ 
ceeds.  This  is  due  to  a  dielectric  relaxation  at  ca.  5  GHz, 
which  results  from  hydrophobic  hydration  (Urry  et  al. 
1997). 


4.  DISCUSSION 

(a)  Relevance  of  proposed  mechanisms  for  the 
entropic  (ideal)  elasticity  of  elastin 

The  kernel  of  the  classical  theory  of  rubber  elasticity 
originated  with  K.  H.  Meyer  in  1932  which,  according  to 
Flory  (1968,  p.  G42),  can  be  stated  as  ‘...  the  orientation 
of  the  molecules  of  a  piece  of  rubber  when  it  is  stretched 
would  entail  a  decrease  in  entropy,  and  thus  could  account 
for  the  retractive  force.’  This  was  restated  by  Flory  (1968) 
in  terms  of  two  subsidiary  hypotheses:  (i)  the  restoring 
force  originates  within  the  chain  molecules  by  distortion 
of  the  distribution  of  their  configurations  and  not  from 
interactions  between  chains;  and  (ii)  the  restoring  force  is 
due  ‘exclusively  to  the  entropy  of  the  chain  configuration’. 
However,  is  the  alignment  of  chains  really  the  most  appro¬ 
priate  description  of  entropic  elasticity?  Even  in  the  case 
of  a  bulk  elastomer,  might  the  most  fundamental  state¬ 
ment  be  the  decrease  in  internal  chain  dynamics  within 
each  chain,  which  would  be  fundamental  to  the  alignment 
of  elongation?  Clearly,  the  above-demonstrated  entropic 
elasticity  exhibited  by  a  single  chain  does  not  arise  from 
the  alignment  of  many  chains  on  extension.  Rather,  the 
observation  of  entropic  elasticity  on  extension  of  a  single 
chain  fixed  at  its  ends  limits  the  source  of  entropy  change 
to  the  damping  of  internal  chain  dynamics  within  the  sin¬ 
gle  chain.  In  this  case  entropic  elasticity  cannot  arise  from 
the  change  in  relative  orientation  of  many  chain  mol¬ 
ecules. 
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(«) 


temperature  (°C) 

Figure  8.  Thermoelasticity  curves  for  (a)  fibrous  elastin  and 
(: b )  20  Mrad  cross-linked  poly(GVGVP).  See  text  for 
discussion.  Whether  or  not  a  near-zero  slope,  i.e.  an  /Jf 
«  1,  is  obtained  depends  on  the  rate  at  which  the 
temperature  is  raised.  This  is  due  to  a  slow  denaturation  at 
temperatures  above  60  °C.  (Times  per  data  point  for  line 
segments  (a),  (b)  and  (c)  are  4  h,  30  min  and  4  h, 
respectively.)  Reproduced,  with  permission,  from  Urry 
(1988a,  b). 


(i)  Relevance  of  the  Floty  random-chain  network  theory  of 
entropic  elasticity 

Flory’s  publications  describe  a  Gaussian  distribution  of 
end-to-end  chain  lengths  that  results  from  a  random-chain 
network,  and  they  assert  that  the  decrease  in  entropy  due 
to  extension,  or  compression,  arises  from  distortion  of  the 
network  from  this  most  probable  random  distribution 
(Flory  1953,  1968).  Hoeve  &  Flory  (1958,  1974)  con¬ 
cluded  that  the  entropic  elasticity  of  elastin  arises  from 
such  a  random-chain  network. 

There  are  four  principal  reasons  why  this  perspective  is 
questionable.  These  are  (i)  a  limited  potential  for  cross¬ 
links  in  the  fixed  sequence  of  elastin,  (ii)  the  above  spec¬ 
troscopic  evidence  for  regular,  non-random,  dynamic 
structural  features  in  components  of  elastin,  (iii)  the  above 
report  of  the  occurrence  of  ideal  reversible  elasticity  exhib¬ 
ited  by  a  single  chain  of  repeat  sequences  of  elastin  with 
calculated  single-chain  elastic  moduli  similar  to  the 
macroscopic  elastic  modulus,  and  (iv)  an  irreversible  loss 
of  force  due  to  thermal  denaturation  when  the  elastomer 
is  extended  and  heated  at  temperatures  above  60  °C. 


Relevance  of  the  probability  of  cross-link  formation  to 
structure 

In  rubber,  each  repeating  unit  is  a  potential  cross-link, 
and,  on  vulcanization,  any  pair  of  repeating  units  can  be 
converted  to  a  cross-link  wherever  two  chain  segments  are 
sufficiently  proximal.  In  elastin,  only  lysine  residues  form 
the  cross-links,  and  yet  there  are  only  some  40  lysine  resi¬ 
dues  for  an  elastin  protein  of  ca.  800  residues,  i.e.  one 
lysine  in  every  20  residues.  Furthermore,  the  lysine  resi¬ 
dues  are  arranged  in  an  orderly  manner  along  the 
sequence.  Most  commonly,  four  lysine  residues  combine 
to  form  a  single  cross-link,  called  a  desmosine,  and  essen¬ 
tially  all  of  the  lysine  residues  in  fibrous  elastin  have  been 
formed  into  cross-links.  Such  cross-linking  cannot  be 
achieved  within  a  random-chain  network.  The  achieve¬ 
ment  of  such  rare  cross-linkages  requires  elements  of  order 
(Urry  1976). 

Low-frequency  motions  provide  much  entropy  to  chain 
structures 

Motions  that  are  due  to  rotations  about  backbone  bonds 
in  a  chain  molecule  that  exhibit  intense  relaxations  limited 
to  localized  frequency  intervals,  that  is,  intense  mechanical 
resonances,  can  only  result  from  polymers  containing  reg¬ 
ular,  non-random,  dynamic  structures.  For  model  elastin 
systems  including  purified  fibrous  elastin,  such  mechanical 
relaxations  are  observed  at  ca.  1  kHz  and  5  MHz  (see 
figures  6  and  7).  As  will  be  demonstrated  below,  these 
relaxations  provide  a  very  ample  source  of  entropy  to 
structured  elastin  systems.  These  contributions  to  entropy 
increase  greatly  as  the  frequency  of  the  relaxation  is  lower, 
and  they  are  contributions  to  the  free  energy  of  protein 
structure  that  have  yet  to  be  adequately  treated  in  the  cur¬ 
rent  computations  of  low-energy  protein  structures.  Most 
importantly  for  entropic  elasticity,  low-frequency  motions 
provide  an  abundant  and  ready  source  of  entropy  decrease 
on  extension. 

Entropic  elastomeric  force  arising  from  a  single  chain 

At  a  fixed  extension  in  the  AFM  single-chain  force- 
extension  curves  of  figure  3,  there  is  a  single  end-to-end 
chain  length.  The  entropic  elastic  force  measured  at  a 
given  extension  cannot  have  arisen  from  displacement  of 
a  network  of  chains  from  a  Gaussian  distribution  of  end- 
to-end  chain  lengths  nor  can  it  arise  from  aligning  chains 
in  parallel.  The  decrease  in  entropy  that  provides  the 
entropic  restoring  force  in  the  single-chain  experiment 
must  arise  from  a  damping  of  internal  chain  dynamics  of 
the  single  chain.  No  other  degrees  of  freedom  can  cause 
a  decrease  in  chain  entropy  to  occur.  Furthermore,  in 
order  that  the  ends  of  the  chain  remain  fixed  in  space,  the 
rotation  about  a  given  bond  must  be  compensated  for  by 
a  corresponding  correlated  rotation  about  one  or  more 
additional  bonds.  The  peptide  librational  mechanism  of 
entropic  elasticity,  which  arises  so  readily  out  of  the 
p-spiral  structure  of  poly(GVGVP)  represents  one  inter¬ 
esting  means  of  achieving  this. 

Studies  of  the  temperature  dependence  of  force  at  a 
fixed  extension  have  yet  to  be  carried  out  on  single  chains 
in  order  to  determine  the  contribution  of  internal  energy 
or  the  enthalpy  component  of  the  force,  i.e.  to  determine 
the  fs/f  ratio  as  given  in  equations  (1.7)  and  (1.8).  In 
respect  of  the  internal  energy  or  enthalpic  contribution  to 
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Figure  9.  Description  of  the  proposed  molecular  structure  of  poly(GVGVP).  (a)  Schematic  representation  of  the  repeating 
Pro-Gly-containing  sequence  inserting  the  p-tum  structure,  which  is  given  in  crystallographic  detail  in  part  (fc).  (c)  Schematic 
helical  representation  of  the  structure  that  forms  on  raising  the  temperature  above  the  inverse  temperature  transition.  ( d ) 
Schematic  representation  of  the  helical  structure,  called  a  p-spiral,  but  with  the  P-turns  included  and  functioning  as  spacers 
between  turns  of  the  helix,  (e)  Stereo  pairs  of  the  p-spiral  structure  in  atomic  detail  in  (i)  end  view  and  (ii)  side  view.  (/) 
Associated  p-spirals  having  formed  multi-stranded  twisted  filaments,  which  is  the  more  accurate  description  of  the 
hydrophobically  folded  and  assembled  state. 


the  elastic  force,  these  would  arise  reversibly  without  bond 
breakage,  by  strain-imposed  distortions  in  bond  lengths 
and  bond  angles.  Interestingly,  Young  (private 
communication)  probed  the  strain  deformation  of  silk 
fibres  (Bombyx  mori  and  spider  dragline)  using  Raman 
spectroscopy  and  observed,  among  other  spectral  changes, 
a  shift  in  the  1095  cm-1  (wavenumber)  band  on  extension 
to  3.5%.  This  band  shift,  assigned  to  the  C-C  stretch 
mode,  represents  an  internal  energy  component  to  the 
force.  Provided  that  the  stretch  is  reversible  with  no  bond 
breakage,  it  can  contribute  to  the  force  of  an  ideal  elasto¬ 
mer. 

At  high  extensions,  where  the  force  resulting  in  bond 
deformations  can  become  significant,  however,  distortions 
in  bond  lengths  and  bond  angles  can  decrease  the  dura¬ 
bility  of  an  elastic  fibre  by  chain  rupture.  In  this  regard, 
the  remarkable  feature  of  the  mammalian  elastic  fibre  is 
its  striking  durability  where,  in  the  aortic  arch  and 
descending  thoracic  aorta,  the  half-life  is  some  70  years 


during  which  period  the  elastic  fibres  will  have  undergone 
more  than  one  billion  stretch-relaxation  cycles.  This  is 
possible  due  to  the  dominantly  entropic  nature  of  its  elas¬ 
tic  force  as  seen  in  figure  8 b  and  due  to  uniform  chain 
extension  arising  out  of  a  regular,  cross-linked  structure 
where  no  single  chain  bears  an  excessive  amount  of  the 
load. 

Loss  of  elastomeric  force  of  elastin-based  systems  due  to 
thermal  denaturation 

In  figure  1,  it  was  noted  that  poly(GVGVP)  exhibits  a 
slow  irreversible  denaturation  when  held  at  temperatures 
above  60  °C.  This  same  phenomenon  is  apparent  in  the 
thermoelasticity  studies  on  fibrous  elastin,  as  originally 
noted  by  Dorrington  &  McCrum  (1977).  In  the  thermo¬ 
elasticity  experiments  illustrated  in  figure  8  (Urry 
1 988a,£>),  the  elastomer  is  equilibrated  at  40  °C  and  at 
that  temperature  is  extended  to  60%.  At  this  fixed  length, 
the  temperature  is  lowered  to  ca.  20  °C  and  slowly  raised 
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to  80  °C.  The  force  dramatically  increases  as  the  tempera¬ 
ture  is  raised  through  the  range  of  the  inverse  temperature 
transition  from  20  to  40  °C.  Then,  the  force/7  (°K)  tends 
to  remain  nearly  constant  from  40  to  60  °C,  especially  for 
X20-poly(GVGVP).  The  value  of  InfIT  begins  to  decrease 
above  60  °C  in  a  time-dependent  manner.  When  using  a 
rate  of  increase  of  4°C  30  min-1,  S^h"1,  there  is  a 
decrease  in  \nfIT ,  and,  when  the  temperature  is  returned 
to  40  °C,  an  irreversible  loss  of  force  occurs.  The  decrease 
in  force  above  60  °C  is  more  dramatic  on  raising  the  tem¬ 
perature  more  slowly,  e.g.  when  raising  the  temperature 
at  a  rate  of  1  °C  h-1.  When  the  temperature  maintained 
above  60  °C  for  a  longer  interval  irreversible  thermal 
denaturation  is  seen  most  dramatically.  Furthermore,  a 
plot  of  ln(elastic  modulus)  versus  time  at  80  °C  is  linear, 
providing  evidence  of  a  first  order  process  (Urry  1988a, 6). 
Interestingly,  while  heating  above  60  °C  decreases  force, 
it  also  shortens  the  length  at  which  zero  force  is  obtained, 
as  is  to  be  expected  when  a  filamentous,  p-spiral-like 
structure  becomes  randomized  by  thermal  denaturation. 

As  discussed  in  relation  to  equations  (1.7)  and  (1.8),  a 
slope  of  zero  for  plots  like  those  in  figure  8  would  mean 
an  ideal  (entirely  entropic)  elastomer.  As  clearly  seen  for 
X20-poly(GVGVP)  in  the  40-60  °C  temperature  range, 
when  data  are  collected  at  the  rate  of  30  min  per  data 
point,  a  near-zero  slope  is  obtained.  Commonly,  an 
ratio  of  0.1  is  obtained  for  X20-poly(GVGVP)  (Luan  etal 
1989).  Because  of  the  steep  positive  slope  between  20  and 
40  °C,  due  to  the  hydrophobic  folding  and  assembly  of 
the  inverse  temperature  transition,  and  because  of  the  slow 
irreversible  thermal  denaturation  above  60  °C,  the  slope, 
obtained  in  the  thermoelasticity  experiment  from  which 
the  entropic  contribution  is  estimated,  depends  on  the  rate 
at  which  the  data  are  collected  and  the  temperature  inter¬ 
val  over  which  they  are  collected. 

A  steep  rise  in  force  in  the  20-40  °C  temperature  inter¬ 
val  for  fibrous  elastin  in  figure  8a,  although  less  steep  than 
that  seen  in  figure  8 b  for  X20-poly (GVGVP),  demon¬ 
strates  the  effect  of  the  inverse  temperature  transition 
within  elastin  itself.  More  dramatic  for  fibrous  elastin  is 
the  loss  of  elastomeric  force  on  raising  the  temperature 
above  60  °C.  This  finding,  of  a  slow  loss  of  force  at  elev¬ 
ated  temperature,  would  not  have  been  so  apparent  if  the 
time  had  been  limited  to  a  single  day  to  perform  the 
experiment  (without  the  advantage  of  computerized 
equipment),  as  was  the  case  when  Hoeve  &  Flory  (1958) 
carried  out  their  original  study  on  fibrous  elastin.  For  the 
thermoelasticity  study  of  fibrous  elastin  shown  in  figure  8, 
however,  the  rate  of  loss  of  force  at  elevated  temperature 
is  even  greater  than  for  X20-poly (GVGVP). 

Additionally,  the  plot  of  ln(elastic  modulus)  versus  time 
for  fibrous  elastin  maintained  at  80  °C,  as  determined  by 
a  stress-strain  curve  every  24  h,  is  a  straight  line  with  a 
half-life  of  ca.  10  days  (Urry  19886).  Another  important 
feature  of  the  data  for  determining  half-life  is  that  the 
length  at  which  zero  force  is  obtained  decreases  with  time 
at  80  °C.  Remarkably,  this  means  that  the  loss  of  force  in 
figure  8,  and  the  progressive  loss  of  elastic  modulus  exhib¬ 
ited  by  a  series  of  stress-strain  curves  after  heating  for  24  h 
at  80  °C,  are  not  due  to  structural  rearrangement  to  an 
extended  state  under  load.  Additionally,  it  cannot  be  due 
to  slipping  of  the  sample  in  the  grip.  Instead,  extended 
time-periods  at  a  high  temperature  result  in  a  denatur¬ 


ation,  which  entails  a  randomization  from  linear  filamen¬ 
tous  structures  towards  denatured  random-coil  structures. 
This  has  an  analogy  to  the  loss  of  percentage  water  by 
weight  for  time  above  60  °C  in  figure  1 . 

Figure  8  clearly  displays,  in  the  thermoelasticity  data  of 
elastin-based  systems,  in  the  half-life  studies  for  the  elastic 
modulus  at  80  °C,  and  in  the  structural  change  above 
60  °C  in  figure  1,  a  slow  thermal  denaturation  of  protein- 
based  elastomers.  As  random-chain  networks  do  not  exhi¬ 
bit  thermal  denaturation,  these  elastin-based  systems  can¬ 
not  properly  be  described  as  random-chain  networks. 
Thus,  the  evidence  is  clear  and  abundant;  elastin-based 
systems  contain  regular  repeating  conformations  and 
those  ordered  structures  denature. 

(ii)  Relevance  of  solvent  entropy  changes  as  a  source  of 
entropic  elastic  force 

Effect  of  removing  solvent  entropy  change  during  the  inverse 
temperature  transition 

The  change  in  solvent  entropy  is  fundamental  to  the 
occurrence  of  an  inverse  temperature  transition.  This 
transition  has  been  extensively  utilized  to  perform  the  fam¬ 
ily  of  energy  conversions  that  sustain  living  organisms 
(Urry  1993,  1997),  and  therefore  the  inverse  temperature 
transition  is  extensively  characterized.  One  approach  to 
testing  whether  the  change  in  solvent  entropy  contributes 
to  the  elastomeric  force  is  to  utilize  an  approach  due  to 
Flory,  that  is,  to  reduce  the  heat  of  the  inverse  temperature 
transition  to  near  zero  by  an  appropriate  solvent  mixture. 
Hoeve  &  Flory  (1958)  used  a  mixture  of  30%  ethylene 
glycol  by  weight  in  water  in  their  thermo  elasticity  study 
that  deduced  a  dominant  entropic  elastic  force  for  fibrous 
elastin.  When  Luan  et  al.  (1989)  compared  the  thermo- 
elasticity-derived  f^f  ratio  of  X20-poly(GVGVP)  in  water 
and  in  30%  ethylene  glycol,  the  values  were  the  same,  0.1. 
Based  on  the  approximations  of  equations  (1.7)  and  (1.8), 
X20-poly(GVGVP)  is  a  dominantly  entropic  elastomer  in 
both  solvent  systems.  The  effect  of  adding  the  ethylene 
glycol  is  to  reduce  the  heat  of  the  transition  to  near  zero. 
Now,  since  the  entropy  of  the  desolvation  process  is  the 
heat  of  the  transition  divided  by  the  temperature  of  the 
transition,  A HJTU  the  solvent  entropy  change  would  be 
approaching  zero,  and  yet  the  solvent  mixture,  instead  of 
reducing  the  elastomeric  force,  remarkably,  did  the 
opposite  and  resulted  in  a  substantial  increase  in  force. 
The  conclusion  is  that  there  is  no  experimental  basis  for 
believing  that  solvent  entropy  change  contributed  to  elas¬ 
tomeric  force. 

The  polymer  chains  bear  the  force ;  solvation  changes  alter 
the  fraction  of  chain  segments  that  sustain  the  force ? 

The  force  of  deformation  must  be  borne  by  the  chains; 
it  cannot  be  borne  by  the  solvent.  At  a  fixed  length,  raising 
the  temperature  of  a  band  of  elastomer  from  below  to 
above  the  temperature  of  the  inverse  temperature  tran¬ 
sition  increases  the  force.  Alternatively,  increasing  the 
chemical  potential  to  lower  the  transition  temperature 
from  above  to  below  the  operating  temperature  at  a  fixed 
length  causes  an  increase  in  the  elastomeric  force  (Urry 
1997;  Urry  et  al.  1988a, 6).  These  results  are  due  to  an 
increase  in  hydrophobic  folding  and  assembly  of  the 
chains.  As  this  is  the  case  at  a  fixed  length,  the  smaller 
fraction  of  a  chain’s  length,  which  is  not  folded,  must 
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(a) 


Figure  10.  Stereo  view  of  the  torsional  oscillations  possible 
within  a  2  kcal  mol-1  residue  cut-off  energy  for  the  central 
pentamer  of  one  turn  (three  pentamers)  of  the  p-spiral 
structure  of  poly(GVGVP)  (a)  in  the  relaxed  state  and  ( b )  in 
the  130%  extended  state.  See  text  for  discussion. 
Reproduced,  with  permission,  partly  from  Urry  et  al.  (1982) 
and  the  remainder  from  Urry  et  al  (1985c). 


necessarily  be  further  extended.  Accordingly,  in  order  to 
maintain  a  fixed  length,  the  unfolded  segments  must  be 
further  extended.  This  provides  for  the  increase  in 
entropic  elastic  force  by  means  of  a  decrease  in  configur¬ 
ational  entropy. 

(iii)  Entropic  elasticity  due  to  damping  of  internal  chain 
dynamic  on  extension 
Molecular  structure  of  poly  (GVGVP) 

In  order  to  understand  the  concept  of  the  damping  of 
internal  chain  dynamics  on  extension,  the  structure  given 
in  figure  9  is  helpful.  The  Pro-Gly-containing  elastomeric 
sequential  polypeptides  exhibit  a  repeating  p-turn 
structure  as  schematically  shown  in  figure  9a,  and  the 
p-tum  is  given  in  crystallographic  detail  in  figure  9b.  A  rise 
in  the  temperature  of  the  inverse  temperature  transition 
results  in  the  optimization  of  intra-  and  intermolecular 
hydrophobic  association.  Intramolecularly,  this  results  in 
a  helical  structure  called  a  p-spiral,  which  is  shown 
schematically  as  a  helix  in  figure  9c,  in  more  detail,  with 
P-tums  schematically  shown  as  spacers  between  turns  of 
the  helix,  in  figure  9 d,  and  in  atomic  stereo  detail  in  figure 
9e.  Intermolecularly,  the  p-spirals  associate  to  form  multi- 
stranded  twisted  filaments,  as  depicted  in  figure  9f 


Because  there  is  water  within  the  P-spiral,  and  as  a  result 
of  the  p-tums  functioning  as  spacers,  the  peptide  segments 
connecting  P-tums  are  free  to  undergo  large  amplitude 
torsional  oscillations. 

Effect  of  extension  on  the  amplitude  of  the  torsional 
oscillations  and  on  chain  entropy 

The  magnitude  of  the  torsional  (dihedral  angle)  excur¬ 
sions  for  a  single  pentamer  sandwiched  within  a  turn  of 
p-spiral  (three  pentamers)  for  a  1.5  kcal  mol-1  residue 
cut-off  energy  is  shown  in  figure  10a  when  using  the  mol¬ 
ecular  mechanics  calculations  of  the  Scheraga  approach 
(Momany  et  al.  1975).  Dihedral  angle  oscillations  of  160° 
are  possible.  The  effect  of  a  130%  extension  of  the  p-spiral 
pitch,  shown  in  figure  106,  is  a  dramatic  decrease  in  the 
amplitude  of  the  torsional  oscillations  allowed  within  a 

1.5  kcal  mol-1  residue  cut-off  energy.  Using  an  enumer¬ 
ation  of  states  approach  with  each  change  of  5°  in  any 
single  torsion  angle  counted  as  a  new  state,  the  decrease 
in  numbers  of  states  for  a  0.6  kcal  mol"1  residue  cut-off 
energy  (i.e.  within  an  energy  cut-off  of  RT)  is  from  a  WT 
of  342  to  a  IVs  of  24.  Using  equation  (1.10),  this  gives  a 
value  for  A S  of  -5.28  cal  mol"1  deg"1  per  pentamer.  If 
the  Boltzmann  summation  over  all  energies  is  used,  the 
value  is  -5.06  cal  mol"1  deg"1  per  pentamer  (Urry  et  al. 
1985 c). 

Alternatively,  when  using  the  molecular  dynamics 
approach  of  Karplus,  the  expression  becomes 

HS  =  R  ]n{nA<t>fA^fiUA<t>\A^y  (3.1) 

where  the  product  IIj  is  over  the  rms  torsion  angle  fluctu¬ 
ations  A</>f  and  A^f  for  the  extended  state  divided  by  the 
product  of  the  rms  torsion  angle  fluctuations  A</>[  and 
Ai for  the  relaxed  state  as  obtained  from  the  molecular 
dynamics  trajectories  for  the  relaxed  and  130%  extended 
states.  The  calculated  decrease  in  entropy  was  - 

5.5  cal  mol"1  deg"1  per  pentamer  for  the  130%  extension 
(Chang  &  Urry  1989).  Whether  the  calculation  of  entropy 
is  by  means  of  an  enumeration  of  states  in  i /*-<!>  configur¬ 
ation  space  or  the  amplitude  of  the  torsional  oscillations 
in  a  molecular  dynamic  simulation,  the  calculated  entropy 
change  for  a  130%  extension  is  the  same. 

Calculation  of  contribution  of  force  due  to  decrease  in 
freedom  of  motion  on  extension 

From  equations  (1.4)  and  (1.5)  we  have 

fs=ndS/dL)VtT.  (3.2) 

At  physiological  temperatures,  310  K,  and  for  an  exten¬ 
sion  from  0.35  to  0.8  nm  for  a  pentamer  within  the  P- 
spiral  structure  of  figure  9e  with  a  calculated  decrease  in 
entropy  of  5  cal  mol"1  deg”1  per  pentamer,  the  calculated 
value  for  the  entropic  contribution  to  the  elastomeric  force 
is  24  pN  for  a  single  chain  of  p-spiral. 

To  put  this  in  terms  of  an  elastic  modulus  in  units  of 
N  m"2,  it  is  necessary  to  estimate  the  number  of  p-spirals 
per  m2.  From  figure  1,  the  functional  state  is  37%  peptide 
by  weight.  Taking  the  density  of  the  peptide  to  be 
1.33  gem"3,  and  a  cross-sectional  area  of  3.2 nm2  per 
p-spiral,  there  would  be  ca.  8.7  x  1016  parallel-aligned 
P-spirals  per  m2.  The  calculated  result  is  2x  106Nm"2, 
which  is  an  order  of  magnitude  larger  than  the 
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experimental  value  of  1 .6  x  105  N  m  2  for  20  Mrad  cross- 
linked  (GVGVP)251. 

Several  obvious  factors  can  be  considered,  that  would 
contribute  to  a  higher  calculated  value  for  the  elastic 
modulus  of  the  isolated  (5-spiral  than  for  the  experimental 
elastic  modulus  of  the  cross-linked  matrix.  First,  in  the 
cross-linked  matrix  the  [5-spirals  would  be  randomly  orien¬ 
tated  with  respect  to  the  direction  of  extension.  Second, 
the  (5-spirals  occur  as  part  of  a  multistranded  twisted 
filament  as  shown  in  figure  9/,  which  would  tend  to 
dampen  the  torsional  oscillations  within  the  triple- 
stranded  filament  and  thereby  result  in  a  lower  entropy 
change  on  extension.  Thus,  it  would  seem  that  the  damp¬ 
ing  of  internal  chain  dynamics  on  extension  provides 
abundant  entropic  restoring  force  in  elastin-based  sys¬ 
tems.  Furthermore,  given  the  entropic  elasticity  of  single 
chains  demonstrated  by  the  AFM  single-chain  force- 
extension  curves,  damping  of  internal  chain  dynamics  on 
extension  should  be  considered,  in  general,  as  a  potential 
fundamental  source  for  the  development  of  entropic  elas¬ 
tic  force  in  chain  molecules. 

(iv)  Dependence  of  entropy  and  structural  free  energy  on 
oscillator  frequency 

One  final  point  for  consideration  in  this  report  is  the 
significance  of  the  occurrence  of  intense  low-frequency 
mechanical  motions  in  elastin-based  systems.  The  insight 
comes  from  the  expression  for  the  dependence  of  entropy 
on  oscillator  frequency.  From  the  expression  for  the  par¬ 
tition  function  of  the  harmonic  oscillator  (Dauber  et  al . 
1981),  it  is  seen  that  the  slope  for  the  frequency  depen¬ 
dence  of  entropy  is  -4.6  EU  log^r  !.  The  extrapolation  of 
the  harmonic  oscillator  treatment  to  low-frequency  tor¬ 
sional  oscillations  is  surely  not  quantitative,  but  it  is  none¬ 
theless  informative.  While  the  contribution  of  a  lower 
classical  vibration  frequency  is  of  the  order  of 
5  cal  mol-1  deg~!  (EU),  giving  a  Gibbs  free  energy  of 
LSkcalmol"1  at  298  K,  the  contribution  of  the  5  MHz 
relaxation,  seen  in  figures  6  and  7,  would  approach 
30  cal  mol-1  deg-1  or  9  kcal  mol-1,  and  that  for  the  mech¬ 
anical  resonance  at  ca.  1  kHz,  reported  in  figure  6a,  would 
approach  47  cal  mol-1  deg-1  or  Hkcalmol-1.  From  this 
it  should  be  apparent  that  these  relaxations  are  fundamen¬ 
tal  in  providing  for  the  entropic  elasticity  of  elastin-based 
systems.  On  a  more  general  note,  computations  of  protein 
structure  that  do  not  include  these  important  low- 
frequency  contributions  to  the  free  energy  would  seem  to 
be  limited  in  their  treatment  of  the  structure-function 
issue. 

The  work  was  supported,  in  part,  by  the  Office  of  Naval 
Research  under  contract  numbers  N00014-00-C-0178  and 
N000 14-00-C-0404. 
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ABSTRACT 

During  the  last  half  century,  identification  of  an  ideal  (predominantly  entropic)  protein 
elastomer  was  generally  thought  to  require  that  the  ideal  protein  elastomer  be  a  random  chain 
network.  Here,  we  report  two  new  sets  of  data  and  review  previous  data.  The  first  set  of  new 
data  utilizes  atomic  force  microscopy  to  report  single-chain  force-extension  curves  for 
(GVGVP)25i  and  (GVGIP)26q>  and  provides  evidence  for  single-chain  ideal  elasticity.  The 
second  class  of  new  data  provides  direct  contrast  between  low  frequency  sound  absorption  (0.1 
to  10  kHz)  exhibited  by  random  chain  network  elastomers  and  by  elastin  protein-based  polymers. 

Earlier  composition,  dielectric  relaxation  (1  to  1000  MHz),  thermoelasticity,  molecular 
mechanics  and  dynamics  calculations  and  thermodynamic  and  statistical  mechanical  analyses  are 
presented  that  combine  with  the  new  data  to  contrast  with  random  chain  network  rubbers  and  to 
detail  the  presence  of  regular  non-random  structural  elements  of  the  elastin-based  systems  that 
lose  entropic  elastomeric  force  on  thermal  denaturation. 

The  data,  analyses  and  calculations  affirm  an  earlier  contrary  argument  that  components  of 
elastin,  the  elastic  protein  of  the  mammalian  elastic  fiber,  and  purified  elastin  fiber  itself  contain 
dynamic,  nonrandom,  regularly  repeating  structures  that  exhibit  dominantly  entropic  elasticity  by 
means  of  a  damping  of  internal  chain  dynamics  on  extension. 
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I 


1.  Introduction 

1.1  Definition  of  ideal  or  perfect  elasticity 

Ideal  elasticity  is  the  property  whereby  the  energy  expended  in  deformation  of  the  elastomer 
is  completely  recovered  on  removal  of  the  deforming  force.  Since  the  energy  expended  in 
deformation  is  given  by  the  area  under  the  force,  f,  versus  increase  in  length,  AL,  curve,  a 
perfectly  reversible  force-extension  curve  means  complete  recovery  on  relaxation  of  the  energy 
expended  on  deformation.  Therefore,  ideal  elastomers  exhibit  perfectly  reversible  force- 
extension  curves. 

Perhaps  our  earliest  perspective  of  the  mechanism  underlying  ideal  elasticity  comes  from  a 
fundamental  observation  concerning  rubber  elasticity.  In  the  mid  nineteenth  century,  Joule  and 
Thomson  noted  a  quantitative  correlation  between  the  increase  in  temperature  of  the  elastomer 
due  to  stretching  and  the  increase  in  force  due  to  increasing  the  temperature  (Flory  1968). 
Thermodynamics  provides  for  the  analysis  underlying  this  correlation,  and  the  Boltzmann 
relation  provides  the  bridge  between  experimental  thermodynamic  quantities  and  statistical 
mechanical  description  of  molecular  structures. 

Continuing  qualitatively  with  the  Joule  and  Thomson  correlation,  heat  produces  motion,  and 
the  energy  represented  by  heat  distributes  into  the  various  available  degrees  of  freedom  in  the 
chain  molecules  comprising  the  elastomer.  Accordingly,  the  release  of  heat  on  stretching 
correlates  with  a  loss  of  motion.  By  means  of  statistical  mechanics,  the  loss  of  motion  is  seen  as 
a  decrease  in  entropy  on  extension.  In  addition,  should  solvent  be  essential  for  elasticity,  this 
requires  explicit  consideration. 

But  there  is  more  to  elastomeric  force  than  entropy  changes  arising  from  changes  in  mobility. 
There  is  also  internal  energy,  and  a  proper  understanding  of  elasticity  requires  delineation  of 
internal  energy  and  entropy  components  of  the  elastomeric  force.  As  the  internal  energy 
component  of  the  elastic  force  increases,  irreversible  processes  such  as  chain  breakage  become 
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more  probable.  So  greater  durability  results  as  elastomers  become  more  dominantly  entropic  and 
durability  becomes  a  feature  of  an  ideal  or  more  perfect  elastomer. 


1.2  Delineation  of  internal  energy  and  entropy  components  of  elastomeric  force 

Delineation  of  internal  energy  and  entropy  components  of  elastomeric  force  can  begin  with 
the  definition  of  the  Helmholtz  free  energy,  A,  also  referred  to  as  the  maximal  work  function, 

A  s  E  -  TS,  (1) 

where  E  is  the  internal  energy,  T  is  the  absolute  temperature  (degrees  K),  and  S  is  the  entropy  for 
the  system.  The  differential  of  the  maximum  work  function,  dA,  can  be  written  to  include  the 
work  done  on  an  elastomer  by  application  of  a  force,  f,  over  the  change  in  length,  dL,  i.e., 

dA  =  -PdV-SdT  +  fdL  (2) 

where  P  is  the  pressure,  and  V  is  the  volume.  Writing  the  partial  differential  of  A  in  Eq.  (1)  with 
respect  to  length  at  constant  V,  T  and  composition,  n,  gives 

0A/0L)VXn  =  OE/aL)VXn  -  T0S/3L)VXn,  (3) 

but  by  Eq.  (2),  (5A/3L)VTn  =  f,  such  that 

f  =  0E/8L)VXn  -  T0S/3L)VXn.  (4) 

Accordingly,  the  force  is  seen  to  be  comprised  of  two  components,  an  internal  energy 
component,  fE,  and  an  entropy  component,  fs,  i.e., 

f  =  fE  +  fs.  (5) 

Now  it  becomes  useful  to  find  an  expression  that  will  allow  experimental  estimation  of  the 
relative  magnitude  of  the  internal  energy  and  entropy  components  of  the  force.  Following  Flory 
et  al  (1960),  when  the  functions  exist  and  are  continuous,  the  order  of  a  partial  differential  does 
not  matter,  and  it  can  be  shown  that  0S/3L)vt  =  (3f/dT)V  L,  such  that, 

f  =  0E/3L)VX„  +  T(9f/3T)VLn.  (6) 

Interestingly,  this  can  be  rewritten  as 

fE/f  =  -  T(3ln[f/T]/9T)VLn,  (7) 
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which  allows  for  an  experimental  estimate  of  the  fE/f  ratio  from  the  slope  of  a  plot  of  ln(f/T) 
versus  temperature  under  conditions  of  constant  V,  L  and  composition,  n. 

By  means  of  an  approximate  correction  term,  the  fE/f  ratio  can  be  estimated  under  more  usual 
experimental  conditions  of  constant  P,  L  and  at  equilibrium,  e.g.,  with  surrounding  solvent,  i.e., 

fE/f  =  -T(9ln[f/T]/3T)p  Leq,-  PeqT/[a3(V/V)  -  1]  (8) 

where  Peq  =  (3lnV/3T)PLeq  is  the  thermal  expansion  coefficient;  a  is  the  fractional  increase  in 
length,  L/Lj,  with  the  subscript  i  indicating  initial  length;  and  Vj  and  V  are  the  volumes  of  the 
elastomer  before  and  after  elongation,  respectively  (Dorrington  &  McCrum  1977;  Hoeve  and 
Flory  1962).  The  correction  term  was  derived  under  the  assumption  of  random  chain  networks 
with  a  Gaussian  distribution  of  end-to-end  chain  lengths. 

As  we  will  see,  a  random  chain  network  comprised  of  a  Gaussian  distribution  of  end-to-end 
chain  lengths  between  cross-links  is  not  an  accurate  structural  description  of  the  elastin  protein 
nor  of  its  simpler  models.  Also,  there  are  important  solvation  changes  that  can  occur  under 
experimental  conditions.  Nonetheless,  the  insight  into  the  magnitude  of  the  entropic  component 
of  the  elastomeric  force  provided  by  Eq.  (8),  with  experimental  estimates  of  the  quantities  in  the 
correction  term  or  when  neglecting  the  correction  term  altogether,  is  informative,  but  only  when 
the  temperature  range  for  determination  of  slope  is  judiciously  chosen. 

1.3  Basic  statistical  mechanical  expression  for  entropy 

1.3.1  The  Boltzmann  relation 

The  starting  point  for  the  statistical  mechanical  expressions  of  entropy  is  the  Boltzmann 
relation, 

S  =  RlnW  (9) 

R  (1.987  cal/degree-mole)  is  the  gas  constant,  R  =  Nk,  with  N  being  Avogadro’s  number  (6.02  x 
1023  /mole)  and  k  the  Boltzmann  constant  (1.38  x  10  16  erg/degree  K),  and  W  is  the  number  of  a 
priori  equally  probable  states  accessible  to  the  system  (Eyring  et  al  1964).  W  is  the  volume  in 
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phase  space  occupied  by  a  particular  state  of  a  molecular  system.  In  practice,  W  becomes  the 
product  of  partition  functions  for  each  of  the  degrees  of  freedom  of  the  molecular  system.  There 
are  3n  -  it  degrees  of  freedom  where  n  is  the  number  of  atoms  in  the  molecule  and  k  are  the 
holonomic  constraints  on  the  system.  In  general,  there  are  three  translational  and  three  rotational 
degrees  of  freedom,  and  the  remaining  are  vibrational  degrees  of  freedom.  For  large  enough 
molecules,  the  vibrational  degrees  of  freedom  also  include  motions  of  rotation  (or  torsional 
oscillations)  about  bonds.  Due  to  the  constraints  of  the  elasticity  measurement  k  is  6;  the  ends  of 
the  molecules  are  effectively  fixed  in  space  such  that  there  are  neither  translational  nor  whole 
molecule  rotational  degrees  of  freedom. 

1.3.2  Fundamental  expression  for  the  change  in  entropy  on  extension 

Fortunately,  when  calculating  entropic  contributions  to  the  elasticity,  interest  centers  on  the 
change  in  entropy  on  extension  from  a  relaxed,  Sr,  to  an  extended,  Se,  state  of  the  elastomer, 

AS  =  (Se  -  Sr)  =  Rln(We/Wr).  (10) 

Since  the  change  in  entropy  is  a  ratio,  it  becomes  relatively  straight-forward  to  calculate  the 
contribution  to  the  change  in  entropy  of  a  particular  expression  of  motion  accessible  to  the 
representation  of  the  molecular  structure  to  be  examined  in  its  relaxed  and  extended  states.  For 
our  purposes  the  P-spiral  structure  of  the  (Gly-Val-Gly-Val-Pro)n,  or  poly(GVGVP),  model  of 
elastin  is  used  and  two  different  means  of  representing  that  structure  are  considered.  The  first  is 
the  enumeration  of  states  in  configuration  (<j)-\|r  torsion  angle)  space  within  a  given  energy  cut-off 
for  relaxed  and  extended  states,  and  the  second  is  to  use  the  molecular  dynamics  approach  to 
determine  the  change  in  rms  torsional  oscillations  that  occur  on  extension  by  the  same  amount  as 
used  in  the  enumeration  of  states  approach.  The  magnitudes  of  the  calculated  entropy  changes 
on  extension  obtained  by  each  approach  can  then  be  compared. 

1.4  Historical  notes  of  proposed  mechanisms  for  protein  entropic  elasticity 
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1.4.1  Classical  (random  chain  network)  theory  of  rubber  elasticity 

During  the  last  half  of  the  last  century,  identification  of  ideal  (dominantly  entropic)  protein 
elasticity  required,  by  the  recognized  authorities  of  the  period,  the  conclusion  that  the  ideal 
elastomer  was  comprised  of  random  chain  networks  with  a  Gaussian  distribution  of  end-to-end 
chain  lengths  between  cross-links.  For  elastin,  the  protein  of  interest  here,  this  point  of  view 
began  with  a  paper  by  Hoeve  and  Flory,  (1958).  This  perspective  became  entrenched  in  the 
minds  of  the  interested  scientific  community  by  the  award  in  the  Fall  of  1974  of  the  Nobel  Prize 
to  Paul  Flory  following  reaffirmation  of  the  random  chain  network  conclusion  (Hoeve  and  Flory 
1974).  The  message  of  the  paper  is  unmistakable.  In  the  Synopsis  it  is  stated,  “4  network  of 
random  chains  within  the  elastic  fibers,  like  that  in  a  typical  rubber,  is  clearly  indicated. " 
Furthermore,  in  Figure  1  of  their  paper  is  a  structural  representation  of  the  chains  between  cross¬ 
links  with  a  statement  in  the  figure  legend  that,  “Configurations  of  chains  between  cross- 
linkages  are  much  more  tortuous  and  irregular  than  shown." 

1.4.2  Decrease  in  solvent  entropy  on  extension 

The  purpose  of  the  1974  Hoeve  and  Flory  paper,  which  presented  no  additional  data,  was  to 
refute  an  earlier  publication  by  Weis-Fogh  and  Anderson  (1970),  who  had  suggested  an  alternate 
mechanism  for  the  elasticity  of  elastin.  The  alternate  mechanism  proposed  that  hydration  of 
hydrophobic  side  chains  of  the  protein  that  become  exposed  to  solvent  on  extension  would  be 
responsible  for  the  stretch-induced  decrease  in  entropy,  that  is,  by  this  proposal  changes  in 
solvent  entropy  became  the  source  of  the  entropic  component  of  the  elastomeric  force. 

1.4.3  Damping  of  internal  chain  dynamics  on  extension 

Just  over  a  decade  later  yet  another  perspective  was  put  forward.  Considering  structural 
studies  and  molecular  mechanics  calculations  of  the  most  prominent  repeating  sequence  of 
bovine  elastin,  (Gly-Val-Gly-Val-Pro)n,  entropic  elasticity  was  described  as  arising  from  a 
decrease  in  available  configuration  space  on  extension  (Urry  et  al  1982).  Equivalently,  treating 


7 


the  experimentally  and  computationally  derived  regularly  repeating  structure  from  the 
perspective  of  molecular  dynamics,  the  damping  of  internal  chain  dynamics  on  extension 
described  the  same  decrease  in  entropy  for  the  same  degree  of  extension  (Chang  &  Urry  1989). 

In  the  Discussion  section  below,  these  three  mechanisms  will  be  considered  in  more  detail, 
following  consideration  of  new  experimental  data  on  models  of  elastin  along  with  previously 
published  relevant  data. 

1.5  Inverse  Temperature  Transition  Behavior  of  Elastin  and  its  Models 

1.5.1  Increase  in  order  with  increase  in  temperature 

1.5.1.1  Filament  Formation  on  Raising  the  Temperature:  When  the  temperature  of  an 
aqueous  solution  of  the  precursor  protein,  tropoelastin,  is  increased,  the  protein  aggregates  to 
form  a  more-dense  viscoelastic  phase,  called  a  coacervate.  When  a  droplet  of  the  aqueous 
suspension  of  incipient  aggregates  is  placed  on  a  carbon  coated  grid,  negatively  stained  with 
uranyl  acetate  and  oxalic  acid  at  the  appropriate  pH,  and  examined  in  the  transmission  electron 
microscope,  fibrils  comprised  of  parallel  aligned  filaments  are  observed  with  a  5  nm  periodicity 
(Cox  et  al  1974).  Similar  results  are  obtained  for  a-elastin,  a  70,000  Da  fragment  of  the  elastic 
fiber  (Cox  et  al  1973)  and  for  high  polymers,  poly(GVGVP)  and  poly(GVGVAP),  of  repeating 
sequences  of  elastin  (Volpin  et  al  1976).  Fibrils  of  similar  dimensions  have  been  reported  for 
fibrous  elastin  itself  using  similar  techniques  (Gotte  et  al  1974). 

1.5.1.2  Crystallization  of  Cyclic  Analogues  on  Raising  the  Temperature  of  the  inverse 
temperature  transition:  Most  strikingly,  cyclic  analogues  of  repeating  sequences  of  elastin 
crystallize  on  raising  the  temperature  (Urry  et  al  1978;  Cook  et  al  1980)  and  redissolve  on 
lowering  the  temperature.  Without  ambiguity,  these  polymers  increase  order  with  increase  in 
temperature  through  a  transition  temperature  range. 
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We  refer  to  this  phase  transition  as  an  inverse  temperature  transition.  While  the  peptide 
component  of  this  water-peptide  system  increases  in  order  on  increasing  the  temperature, 
ordered  water  molecules  surrounding  hydrophobic  residues  become  less  ordered  bulk  water  as 
hydrophobic  groups  separate  from  solution.  Accordingly,  the  overall  change  effected  by  the 
phase  transition  is  to  less  order  in  keeping  with  the  second  law  of  thermodynamics. 

Indeed  the  phase  diagram  for  high  molecular  weight  polymers  of  repeating  sequences  of 
elastin  is  inverted  with  the  soluble  phase  below,  at  lower  temperature,  and  the  insoluble  phase 
above  the  binodal  or  coexistence  line  and  with  an  inverted  coexistence  line,  i.e.,  with  the 
curvature  of  the  coexistence  line  convex  to  the  volume  fraction  axis  (Sciortino  et  al  1993;  Manno 
et  al  2001).  On  the  other  hand,  the  usual  circumstance  for  polymers  is  for  the  coexistence  line  to 
be  concave  to  the  volume  fraction  axis  and  for  the  polymer  to  be  insoluble  below  and  soluble 
above  the  coexistence  line  (Flory  1953). 

1.5.2  Composition  of  poly(GVGVP)  in  water  as  a  function  of  temperature 

The  temperature  dependence  of  composition  of  the  poly(GVGVP)-water  system  provides  a 
clear  visualization  of  a  water-containing,  structured  state  at  intermediate  temperatures  (Urry  et  al 
1985b).  A  plot  of  percent  water  by  weight  as  a  function  of  temperature  is  schematically  shown 
in  Figure  1.  As  the  polymer  is  miscible  with  water  in  all  proportions  below  20°C,  the  plot 
arbitrarily  starts  at  about  90%  water  by  weight  near  0°C.  On  raising  the  temperature  from  20  to 
30°C,  there  occurs  a  reversible  phase  separation  to  form  a  state  that  is  63%  water  by  weight.  The 
percent  water  by  weight  of  this  intermediate  water-retaining  state  remains  essentially  unchanged 
until  the  temperature  is  raised  above  60°C.  Then  a  veiy  slow  irreversible  transition  occurs  on 
raising  the  temperature  from  60  to  80°C  to  form  a  state  that  is  32%  water  by  weight.  Dialysis  of 
the  polymer  against  100,000  Da  cut-off  membranes,  after  the  prolonged  heating  at  80°C,  resulted 
in  less  than  0.1%  loss  of  polymer,  that  is,  prolonged  heating  did  not  result  in  chain  breakage. 
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The  results  on  poly(GVGVP)  are  compared  in  Figure  1  to  those  of  Grinberg  et  al  (1999)  on 
the  petroleum-based  polymer,  poly(N-isopropylacrylamide),  (PNIPAM).  This  polymer  shows  an 
inverted  phase  diagram  like  that  of  the  elastin-based  polymers.  Though  with  a  composition  more 
hydrophobic  than  that  of  poly(GVGVP),  the  transition  for  this  LCST  (lower  critical  solution 
temperature)  polymer  occurs  at  a  higher  temperature  (34°C),  is  much  sharper,  and  directly  results 
in  the  formation  of  the  state  with  approximately  30%  water  by  weight.  This  PNIPAM  state  of 
30%  water  by  weight  may  reasonably  be  described  as  a  random  chain  network.  It  is  clear  that 
poly(GVGVP)  forms  an  intermediate,  water  containing  state  that  holds  a  fixed  amount  of  water 
until  irreversible  denaturation  occurs  at  high  temperature.  This  state  of  63%  water  by  weight  is 
an  elementary  property  that  simply  is  not  consistent  with  the  recently  computed  “compact 
amorphous  globule”  for  (GVGVP)jg  due  to  Li  et  al  (2000). 

1.5.3  A  structured  state  at  intermediate  temperatures  for  elastin-based  systems 

Both  the  elastin  fiber  and  a-elastin  at  intermediate  temperatures  form  similar  water- 
containing  states  of  about  50%  water  by  weight  (Partridge  et  al  1955;  Partridge  &  Davis  1955; 
Partridge  1966,  1967).  It  will  also  be  shown  below  (see  Figure  9)  that  both  cross-linked 
poly(GVGVP)  and  fibrous  elastin  exhibit  a  slow,  irreversible,  loss  of  elastic  force  when  at  fixed 
extension  after  prolonged  exposure  to  temperatures  above  60°C.  These  are  further 
demonstrations  of  the  slow  irreversible  denaturation  observed  in  Figure  1  above  60°C  for 
poly(GVGVP).  Furthermore,  detailed  dielectric  relaxation  studies  on  the  elastin-based  systems  - 
fibrous  elastin,  a-elastin,  and  poly(GVGVP)  -  demonstrate  molecular  motions,  relaxations, 
limited  to  localized  frequency  ranges  (See  Figures  6,  7  and  8),  and  these  could  only  be  the  result 
of  non-random,  regularly  structured,  albeit  dynamic  conformational  states.  These  experimental 
results  are  consistent  with  neither  random  chain  networks  nor  compact  amorphous  globules. 
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2.  MATERIALS 


2.1  Preparation  of  elastin  models 

In  its  native  state  elastin,  a  protein  of  about  70,000  Da,  occurs  as  the  overwhelming  dominant 
component  (90+%  by  weight)  of  the  mammalian  elastic  fiber,  which  resides  within  the 
extracellular  matrix,  for  example,  as  a  structure  several  microns  in  diameter  in  the  vascular  wall. 
Non-elastin  components  are  trapped  enzymatic  protein  involved  in  the  cross-linking  and  a  fine 
fibrillar  coating  of  some  hundreds  of  nanometer  thickness.  The  elastic  fiber  is  a  substantially 
cross-linked,  insoluble  material  of  the  extracellular  matrix.  Efforts  to  isolate  it  for  purposes  of 
defining  its  properties  require  draconian  methods,  for  example,  refluxing  in  sodium  hydroxide,  to 
remove  all  of  the  surrounding  matrix  components.  This  product  is  called  purified  fibrous  elastin. 
Alternatively,  model  approaches  are  used  such  as  isolating  and  characterizing  (1)  the  precursor 
protein,  tropoelastin,  which  may  be  obtained  in  ill-defined  intermediate  states  of  cross-linking, 
(2)  the  pure  precursor  protein  prepared  by  microbial  biosynthesis,  (3)  chemical  degradation 
products  of  purified  fibrous  elastin,  and  (4)  chemically  and  microbially  synthesized  repeating 
sequences  found  within  the  protein  sequence  as  well  as  interesting  analogues. 

2.1.1  Synthesis  of  model  systems 

2.1.1.1  Chemical  Synthesis:  Initially  several  thousand  protein-based  polymers  were 
synthesized  in  order  to  determine  the  conformation  and  function  of  these  model  proteins.  This 
involved  the  synthesis  of  oligomers  and  high  polymers  of  repeating  tetra-,  penta-,  hexa-,  nona- 
and  even  decapeptides  reported  in  the  sequence  of  elastin.  By  means  of  chemical  synthesis, 
many  compositions  could  be  prepared  and  characterized,  and  axioms  were  developed  for  the 
function  of  these  elastic  model  proteins  in  energy  conversion  and  for  their  use  as  biomaterials  for 
medical  and  non-medical  applications. 
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2.1.1.2  Microbial  Biosynthesis  Once  a  specific  composition  was  identified  as  being  of 
sufficient  interest,  it  was  prepared  by  means  of  recombinant  DNA  technology  where  hundreds  of 
grams  could  be  prepared  in  a  single  fermentation  that  by  means  of  bench  chemical  synthesis 
would  require  much  of  a  year  to  prepare. 

2.1.2  Preparation  of  natural  materials 

2.1.2.1  Isolation  of  the  precursor  protein,  tropoelastin :  Initially,  the  precursor  protein 
was  prepared  from  animals  made  copper  deficient  and  treated  with  inhibitors  of  enzymatic  cross- 
linking.  The  product  was  meager  in  amount  and  a  mixture  of  partially  cross-linked  and  oxidized 
lysyl  side  chains.  It  was  nonetheless  helpful  in  determining  many  of  the  fundamental  physical 
properties  relating  to  structure  and  elastogenesis. 

2.1.2.2  Preparation  of  tropoelastin  by  microbial  biosynthesis  In  an  important 
development,  Weiss  and  coworkers  utilized  recombinant  DNA  technology  to  prepare  the 
precursor  protein,  tropoelastin,  in  E.  coli,  which  allowed  preparation  of  pure  protein  in  large 
quantities  (Vrhovski  et  al  1997;  Wu  et  al  1999).  Important  aspects  deduced  by  the  model 
systems  and  by  preparations  of  model  structures  derived  from  the  elastin  and  from  knowledge  of 
sequence  have  been  confirmed  and  extended  by  recombinantly  prepared  tropoelastin. 

2.1.2.3  Preparation  ofa-elastin  A  model  system  useful  in  stepping  systematically  from 
polymers  of  model  repeating  sequences  to  the  natural  fibrous  elastin  has  been  a  chemical 
fragmentation  product  from  fibrous  elastin  called  ot-elastin.  It  is  a  70,000  molecular  weight 
fragment  comprised  of  16  cross-linked  chains  (Partridge  et  al  1955;  Partridge  &  Davis  1955; 
Partridge  1966,  1967). 

2.1.3  Purification  using  phase  transitional  behavior 
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In  all  cases  these  protein-based  polymers  are  soluble  in  water  at  a  lower  temperature  and 
phase  separate  on  raising  the  temperature  above  that  for  the  onset  of  a  hydrophobic  folding  and 
assembly  transition.  This  phase  separation  process  has  been  called  coacervation.  It  is  a 
fundamental  property  of  these  polymers  with  the  correct  balance  of  apolar  and  more  polar 
residues. 

Cross-linking  of  Elastin  Models 

2.2  Cross-linking  of  elastin  models 

2.2.1  General  cross-linking  procedure 

Cobalt  60  y-irradiation  has  been  used  to  cross-link  the  phase  separated,  coacervated  and 
structured  state  of  the  elastic  protein-based  polymers  with  interpenetrating  polymer  chains.  The 
effective  cross-linking  doses  have  ranged  from  6  to  30  Mrads  (million  radiation  absorbed  dose) 
with  20  Mrads  chosen  as  the  usual  circumstance  for  the  elastic  model  proteins. 

2.2.2  Efforts  to  cross-link  a-elastin 

Even  though  a-elastin  undergoes  phase  separation,  it  does  not  form  an  elastic  matrix  on  y- 
irradiation.  As  a  cluster  of  some  16  cross-linked  chains  of  a  total  of  70  kDa  molecular  weight  it 
does  not  cross-link  because  there  are  no  chains  of  sufficient  length  to  interpenetrate. 

2.2.3  Efforts  to  cross-link  heat  denatured  poly(GVGVP) 

While  the  phase  separated  state  of  poly(GVGVP)  with  63%  water  by  weight  (Figure  1)  cross¬ 
links  very  effectively,  the  heat  denatured  state  obtained  on  prolonged  heating  above  60°C  does 
not.  Even  though  the  density  of  chains  is  twice  as  great  for  the  denatured  state  formed  above 
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60°C  (a  state  of  32%  water  by  weight  as  seen  in  Figure  1),  the  chains  no  longer  form  the 
intermediate  structured  state  of  Figure  1  and  the  chains  no  longer  sufficiently  interpenetrate  to 
result  in  an  elastic  matrix  on  y-irradiation.  This  is  simply  another  demonstration  that  the 
intermediate  between  30  and  60°C  is  a  structured  state  of  interpenetrating  chains  and  filaments. 

3.  Specialized  methodologies  and  analysis  of  results 

3.1  Atomic  Force  Microscopy  (AFM)  in  the  force-extension  mode 

3.1.1  Preparation  of  the  sample  for  AFM 

Glass  microscope  slides  were  cleaned  and  coated  with  about  30  nm  of  gold  in  a  home  built 
machine.  Si3N4  AFM  tips  were  used  (Microlevers,  Park  Scientific  Instruments,  Sunnyvale,  CA), 
and  also  coated  with  gold.  Quantities  of  1  mg  of  the  polypentapeptides  and  0.5  mg  methoxy- 
PEG-thiol  (Mw  5000)  were  dissolved  in  1  ml  Milli-Q  water  unless  otherwise  noted.  20  jd  of  this 
solution  were  incubated  on  the  gold  coated  slide  for  about  30  min  at  3°C  and  then  rinsed  with 
Milli-Q  water. 

3.1.2  The  AFM  Instrument 

A  detailed  description  of  the  AFM  force  spectroscopic  technique  and  details  of  the 
instrumental  set-up  have  been  given  elsewhere  (Oesterhelt  et  al  1999;  Clausen-Schaumann  et  al 
2000).  Briefly,  the  tip  of  a  cantilever  is  brought  into  contact  with  molecules  on  the  surface  and 
then  retracted,  its  deflection  -  and  therefore  the  force  -  are  detected  with  a  laser  by  optical  lever 
detection  (Cleveland  et  al  1993),  while  the  z-distance  is  controlled  by  a  strain  gauge  (Figure  2). 
The  nominal  spring  constants  of  cantilevers  used  in  the  experiments  were  10  mN/ra.  Prior  to  the 
first  approach  of  the  AFM  tip  to  the  surface,  the  spring  constants  of  each  lever  were  individually 
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calibrated  by  measuring  the  amplitude  of  its  thermal  oscillations  (Butt  &  Jaschke  1995).  The 
sensitivity  of  the  optical  lever  detection  was  measured  by  indenting  the  AFM  tip  into  a  hard 
surface.  All  experiments  were  conducted  in  Milli-Q  water  at  room  temperature  (21°C),  unless 
otherwise  specified. 

3.1.3  Obtaining  the  single-chain  force-extension  curve 

For  the  characterization  of  elastic  properties,  the  gold-coated  AFM  tip  was  brought  into 
contact  with  the  polypentapeptides  on  the  microscope  slide  by  manual  control.  The  AFM  tip  and 
the  polymer  layer  were  kept  in  contact  under  a  contact  force  of  several  nanoNewton  (nN)  for 
approximately  30  seconds  to  allow  for  an  attachment  of  the  chain  to  the  tip.  Usually  the  chain 
adsorbed  non-covalently  at  some  position  on  the  tip,  maybe  sometimes  a  chemical  reaction 
between  the  gold-coated  cantilever  and  the  cysteine  on  the  end  of  a  polypentapeptide  allowed  for 
really  high  rupture  forces.  Upon  retraction  of  the  cantilever,  individual  polypentapeptides  were 
stretched  between  the  surface  and  the  AFM  tip. 

Note  that  the  first  force-distance  profile  recorded  after  tip-substrate  contact  can  be  rather 
complex,  consisting  of  contributions  from  stretching  several  polypentapeptides,  desorption  from 
the  substrate  and/or  cantilever,  bond  rupture  of  short  strands,  as  well  as  interchain  aggregation 
and  entanglements.  Therefore,  in  each  measurement  the  cantilever  was  first  retracted  from  the 
substrate  to  a  distance  at  which  unspecific  adhesion  was  no  longer  observed.  Then,  in  successive 
retraction-approaching  cycles  the  distance  range  is  increased  while  trying  to  avoid  contact 
between  the  tip  and  additional  polypentapeptide  strands  at  the  substrate  surface,  until  only  one 
polymer  strand  remained  between  tip  and  substrate.  The  force-distance  profile  of  this  strand  was 
then  measured  repeatedly  until  rupture. 

3.1.4  Analysis  of  the  results 
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The  experimental  traces  were  fit  by  an  extended  wormlike  chain  (WLC)  model  including 
linear  elastic  contributions  arising  from  the  stretching  of  bond  angles  and  covalent  bonds: 

f.±jl=!L_L+ - l - -l  (ii) 

k„T  L  Ka  4(1  -RJL+FIKJ2  4 

In  this  expression,  Rz  is  the  measured  end-to-end-distance  at  any  given  force,  F,  and  L  is  the 
contour  length  of  the  stretched  chain  (polypentapeptide)  under  zero  force  (F  =  0).  The  peptides 
bending  rigidity  is  expressed  by  the  chain’s  persistence  length,  Lp.  Finally,  the  chain's 
extensibility  upon  stretching  is  described  by  the  segment  elasticity,  Kq,  which  is  introduced  into 
Eq.  (9)  as  a  linear  term  (hereby,  K0  can  be  understood  as  the  inverse  of  the  normalized 
compliance  of  a  Hookean  spring;  the  spring  constant  of  the  polymer  chain  is  given  by  Kq/L).  (It 
should  be  appreciated  that  the  WLC  model  is  not  modeled  at  the  molecular  structure  level  but 
rather  is  based  on  more  macroscopic  considerations.) 

This  resulted  in  values  for  the  persistence  length  Lp  of  (GVGVP)nX25i  :  Lp  =  0.4  nm  in  the 
low  force  regime  and  Lp  =  0.6  nm  when  fitted  to  higher  forces.  The  reversible  traces  provide 
evidence  for  single  chain  ideal  elasticity  (entropic  and  therefore  reversible  -  most  probably  from 
backbone  torsional  movements  (rocking),  wherein  the  contribution  to  the  entropic  component  of 
the  force  increases  with  decrease  in  the  frequency  of  the  oscillation).  The  small  deviations  from 
perfect  reversibility  in  traces  2,  3  and  5  of  Figure  3  occurred  when  the  chain  was  held  at  least  30 
s  in  the  relaxed  state.  These  deviations  could  have  several  origins.  There  could  be  a 
configurational  transition,  e.g.,  an  opening  of  P -turns  and/or  hydrophobic  interaction,  a 
backfolding,  between  different  parts  of  the  chain.  The  latter  is  favored  by  the  observed  rate 
dependence  of  this  non-equilibrium  contribution. 

The  persistence  length  Lp  of  (GVGIP)nx26o  is  difficult  to  determine,  when  intra  and/or  inter 
chain  aggregation  gives  rise  to  a  strong  hysteresis  even  when  being  far  away  from  the  surface 
within  two  successive  cycles  (Figure  4).  It  is  determined  to  exhibit  an  Lp  =  0.7  nm,  which  might 
be  a  little  higher  than  with  (GVGVP)nx25i,  but  this  value  lies  well  within  the  error  range  of  about 
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20  %.  This  strong  aggregation  of  (GVGIP)nx260  happens  at  the  same  concentration,  at  which  the 
aggregation  for  (GVGVP)nx25i  is  much  smaller.  This  reflects  the  increase  in  hydrophobicity  due 
to  an  additional  CH2  group  in  GVGIP,  which  is  the  only  difference  with  GVGVP.  Unfortunately, 
it  is  not  yet  clear,  whether  the  aggregation  is  mainly  due  to  backfolding  within  a  single  chain  or 
due  to  association  between  different  chains.  This  is  under  investigation. 

3.1.5  Comparison  of  AFM  single-chain  and  macroscopic  elastic  moduli 

For  the  characterization  of  elastomers  at  the  macroscopic  scale,  force  is  plotted  versus 

elongation,  but  instead  of  being  given  as  actual  increase  in  length  as  in  Figures  3  and  4, 

elongation  is  given  in  terms  of  the  ratio  of  the  initial  length  to  the  extended  length.  Furthermore, 

the  deformation  is  generally  uniform  so  that  the  initial  cross-sectional  area.  A,  is  used.  So  the 

Young’s  elastic  modulus,  Y,  is  the  force,  f,  times  the  initial  length,  Lj,  divided  by  the  elongation, 

2  2 

L,  and  the  cross  section  of  the  elastomer,  that  is,  Y  =  fLj  /AL.  The  units  are  N/m  or  dynes/cm 
where  N/m2  =  10  dynes/cm2  and  N  stands  for  Newtons.  The  elastic  modulus  of  20  Mrad  y- 
irradiation  cross-linked  (GVGVP)251 ,  i.e.,  X20-(GVGVP)251,  is  1.6  x  105  N/m2. 

In  order  to  utilize  the  data  in  Figure  3  to  obtain  an  estimate  of  single  chain  elastic  modulus, 
estimates  of  the  Lj  /L  ratio  and  of  A  are  required.  Initially,  the  total  number  of  repeats  being 
extended  must  be  deduced.  In  Figure  3,  the  length  to  rupture  is  just  greater  than  700  nm,  which 
we  take  to  be  the  fully  extended  chain.  This  is  too  long  for  a  chain  of  251  pentamers,  but  appears 
to  be  too  short  for  a  chain  of  753  pentamers.  By  elimination  of  chain  lengths  with  251  and  753 
pentamers,  this  indicates  that  the  structure  would  be  (GVGVP)2x25j  giving  502  pentamers,  i.e., 
2510  residues  with  a  value  of  0.29  nm  per  residue.  This  is  shorter  than  the  usual  value  of  0.35 
nm  per  residue  due  to  the  presence  of  a  proline  residue  in  every  pentamer.  This  is  considered  to 
be  due  to  the  cis  Val-Pro  peptide  bond,  which  introduces  a  kink  in  the  extended  backbone  for 
every  pentamer.  It  should  be  noted,  however,  that  if  all  of  the  P-tums  were  retained  at  rupture, 
the  3x251  chain  length  could  apply  with  approximately  1  nm  per  pentamer,  or  just  greater  than 
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700  nm.  This  was  not  considered  because  the  forces  for  detachment  can  be  greater  than  500  pN, 
and  retention  of  the  (5- turn  was  considered  less  likely  at  such  forces. 

Perhaps  the  most  significant  factor  in  obtaining  a  relevant  value  for  the  single  chain  elastic 
modulus  is  the  choice  of  the  structure  before  elongation  in  order  to  provide  an  estimate  of  an 
initial  length  and  cross  sectional  area.  In  this  regard  two  quite  different  initial  lengths  may  be 
considered,  one  for  the  random  coil  state  below  the  hydrophobic  folding  transition  and  another 
for  the  P-spiral  structure  that  could  only  occur  above  the  onset  temperature  for  the  hydrophobic 
folding  transition. 

3.1.5.1  Calculation  of  elastic  modulus  from  AFM  data  assuming  a  random  coil  state  for 

(GVGVP^si  At  temperatures  lower  than  the  onset  temperature  for  the  hydrophobic  folding 

and  assembly  transition,  the  more  correct  representation  of  the  unstretched  single  chain  would  be 

the  random  coil  state.  Based  on  photon  correlation  spectroscopic  studies,  the  hydrodynamic 

diameters  of  disordered  polymeric  coils  resulting  from  502  pentamers  (some  200,000  Da 

molecular  weight)  below  the  transition  would  be  about  18  nm  (San  Biagio  et  al  1988).  By  way 

of  example,  at  a  distance  of  378  nm  in  Figure  3,  the  Lj  /L  ratio  would  be  1/20.  With  a  slope  of 

0.1  pN/nm  from  Figure  3,  at  a  distance  of  378  nm  from  the  substrate  the  increase  in  force  would 

be  36  pN,  i.e.,  due  to  elongation  from  18  to  378  nm,  and  the  change  in  length,  AL,  would  be  a 

factor  of  20  or  a  slope  of  1.8  pN  for  a  100%  extension.  Next  this  quantity  should  be  divided  by 

an  appropriate  cross  sectional  area,  that  is,  n  times  the  hydrodynamic  radius  squared  of  the 

random  coil  prior  to  extension  of  the  single  chain,  i.e.,  nr  =  71(9  nm)  =  254  nm  .  Given  the 

4  2 

foregoing  assumptions,  the  calculated  single  chain  elastic  modulus  would  be  7  x  10  N/m  . 

Somewhat  surprisingly,  even  without  taking  cross-link  density  into  consideration  (Urry  et  al 

20 

1984),  this  value  is  only  a  factor  of  two  from  the  experimental  value  at  37°C  for  X  - 
(GVGVP)25i  of  1.6  x  105  N/m2. 
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3. 1.5.2  Calculation  of  elastic  modulus  from  AFM  data  assuming  a  single  chain  |3- spiral 

structure  for  (GVGVP)nx25 1  F°r  comparison  it  is  useful  to  consider  the  structure  representative 

of  the  state  above  the  temperature  of  the  inverse  temperature  for  hydrophobic  folding  and 

assembly.  If  the  single  P-spiral  structure  were  assumed  before  extension  with  three  pentamers 

per  turn  and  a  translation  along  the  axis  of  1  nm  per  turn,  the  length  of  a  502  mer  would  be  about 

170  nm.  This  means  that  an  elongation  of  200%  would  occur  at  510  nm.  With  this  information 

and  a  cross-sectional  area  of  the  P-spiral  structure  of  about  2.5  nm  ,  the  reversible  curves  in 

6  2 

Figure  3  with  a  slope  of  0.1  pN/nm  would  indicate  an  elastic  modulus  of  about  7  x  10  N/m  . 
This  is  a  factor  of  100  greater  than  that  obtained  under  the  assumption  of  a  random  coil.  It  will 
indeed  be  interesting  to  determine  the  elastic  modulus  of  the  fundamental  filament,  once 
synthetically  prepared,  which  is  thought  to  be  a  triple  stranded  twisted  filament,  as  depicted 
below  in  Figure  9F. 

3.2  The  acoustic  absorption  experiment 

3.2.1  Sample  preparation  and  experimental  set-up 

The  polymer  is  prepared  in  the  phase  separated  state,  i.e.,  the  intermediate  structured  state  of 
Figure  1,  and  20  Mrad  y-irradiation  cross-linked  in  the  form  of  the  desired  cylindrical  shape  and 
with  path  lengths  that  differ  by  a  factor  of  two,  as  shown  in  the  lower  part  of  Figure  5. 

The  Naval  Research  Laboratory  at  SSC-SD  has  the  capability  to  measure  very  accurately  the 
elastic  properties  of  materials  by  using  a  broadband  technique  developed  by  Sheiba  (1996).  The 
technique  is  under  US  Patent  Pending.  Older  techniques  are  more  cumbersome  and  less  reliable. 
The  frequency  of  measurement  ranges  from  sub-sonic  to  about  15  kHz.  The  measurements  are 
performed  in  air  inside  an  environmental  chamber.  Figure  5  also  shows  the  block  diagram  of  the 
instrumental  set-up  for  these  measurements. 
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The  network  analyzer  sends  sweeping  signals  to  the  shaker.  The  transducer  (sometimes 
called  an  accelerometer)  and  impedance  head  record  the  velocities  at  the  base  and  at  the  top  of 
the  sample.  After  recording  data  during  this  frequency  sweep,  a  new  set  of  data  is  then  obtained 
at  a  new  temperature  setting.  Commonly  temperatures  can  range  from  -30°  to  +70°C.  The 
velocity  record  constitutes  the  transfer  matrix.  By  mathematically  inverting  the  transfer  matrix, 
it  is  possible  to  extract  the  complex  modulus  and  Poisson  ratio.  These  data  provide  a  complete 
description  of  a  material’s  elastic  properties.  Further  analysis  yields  the  shear  modulus,  p,  the 
Poisson  ratio,  v,  and  the  loss  factor,  r|,  each  as  a  function  of  frequency  and  temperature.  The  loss 
factor,  T),  is  the  main  parameter  of  interest.  It  represents  the  energy  absorption  per  unit  volume  at 
a  given  frequency.  The  loss  factor  versus  frequency  and  temperature  were  found  for  several 
different  elastic  protein-based  polymer  compositions.  The  data  for  (GVGIP)nx260  is  given  in 
Figure  6A.  MATHCAD  7.0  was  used  to  perform  the  mathematical  analysis. 

3.2.2  Comparison  of  the  acoustic  absorptions  of  (GVGIP)  and  natural  rubber 

Natural  rubber  is  well  characterized  as  a  random  chain  network.  This  means  that  there  are  no 
regular  repeating  conformational  features  in  the  polymer  chain.  The  consequence  of  a  random 
chain  is  a  broad,  slowly  changing  dependence  on  frequency  of  physical  properties,  such  as 
absorption.  Simply  stated  the  absence  of  structure  means  that  the  barrier  to  rotation  about  each 
bond  will  tend  to  be  different.  This  is  what  is  observed  in  Figure  6A  for  natural  rubber  and  a 
synthetic  petroleum-based  polymer,  polyurethane,  chosen  for  its  relatively  high  loss  factor  in  the 
frequency  range  reported. 

On  the  other  hand  at  frequencies  lower  than  those  of  the  infrared  spectroscopic  range,  any 
material  that  exhibits  a  localized  absorption  band  as  a  function  of  frequency  must  do  so  by  virtue 
of  a  regular  repeating  structure.  Interestingly,  this  is  seen  in  Figure  6A  for  X  -(GVGIP)26o in 
the  range  of  200  Hz  to  7  kHz,  and  it  is  seen  that  the  intensity  of  the  absorption  per  unit  volume 
increases  as  the  temperature  is  raised  through  the  temperature  interval  of  the  inverse  temperature 
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transition.  Besides  demonstrating  the  remarkable  acoustic  absorption  property  of  this  material, 
the  acoustic  absorption  data  provide  additional  demonstration  of  the  existence  of  a  regular 
structure  for  poly(GVGIP)  and  for  related  polymers.  Actually,  this  is  simply  further 
confirmation  of  earlier  dielectric  relaxation  data,  which  demonstrated  the  development  of  intense 
localized  relaxations  near  5  MHz  as  the  temperature  is  raised  through  the  interval  of  the  inverse 
temperature  transition,  as  reviewed  below. 

3.3  Dielectric  Relaxation  Studies  on  Elastin-Related  Systems 

3.3.1  Dielectric  relaxation  and  acoustic  absorption  data  for  (GVGIP)„  over  the 
frequency  range,  0.1  to  100  kHz. 

In  the  acoustic  absorption  measurement  of  Figure  6A,  an  oscillating  compressional  wave 
activates  the  mechanical  relaxation  in  the  vicinity  of  one  kHz,  and  the  relaxation  requires 
reorientation  of  peptide  moieties  in  order  for  the  polypentapeptide  to  respond  to  the  acoustical 
frequency  of  the  vibrational  wave.  As  peptide  moieties  have  very  large  dipole  moments, 
poly(GVGIP)  could  be  expected  to  function  as  a  transducer  in  which  the  acoustically  activated 
mechanical  relaxation  would  be  transduced  into  a  dielectric  relaxation  signal  of  the  same 
frequency.  The  experimental  equipment  to  demonstrate  this  free  energy  transduction  has  yet  to 
be  designed. 

The  existence  of  the  acoustic  absorption,  however,  indicates  that  a  dielectric  relaxation 
should  be  observed  in  a  not  too  distant  frequency  range.  An  exact  overlap  is  not  expected  in  the 
two  different  experiments  because  in  the  dielectric  relaxation  experiment  that  alternates  positive 
and  negative  charge  on  a  pair  of  plates  directly  drives  the  peptide  reorientation,  and  the  particular 
repeating  structure  determines  the  frequency  range  and  intensity  of  the  relaxation.  On  the  other 
hand,  in  the  acoustically  driven  relaxation  the  peptide  reorientations  are  a  secondary  result  that  is 
required  in  order  for  the  repeating  structure  to  respond  to  the  hammering  frequency  of  the 
acoustic  wave.  Of  course,  the  fact  that  a  mechanical  resonance  grows  in  intensity  as  the  inverse 
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temperature  transition  progresses  requires  progressive  formation  of  a  regular  repeating  structural 
feature  on  raising  the  temperature. 

It  is  to  be  expected,  therefore,  that  a  dielectric  relaxation  be  observed  in  the  frequency  range 
of  0.1  to  100  kHz  as  reported  in  Figure  6B,  which  plots  the  imaginary  (absorption)  part  of  the 
dielectric  permittivity  of  poly(GVGIP)  as  a  function  of  temperature.  Just  as  for  the  acoustic 
absorption  study  in  Figure  6A,  there  is  an  increase  in  the  absorption  (imaginary)  component  of 
the  dielectric  relaxation  as  the  temperature  is  raised  through  the  temperature  interval  of  the 
inverse  temperature  transition  for  hydrophobic  folding  and  assembly.  Interestingly  intensity  of 
the  relaxation  continues  to  grow  above  20  to  40°C  suggesting  that  an  acoustic  absorption  even 
more  intense  than  reported  in  Figure  6A  might  be  expected. 

3.3.2  A  related  dielectric  relaxation  (mechanical  resonance)  near  5  MHz  for 
poly(GVGIP). 

For  poly(GVGIP)  in  the  frequency  range  of  1  MHz  to  1000  MHz  and  over  the  temperature 
range  of  0  to  60°C,  the  real  part  of  the  dielectric  permittivity  is  given  in  Figure  7A  and  the 
imaginary  (absorption)  part  in  Figure  7B.  As  with  the  data  in  the  1  kHz  range  of  Figure  6,  there 
is  an  increase  in  regular  structure  as  the  temperature  is  raised  from  below  to  above  the  inverse 
temperature  transition.  This  type  of  data  was  used  more  than  a  decade  ago  to  demonstrate  the 
formation  of  regular  structure;  it  is  equally  valid  today,  and  it  is  what  allowed  anticipation  of  the 
data  recently  obtained  in  Figure  6,  that  is,  anticipation  of  the  remarkable  acoustic  absorption 
properties  of  these  materials.  This  successful  trail  of  proposed  structural  regularity,  testing  of  the 
proposal,  and  predicting  of  additional  properties  continues  unbroken  for  much  of  two  decades. 

The  predictive  success  and  usefulness  of  a  perspective  is  the  truest  test  of  its  validity.  The 
perspectives  of  disorder  and  randomness  have  given  rise  to  no  such  progress.  On  the  other  hand 
the  perspective  of  a  regular  dynamic  structure  has  led  to  properties  such  as  numerous  free  energy 
transductions  otherwise  unimaginable.  Now,  the  data  on  the  elastin  models  can  be  directly 
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compared  to  a-elastin  and  fibrous  elastin,  itself.  Mechanical  resonances  due  to  regularly 
repeating  structures  are  observed  both  in  the  vicinities  of  1  kHz  and  of  5  MHz. 

3.3.3  Comparison  of  poly(GVGVP)  and  a-elastin 

Plotted  in  Figure  8A  is  the  real  part  of  the  dielectric  permittivity  as  a  function  of  temperature 
for  a-elastin  and  for  poly(GVGVP)  in  the  inset  (Urry  et  al,  1985a).  In  both  cases,  there  is  the 
development  of  an  intense  relaxation  near  5  MHz  as  the  protein-based  elastomers 
hydrophobically  fold  and  assemble  on  increasing  the  temperature  through  the  range  of  the 
inverse  temperature  transition.  Both  elastomeric  molecular  systems  exhibit  the  same 
development  of  regular  structure  on  raising  the  temperature  through  the  range  of  the  inverse 
temperature  transition. 

3.3.4  Comparison  of  cross-linked  poly(GVGVP)  to  fibrous  elastin 

In  Figure  8B  is  plotted  the  imaginary  (absorption)  part  of  the  dielectric  relaxation  for  purified 
natural  fibrous  elastin  (from  bovine  ligamentum  nuchae)  as  a  function  of  temperature  from  23  to 
58°C  (Luan  et  al  1988).  Also,  in  Figure  8B  is  plotted  the  same  data  at  58°C  for  the  elastomer,  20 
Mrad  cross-linked  poly(GVGVP).  While  there  is  a  small  increase  in  mean  frequency  for  the 
X20-(GVGVP)25i  relaxation  in  Figure  8B  as  in  Figure  8A,  it  is  unmistakable  that  there  is 
hydrophobic  folding  and  assembly  resulting  in  the  formation  of  regular,  dynamic,  non-random 
structures  in  the  natural  fibrous  elastin  that  are  similar  to  those  of  cross-linked  poly(GVGVP) 
and  poly(GVGIP).  Incidentally,  the  decrease  in  intensity  at  1  GHz  for  fibrous  elastin  in  Figure 
8A  for  a-elastin  and  in  Figure  8B  for  fibrous  elastin  reflects  the  loss  of  hydrophobic  hydration  as 
the  hydrophobic  folding  proceeds.  This  is  due  to  a  dielectric  relaxation  near  5  GHz,  which  arises 
from  hydrophobic  hydration  (Urry  et  al  1997). 
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To  be  brief,  the  experimental  data  clearly  reflect  dynamic  regular  structures  in  the  elastic 
protein-based  polymers  with  related  structures  and  corresponding  behaviours  in  the  natural 
elastin  protein. 

4.  DISCUSSION 

4.1  Relevance  of  proposed  mechanisms  for  protein  entropic  (ideal)  elasticity 

The  kernel  for  the  classical  theory  of  rubber  elasticity  originated  with  K.  H.  Meyer  in  1932, 
which  according  to  Flory  (1968)  can  be  stated  as  "...  the  orientation  of  the  molecules  of  apiece  of 
rubber  when  it  is  stretched  would  entail  a  decrease  in  entropy,  and  thus  could  account  for  the 
retractive  force. "  A  common  restatement  has  been  that  a  collection  of  disordered  chains 
becomes  more  aligned  as  the  elastomer  is  stretched,  and  that  this  alignment  of  chains  constitutes 
the  decrease  in  entropy  on  extension.  Flory’ s  restatement  in  1968  was  given  in  terms  of  two 
more  fundamental  subsidiary  hypotheses:  (1)  The  restoring  force  originates  within  the  chain 
molecules  by  distortion  of  the  distribution  of  their  configurations  and  not  from  interactions 
between  chains.  (2)  The  restoring  force  is  due  "exclusively  to  the  entropy  of  the  chain 
configuration."  For  ideal  or  perfect  elastomers,  we  have  no  quarrel  with  these  subsidiary 
hypotheses,  but  those  elastomers  exhibiting  hysteresis  appear  not  to  be  consistent  with  the  first 
hypothesis. 

Obviously  with  respect  to  protein  elasticity,  we  do  disagree  with  the  derivative  dogma  (Flory, 
1968;  Hoeve  and  Flory,  1958,  1974)  that  only  random  chain  networks  comprised  of  a  Gaussian 
distribution  of  end-to-end  chain  lengths  can  satisfy  the  two  hypotheses,  which  has  been  a  tenet  of 
the  classical  theory  of  rubber  elasticity.  Clearly,  entropic  elasticity  exhibited  by  a  single  chain 
limits  the  source  of  entropy  change  to  the  damping  of  internal  chain  dynamics  within  a  single 
chain  and  can  not  be  due  to  the  relative  orientation  of  many  chain  molecules  within  a  random 
chain  network. 
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4.1.1  The  Flory  random  chain  network  theory  of  entropic  elasticity? 


Flory’s  publications  describe  a  Gaussian  distribution  of  end-to-end  chain  lengths  that  results 
from  a  random  chain  network,  and  they  assert  that  the  decrease  in  entropy  due  to  extension  or 
compression  arises  from  distortion  of  the  network  from  this  most  probable  Gaussian  distribution 
(Flory  1953,  1968).  Hoeve  and  Flory  (1958,  1974)  concluded  that  the  entropic  elasticity  of 
elastin  arises  from  such  a  random  chain  network.  There  are  four  principal  reasons  why  this 
perspective  is  to  be  questioned.  These  are  (1)  a  limited  potential  for  cross-links  in  the  fixed 
sequence  of  elastin,  (2)  the  above  spectroscopic  evidence  for  regular,  non-random,  dynamic 
structural  features  in  components  of  elastin,  (3)  the  above  report  that  single  chains  of  repeating 
sequences  of  elastin  exhibit  entropic  elasticity,  and  (4)  irreversible  loss  of  elastic  force  due  to 
heating  above  60°C  due  to  thermal  denaturation. 

4.1.1.1  Relevance  of  the  probability  of  cross-link  formation  to  structure:  In  rubber,  each 
repeating  unit  is  a  potential  cross-link,  and,  on  vulcanization,  any  pair  of  repeating  units  can  be 
made  into  a  cross-link  wherever  two  chain  segments  are  sufficiently  proximal.  In  elastin,  only 
lysine  residues  form  the  cross-links,  and  yet  there  are  only  some  forty  lysine  residues  for  an 
elastin  protein  of  about  800  residues,  i.e.,  one  lysine  in  every  twenty  residues.  Furthermore,  the 
lysine  residues  are  arranged  in  an  orderly  manner  along  the  sequence.  Most  commonly,  four 
lysine  residues  combine  to  form  a  single  cross-link,  called  a  desmosine,  and  essentially  all  of  the 
lysine  residues  in  fibrous  elastin  have  been  formed  into  cross-links.  Such  cross-linking  cannot 
be  achieved  within  a  random  chain  network.  Achievement  of  such  rare  cross-linkages  require 
elements  of  order  (Urry,  1976). 

4.1.1.2  Low  frequency  motions  provide  much  entropy  to  chain  structures:  Motions  due 
to  rotations  about  backbone  bonds  in  a  chain  molecule  that  exhibit  intense  relaxations  limited  to 
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localized  frequency  intervals,  that  is,  intense  mechanical  resonances,  can  only  result  from 
polymers  containing  regular,  non-random,  albeit  dynamic  structures.  For  model  elastin  systems 
including  purified  fibrous  elastin,  such  mechanical  relaxations  are  observed  near  1  kHz  and  5 
MHz  (see  Figures  6, 7  and  8,  Buchet,  et  al,  1988;  Urry  et  al,  1985a,  Luan  et  al,  1988).  As  will  be 
treated  below,  these  relaxations  provide  an  ample  source  of  entropy  to  structured  elastin  systems. 
These  contributions  to  entropy  increase  greatly  as  the  frequency  of  the  relaxation  is  lower,  and 
they  are  contributions  to  the  free  energy  of  protein  structure  that  have  yet  to  be  adequately 
treated  in  the  current  computations  of  low  energy  protein  structures.  Most  importantly  for 
entropic  elasticity,  low  frequency  motions  provide  an  abundant  and  ready  source  of  entropy 
decrease  on  extension. 

4.1.1.3  Entropic  elastomeric  force  arising  from  a  single  chain:  At  a  fixed  extension  in 
the  AFM  single-chain  force-extension  curves  of  Figure  3,  there  is  only  one  end-to-end  chain 
length.  The  entropic  elastic  force  measured  at  a  given  extension  cannot  have  arisen  from 
displacement  of  a  network  of  chains  from  a  Gaussian  distribution  of  end-to-end  chain  lengths. 
The  decrease  in  entropy  that  provides  the  entropic  restoring  force  in  the  single  chain  experiment 
must  arise  from  a  damping  of  internal  chain  dynamics  of  the  single  chain.  There  are  no  other 
degrees  of  freedom  from  which  the  decrease  in  chain  entropy  can  arise. 

4.1.1.4  Loss  of  elastomeric  force  of  elastin-based  systems  due  to  thermal  denaturation:  In 
Figure  1  above,  it  was  noted  that  poly(GVGVP)  exhibits  a  slow  irreversible  denaturation  when 
held  at  temperatures  above  60°C.  This  same  phenomenon  is  apparent  in  the  thermoelasticity 
studies  on  fibrous  elastin,  as  originally  noted  by  Dorrington  and  McCrum  (1977).  In  the 
thermoelasticity  experiments  of  Figure  9  (Urry  1988a,  1988b),  the  elastomer  is  equilibrated  at 
40°C  and  at  that  temperature  is  extended  to  60%.  At  this  fixed  length,  the  temperature  is 
lowered  to  about  20°C  and  slowly  raised  to  80°C.  The  force  dramatically  increases  as  the 
temperature  is  raised  through  the  range  of  the  inverse  temperature  transition  from  20  to  40°C. 
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Then,  the  force/T(°K)  tends  to  remain  nearly  constant  from  40  to  60°C,  especially  for 
X20-poly(GVGVP).  The  value  of  lnf/T  begins  to  decrease  above  60°C  in  a  time  dependent 
manner.  When  using  a  rate  of  increase  of  4°C  per  30  min,  8°C/hr,  there  is  a  decrease  in  lnf/T, 
and,  when  the  temperature  is  returned  to  40°C,  there  has  occurred  an  irreversible  loss  of  force. 
The  decrease  in  force  above  60°C  is  more  dramatic  on  raising  the  temperature  more  slowly,  e.g., 
when  raising  the  temperature  at  a  rate  of  l°C/hr.  Irreversible  thermal  denaturation  is  most 
dramatically  seen  when  kept  above  60°C  for  a  longer  time.  Furthermore,  a  plot  of  ln(elastic 
modulus)  versus  time  at  80°C  is  linear  in  evidence  of  a  first  order  process  (Urry  1988a,  1988b). 
Quite  interestingly  while  heating  above  60°C  decreases  force,  it  also  shortens  the  length  at  which 
zero  force  is  obtained,  as  to  be  expected  for  a  filamentous,  P-spiral-like  structure,  that  became 
randomized  by  thermal  denaturation. 

As  discussed  in  relation  to  Eqns.  (7)  and  (8)  above,  a  slope  of  zero  for  plots  like  those  in 
Figure  9  would  mean  an  ideal  (entirely  entropic)  elastomer.  As  clearly  seen  for 
X20-poly(GVGVP)  in  the  40  to  60°C  temperature  range  when  the  data  is  collected  at  the  rate  of 
30  min  per  data  point,  a  near  zero  slope  is  obtained.  Commonly  an  fE/f  ratio  of  0.1  is  obtained 
for  X20-poly(GVGVP)  (Luan  et  al  1989).  Because  of  the  steep  positive  slope  between  20  and 
40°C,  due  to  the  hydrophobic  folding  and  assembly  of  the  inverse  temperature  transition,  and 
because  of  the  slow  irreversible  thermal  denaturation  above  60°C,  the  slope,  obtained  in  the 
thermoelasticty  experiment  from  which  the  entropic  contribution  is  estimated,  depends  on  the 
rate  at  which  the  data  is  collected  and  the  temperature  interval  over  which  it  is  collected. 

A  steep  rise  in  force  in  the  20  to  40°C  temperature  interval  for  fibrous  elastin  in  Figure  9A, 
though  less  steep  than  seen  in  Figure  9B  for  X  -poly(GVGVP),  demonstrates  the  effect  of  the 
inverse  temperature  transition  within  elastin  itself.  More  dramatic  for  fibrous  elastin  is  the  loss 
of  elastomeric  force  on  raising  the  temperature  above  60°C.  This  finding  of  a  slow  loss  of  force 
at  elevated  temperature  would  not  have  been  so  apparent  in  the  limited  time  of  a  single  day  to 
perform  the  experiment  (without  the  advantage  of  computerized  equipment),  as  was  the  case 
when  Hoeve  et  al  (1958)  carried  out  their  original  study  on  fibrous  elastin.  For  the 
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thermoelasticity  study  of  fibrous  elastin  shown  in  Figure  9,  however,  the  rate  of  loss  of  force  at 
elevated  temperature  is  even  greater  than  for  X  -poly(GVGVP). 

Additionally,  the  plot  of  ln(elastic  modulus)  versus  time  for  fibrous  elastin  maintained  at 
80°C,  as  determined  by  a  stress-strain  curve  every  24  hrs,  is  a  straight  line  with  a  half-life  of 
approximately  10  days  (Urry  1988b).  Another  important  feature  of  the  data  for  determining  half- 
life  is  that  the  length  at  which  zero  force  is  obtained  decreases  with  time  at  80°C.  Remarkably, 
this  means  that  the  loss  of  force  in  Figure  9  and  the  progressive  loss  of  elastic  modulus  exhibited 
by  a  series  of  stress-strain  curves  at  80°C  is  not  due  to  structural  rearrangement  to  an  extended 
state  under  load  or  for  that  matter  due  to  slipping  in  the  grip.  Instead,  extended  periods  of  time  at 
high  temperature  results  in  a  denaturation,  which  entails  a  randomization  of  linear  filamentous 
structures  to  denatured,  more  nearly,  globular  structures.  This  has  analogy  to  the  loss  of  percent 
water  by  weight  for  time  above  60°C  in  Figure  1. 

What  is  clearly  being  observed  in  the  thermoelasticity  data  of  elastin-based  systems  in  Figure 
9,  in  the  half-life  studies  for  the  elastic  modulus  at  80°C,  and  in  the  structural  change  above  60°C 
in  Figure  1  is  a  slow  thermal  denaturation  of  protein-based  elastomers.  Since  random  chain 
networks  do  not  exhibit  thermal  denaturation,  these  elastin-based  systems  are  not  properly 
described  as  random  chain  networks.  Thus,  the  evidence  is  clear  and  abundant;  elastin-based 
systems  contain  regular  repeating  conformations  and  those  ordered  structures  denature. 

4.1.2  Solvent  entropy  changes  as  a  source  of  entropic  elastic  force? 

4.1.2.1  Effect  of  removal  of  solvent  entropy  change  during  inverse  temperature 
transition:  The  change  in  solvent  entropy  is  fundamental  to  the  occurrence  of  an  inverse 
temperature  transition.  This  transition  has  been  extensively  utilized  to  perform  the  family  of 
energy  conversions  that  sustain  living  organisms  (Urry  1993,  1997),  and  therefore  the  inverse 
temperature  transition  is  extensively  characterized.  One  approach  to  testing  whether  the  change 
in  solvent  entropy  contributes  to  the  elastomeric  force  is  to  utilize  an  approach  due  to  Flory,  that 
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is  to  reduce  the  heat  of  the  inverse  temperature  transition  to  near  zero  by  an  appropriate  solvent 
mixture.  Hoeve  and  Flory  (1958)  had  used  a  mixture  of  30%  ethylene  glycol  by  weight  in  water 
in  their  thermoelasticity  study  that  deduced  a  dominant  entropic  elastic  force  for  fibrous  elastin. 
When  Luan  et  al  (1989)  compared  the  thermoelasticity-derived  fE/f  ratio  of  X  -poly(GVGVP) 
in  water  and  in  30%  ethylene  glycol,  the  values  were  the  same,  0.1.  Based  on  the 
approximations  of  Eqns  (7)  and  (8),  X20-poly(GVGVP)  is  a  dominantly  entropic  elastomer  in 
both  solvent  systems.  The  effect  of  adding  the  ethylene  glycol  is  to  reduce  the  heat  of  the 
transition  to  near  zero.  Now,  since  the  entropy  of  the  desolvation  process  is  the  heat  of  the 
transition  divided  by  the  temperature  of  the  transition,  AHt/Tt,  the  solvent  entropy  change  would 
be  approaching  zero,  and  yet  the  solvent  mixture,  instead  of  reducing  the  elastomeric  force, 
remarkably  resulted  in  a  substantial  increase  in  force.  The  conclusion  is  that  there  is  no 
experimental  basis  for  believing  that  solvent  entropy  change  contributed  to  elastomeric  force. 

4.1.2.2  The  polymer  chains  bear  the  force;  solvation  changes  alter  the  fraction  of  chain 
segments  that  sustain  the  force :  The  force  of  deformation  must  be  borne  by  the  chains;  it  cannot 
be  borne  by  solvent.  At  fixed  length,  raising  the  temperature  of  a  band  of  elastomer  from  below 
to  above  the  temperature  of  the  inverse  temperature  transition  increases  the  force.  Alternatively, 
increasing  chemical  potential  to  lower  the  transition  temperature  from  above  to  below  the 
operating  temperature  at  fixed  length  causes  an  increase  in  elastomeric  force  (Urry  et  al  1997, 
1988a,  1988b).  These  results  are  due  to  an  increase  in  hydrophobic  folding  and  assembly  of  the 
chains.  Since  this  is  the  case  at  fixed  length,  the  fraction  of  a  chain  length  that  is  not  folded  must 
be  further  extended.  Accordingly,  in  order  to  maintain  a  fixed  length  the  unfolded  segments 
must  be  further  extended.  This  provides  for  the  increase  in  entropic  elastic  force  by  means  of  a 
decrease  in  configurational  entropy  and  not  directly  due  to  a  change  in  solvent  entropy.. 

4.1.3  Entropic  elasticity  due  to  damping  of  internal  chain  dynamic  on  extension 
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4.1.3.1  Molecular  structure  of  poly(GVGVP):  In  order  to  understand  the  concept  of  the 
damping  of  internal  chain  dynamics  on  extension,  the  structure  given  in  Figure  10  is  helpful.  The 
Pro-Gly  containing  elastomeric  sequential  polypeptides  exhibit  a  repeating  P-tum  structure  as 
schematically  shown  in  Figure  10A,  and  the  P-tum  is  given  in  crystallographic  detail  in  Figure 
10B.  On  raising  the  temperature  the  inverse  temperature  transition  results  in  the  optimization  of 
intramolecular  and  intermolecular  hydrophobic  association.  Intramolecularly,  this  results  in  a 
helical  structure  called  a  P-spiral,  which  is  shown  schematically  as  a  helix  in  Figure  10C,  in  more 
detail  with  P-tums  schematically  shown  as  spacers  between  turns  of  the  helix  in  Figure  10D,  and 
in  atomic  stereo  detail  in  Figure  10E.  Intramolecularly,  the  P-spirals  associate  to  form  multi- 
stranded  twisted  filaments,  as  depicted  in  Figure  10F.  As  there  is  water  within  the  p-spiral  and 
as  a  result  of  the  P-tums  functioning  as  spacers,  the  peptide  segments  connecting  P-tums  are  free 
to  undergo  large  amplitude  torsional  oscillations. 


4.1.3.2  Unrolled  perspective  of  P -spiral  showing  suspended  segment  between  p -turns: 
Shown  in  Figure  11  is  an  unrolled  perspective  of  the  polypentapeptide  P-spiral  showing  a 
suspended  segment  between  P-tums  where  peptide  moieties  are  free  to  undergo  large  amplitude 
torsional  oscillations.  In  particular,  the  y4-<t>5  pair  of  torsion  angles  and  the  \|/5-<|)j  pair  of  torsion 
angles  couple  in  each  case  to  result  in  rocking  motions  of  the  respective  peptides,  i.e.,  in  peptide 
librations.  The  point  to  be  made  is  that  the  peptide  moieties  can  undergo  large  torsional 
oscillations  due  to  the  absence  of  steric  interactions  in  the  hydrophobically  folded  P-spiral 
structure.  This  freedom  of  motion  of  the  peptide  moieties  of  the  suspended  segment  within  the 
P-spiral  structure  gives  much  entropy  to  the  structure  and  is  considered  to  be  a  primary  source  of 
the  mechanical  resonances  of  Figures  6, 7  and  8,  as  will  be  treated  below. 

4.1.3.3  Effect  of  extension  on  the  amplitude  of  the  torsional  oscillations  and  on  chain 
entropy.  The  magnitude  of  the  torsional  (dihedral  angle)  excursions  for  a  single  pentamer 
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sandwiched  within  a  turn  of  (3-spiral  (three  pentamers)  for  a  1.5  kcal/mol-residue  cut-off  energy 
is  shown  in  Figure  12A  when  using  the  ECEPP  molecular  mechanics  calculations  of  the 
Scheraga  approach.  Dihedral  angle  oscillations  of  160°  are  possible.  Lambda  oscillational  plots 
of  \(/j— i  for  the  two  peptide  moieties  of  Figure  1 1  are  given  in  Figure  13  for  the  relaxed  and 
extended  (130%)  states  of  poly(GVGVP)  where  the  correlated  motions  of  the  pair  of  angles  is 
apparent  by  the  data  points  falling  on  a  line  with  a  slope  of-1. 

The  effect  of  a  130%  extension  of  the  (3-spiral  pitch,  shown  structurally  in  Figure  12B  and  in 
the  plots  of  Figure  13,  is  a  dramatic  decrease  in  the  amplitude  of  the  torsional  oscillations 
allowed  within  a  1.5  kcal/mol-residue  cut-off  energy.  Using  an  enumeration  of  states  approach 
with  each  change  of  5°  in  any  torsion  angle  counted  as  a  new  state,  the  decrease  in  numbers  of 
states  for  a  0.6  kcal/mol-residue  cut-off  energy  (i.e.,  within  an  RT)  is  from  a  Wr  of  342  to  a  We 
of  24.  Using  Eqn  (10),  this  gives  for  AS  a  value  of  -5.28  cal/mol-deg/pentamer.  If  the 
Boltzmann  summation  over  all  energies  is  used,  the  value  is  -5.06  cal/mol-deg/pentamer  (Urry  et 
al  1985c). 

The  plots  of  Figure  13C  show  similar  reductions  in  freedom  of  motion  to  the  130%  extension 
on  replacement  of  the  Gly5-residue  by  a  D-Ala5-residue.  An  expected  reduction  in  the  intensity 
of  the  5  MHz  relaxation  has  been  observed  and  both  are  calculable  using  the  peptide  dipole 
reorientations  of  the  pentapeptide  moiety  in  the  (3-spiral  structure  (Venkatachalam  and  Urry, 
1986).  Furthermore,  replacement  of  the  Gly5-residue  by  an  L-Ala5-residue  resulted  in  complete 
loss  of  elasticity  (Urry  et  al,  1983). 

Alternatively,  when  using  the  molecular  dynamics  approach  of  Karplus,  the  expression 
becomes, 

AS  =  R  lntrijA^j6  Ay  ,e  /TljA^  Ay  (11) 

where  the  product  TTj  is  over  the  rms  (root  mean  square)  torsion  angle  fluctuations  A^6  and  Ay  j6 
for  the  extended  state  divided  by  the  product  of  the  rms  torsion  angle  fluctuations  A<t>jr  and  Ay  ,r 
for  the  relaxed  state  as  obtained  from  the  molecular  dynamics  trajectories  for  the  relaxed  and 
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130%  extended  states.  The  set  of  rms  torsion  angle  fluctuations  for  the  relaxed  and  extended 
states  are  given  in  Table  1,  where  the  decreases  in  amplitude  of  the  torsional  oscillations  as  the 
result  of  extension  are  apparent.  Using  Eqn  (11),  the  calculated  decrease  in  entropy  is  -5.5 
cal/mol-deg/pentamer  for  the  130%  extension  (Chang  and  Urry,  1989).  Whether  the  calculation 
of  entropy  is  by  means  of  an  enumeration  of  states  in  configuration  space  or  the  amplitude  of  the 
torsional  oscillations  in  a  molecular  dynamic  simulation,  the  calculated  entropy  change  for  a 
130%  extension  is  the  same. 

The  A<|>-A\j/  representation  of  entropy  change  on  extension  of  Eqn  (11)  provides  an 
opportunity  for  a  particularly  insightful  sense  of  entropy.  Figure  14  plots  in  three  dimensions  the 
magnitude  of  the  A(J),  A\|/  values  for  a  chain  segment  of  three  torsion  angles.  As  such,  the  volume 
is  recognized  by  Eqn  (11)  as  being  proportional  to  the  configurational  entropy  of  the  chain 
segment.  The  decrease  in  magnitude  of  the  torsion  angle  oscillations  on  extension  is  seen  as  a 
decrease  in  volume,  which  is  proportional  to  the  decrease  in  the  configurational  entropy  of  the 
chain  segment  due  to  extension.  Volume  in  configuration  space  provides  a  representation  of 
chain  entropy  and  a  visualization  of  the  decrease  in  entropy  on  extension.  Obviously  for  a  single 
chain  of  N-residues  fixed  at  its  ends,  entropy  of  the  chain  can  be  represented  as  a  volume  in  2N- 
dimensional  configurational  space,  and  extension  decreases  the  volume  in  N-dimensional 
configuration  space,  in  a  manner  represented  in  Figure  14  for  three  dimensions.  In  particular,  the 
values  in  Table  1  are  actually  the  rms  A<|>  and  A\|/  for  the  relaxed  state  and  the  extended  state. 
Especially  for  the  suspended  segment  the  decrease  in  rms  A<t>  and  A\| /  values  are  substantial  on 
extension,  leading  to  the  calculated  decrease  in  entropy  on  extension. 

4.1.3.3  Calculation  of  contribution  of  force  due  to  decrease  in  freedom  of  motion  on 
extension :  From  Eqns  (4)  and  (5)  we  have  that 

fs  =  -  T(0S/3L)vt.  (12) 

At  physiological  temperatures,  310K,  and  for  an  extension  from  .35  to  .8  nm  for  a  pentamer 
within  the  P-spiral  structure  of  Figure  10E  with  a  calculated  decrease  in  entropy  of  5.5  cal/mol- 
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deg/pentamer,  the  calculated  value  for  the  entropic  contribution  to  the  elastomeric  force  is  24  pN 
for  a  single  chain  of  P-spiral. 

2 

To  put  this  in  terms  of  an  elastic  modulus  in  units  of  N/m  ,  it  is  necessary  to  estimate  the 

number  of  P-spirals  per  m  .  From  Figure  1  the  functional  state  is  37%  peptide  by  weight. 

3  17  2 

Taking  the  density  of  the  peptide  to  be  1.3  gm/cm  ,  there  would  be  some  10  (J-spirals/m  .  The 
result  would  be  24  x  10  N/m  .  The  calculated  result  is  an  order  of  magnitude  higher  than  the 
experimental  value  of  1.6  x  105  N/m2  for  20  Mrad  cross-linked  (GVGVP)25i.  Several  factors 
would  be  expected  to  result  in  a  lower  experimental  elastic  modulus  for  the  matrix.  Firstly,  in 
the  cross-linked  matrix  the  P-spirals  would  be  randomly  oriented  with  respect  to  direction  of 
extension.  Secondly,  the  P-spirals  occur  as  part  of  a  multi-stranded  twisted  filament  as  shown  in 
Figure  10F,  which  would  somewhat  restrict  the  freedom  of  motion  of  the  torsional  oscillations 
and  result  in  a  lower  entropy  change  on  extension.  What  has  been  demonstrated,  however,  is 
that  the  damping  of  internal  chain  dynamics  on  extension  provides  an  abundant  source  of 
entropic  restoring  force  in  elastin-based  systems.  Furthermore,  given  the  entropic  elasticity  of 
single  chains  demonstrated  by  the  AFM  single-chain  force-extension  curves  for  many  different 
polymers,  damping  of  internal  chain  dynamics  on  extension  would  seem  to  be  a  fundamental 
source  for  the  development  of  entropic  elastic  force  in  chain  molecules. 

4.2  Dependence  of  entropy  and  structural  free  energy  on  oscillator  frequency 

One  final  point  should  be  noted  in  this  report.  It  is  the  significance  of  the  occurrence  of 
intense  low  frequency  mechanical  motions  in  elastin-based  systems.  The  insight  comes  from  the 
expression  for  the  dependence  of  entropy  on  frequency  of  a  harmonic  oscillator  as  given  in 
Figure  15.  From  the  expression  for  the  partition  function  of  the  harmonic  oscillator  (Dauber  et  al 
1981),  it  is  seen  that  the  slope  for  the  frequency  dependence  of  entropy  is  —4.6  EU/logVj.  The 
extrapolation  of  the  harmonic  oscillator  treatment  to  low  frequency  motions  is  surely  not 
quantitative,  but  it  is  nonetheless  informative.  While  the  contribution  of  a  classical  vibration  is 
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the  order  of  5  cal/mol-deg  (EU),  giving  a  Gibbs  free  energy  of  1.5  kcal/mol  at  298K,  the 
contribution  of  the  5  MHz  relaxation,  seen  in  Figures  7  and  8,  would  be  some  30  cal/mol-deg  or 
9  kcal/mol,  and  that  for  the  acoustically  activated  mechanical  resonance  near  1  kHz,  reported  in 
Figure  6A,  would  approach  47  cal/mol-deg  or  14  kcal/mol.  From  this  it  should  be  apparent  that 
the  mechanical  relaxations  of  Figures  6,  7  and  8  arising  from  torsional  oscillations  are 
fundamental  in  providing  for  the  entropic  elasticity  of  elastin-based  systems.  On  a  more  general 
note,  computations  of  protein  structure  that  do  not  include  these  important  low  frequency 
contributions  to  the  free  energy  would  seem  to  be  limited  in  their  treatment  of  the 
structure/function  issue. 
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FIGURE  LEGENDS 


Figure  1:  Plot  of  the  percent  water  by  weight  as  a  function  of  temperature  for  poly(GVGVP) 
showing  a  stable  intermediate  structured  state  of  63%  water  by  weight  between  30  and  60°C,  a 
state  that  slowly  denatures  above  60°C  to  form  a  state  containing  32%  water  by  weight.  Plotted 
from  the  data  of  Urry  et  al  1985b.  Amphiphilic  (LCST,  lower  critical  solution  temperature) 
polymers  such  as  PNIPAM,  poly(N-isopropyl  acrylamide),  exhibit  only  a  single,  sharper 
transition,  and  it  is  to  what  is  generally  considered  a  disordered  state  with  about  30%  water  by 
weight.  Plotted  from  the  data  of  Grinberg  et  al  1999. 

Figure  2:  Schematic  diagram  of  an  atomic  force  microscope  apparatus  adapted  for  varying  the  z- 
direction  for  obtaining  force-extension  curves  of  single-chains. 

Figure  3:  Single-chain  force-extension  curves  for  Cys-(GVGVP)nx25i-Cys  at  a  temperature 
below  the  onset  temperature  for  hydrophobic  folding  and  assembly  for  this  composition.  The 
initial  trace  is  given  at  the  bottom  and  subsequent  traces  of  the  same  chain  without  intervening 
detachment  are  displaced  vertically  by  250  pN.  The  second  and  fifth  curves  exhibit  perfect 
reversibility.  The  curves  marked  by  an  asterisk  and  the  initial  curve  were  the  result  of  a  period 
of  at  least  30  s  in  the  relaxed  state,  where  one  possibility  is  that  some  hydrophobic  back-folding 
occurs.  Reproduced  with  permission  Urry  et  al,  2002. 

Figure  4:  Single-chain  force-extension  curves  for  Cys-(GVGIP)nx260'Cys  at  a  temperature  above 
the  onset  temperature  for  hydrophobic  folding  and  assembly  for  this  composition.  In  each  case 
the  energy  expended  in  deformation  is  greater  than  that  recovered  on  relaxation.  In  addition 
there  appears  to  be  a  complex  relaxation  curve.  Reproduced  with  permission  Urry  et  al,  2002. 
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Figure  5:  A.  Block  diagram  of  the  instrumentation  used  to  determine  the  absorption  per  unit 
volume  of  elastic  samples  in  the  acoustic  range.  B.  Example  of  the  use  of  two  samples  that  differ 
in  length  by  a  factor  of  two  as  required  for  the  analysis  of  Sheiba  (1996)  to  be  utilized. 

Figure  6:  A.  Plots  of  the  loss  factor  (absorption  per  unit  volume)  for  a  series  of  elastomers.  The 
bottom  curve  is  of  natural  rubber  and  the  slightly  higher  structure-less  curve  is  for  polyurethane, 
which  is  recognized  as  a  rubber  with  a  high  loss  factor  in  the  acoustic  absorption  range.  The 
upper  four  curves  are  of  20  Mrad  cross-linked  (GVGVP^o.  As  the  temperature  is  raised  above 
the  temperature  for  the  inverse  temperature  transition  of  hydrophobic  folding  and  assembly  for 
this  composition,  the  intensity  for  a  localized  relaxation  grows  dramatically.  To  the  best  of  our 
knowledge,  these  model  elastic  protein-based  polymers  exhibit  the  greatest  sound  absorption  for 
the  100  Hz  to  10,000  Hz  frequency  range.  This  demonstrates  the  formation  of  a  regular  structure 
for  (GVGIP)n  in  contrast  to  the  lack  of  a  regular  structure  for  the  classical  rubbers.  Reproduced 
with  permission  from  Urry  et  al,  2002. 

B.  The  imaginary  (absorptive)  part  of  the  dielectric  permittivity  of  poly(GVGIP).  The  same 
phenomenon  of  increased  absorption  on  passing  through  the  temperature  range  of  the  inverse 
temperature  transition  is  observed  in  the  100  Hz  to  10,000  Hz  frequency  range  as  exhibited  by 
the  same  composition  in  the  acoustic  absorption  range.  In  this  case  the  temperature  is  increased 
to  40°C  and  the  intensity  of  the  relaxation  continues  to  rise. 

Figure  7:  A.  The  real  part  of  the  dielectric  permittivity  for  poly(GVGIP). 

B.  The  imaginary  (absorptive)  part  of  the  dielectric  permittivity  of  poly(GVGIP).  The  same 
phenomenon  of  increased  absorption  on  passing  through  the  temperature  range  of  the  inverse 
temperature  transition  is  observed  in  the  1  MHz  to  100  MHz  range  as  exhibited  by  the  same 
composition  in  the  acoustic  absorption  range  in  Figure  6.  Reproduced  with  permission  from 
Buchet  et  al  1988.  Such  low  frequency  motions  provide  a  source  of  entropic  elasticity  and  low 
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energy  for  the  structural  state  of  the  elastic  protein-based  polymers.  Reproduced  with  permission 
from  Buchet  et  al,  1988. 

Figure  8:  A.  The  real  part  of  the  dielectric  permittivity  of  a-elastin  showing  an  increasing 
relaxation  near  5  MHz  with  increasing  temperature  from  below  to  above  the  temperature  range 
for  the  inverse  temperature  transition.  The  data  for  a-elastin  is  nearly  identical  to  that  found  for 
poly(GVGVP)  shown  in  the  insert  for  a  single  temperature  below  and  one  above  that  of  the 
inverse  temperature  transition.  This  demonstrates  that  the  same  relaxation  occurs  in  a  fragment 
of  the  natural  elastic  fiber  as  has  been  so  thoroughly  characterized  in  the  model  elastic  proteins, 
poly(GVGVP)  and  poly(GVGIP).  Reproduced  with  permission  from  Urry  et  al,  1985a. 

B.  The  imaginary  (absorptive)  part  of  the  dielectric  permittivity  of  natural  fibrous  elastin 
showing  an  increasing  relaxation  near  5  MHz  with  increasing  temperature  from  below  (23°C)  to 
above  (58°C)  the  temperature  range  for  the  inverse  temperature  transition.  The  same  relaxation 
is  seen  at  58 °C  with  a  somewhat  higher  intensity  and  a  slightly  higher  frequency  is  found  for  the 
20  Mrad  cross-linked  poly(GVGVP).  It  seems  quite  evident  that  dynamic  structured  states  occur 
in  the  natural  elastic  fiber  as  demonstrated  for  the  model  elastic  protein-based  polymer.  Adapted 
with  permission  from  Luan  et  al,  1988. 

Figure  9:  Thermoelasticity  curves  for  A.  fibrous  elastin  and  B.  20  Mrad  cross-linked 
poly(GVGVP).  See  text  for  discussion.  Whether  or  not  a  near  zero  slope,  i.e.,  an  fE/f  «  1,  is 
obtained  depends  on  the  rate  at  which  the  temperature  is  raised.  This  is  due  to  a  slow 
denaturation  at  temperatures  above  60°C.  Reproduced  with  permission  from  Urry  1988a,  1988b. 

Figure  10:  Description  of  the  proposed  molecular  structure  of  poly(GVGVP).  See  text  for 
discussion. 
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Figure  11:  Unrolled  perspective  of  one-half  turn  of  a  P-spiral  showing  the  suspended  segment 
between  P-tums.  Due  to  the  near  collinearity  of  the  bonds  on  each  side  of  a  peptide  moiety,  the 
peptide  moiety  can  undergo  rocking  motions.  In  particular,  the  indicated  suspended  segment 
there  are  two  peptide  moieties  with  the  X|/4-<t>5  pair  of  torsion  angles  and  the  X|/5-<t>i  pair  of  torsion 
angles.  These  pairs  of  torsion  angles  couple  in  each  case  to  result  in  rocking  motions  of  the 
respective  peptides,  referred  to  as  peptide  librations.  Such  motions  provide  much  entropy  to  the 
P-spiral  structure  and  to  all  protein  structures  that  allow  for  these  coupled  motions  and  it  is  a 
decrease  in  this  chain  dynamics  on  extension  that  is  a  dominant  source  of  entropic  elasticity 
exhibited  by  single  chains. 

Figure  12:  Stereo  view  of  the  torsional  oscillations  possible  within  a  2  kcal/mol-residue  cut-off 
energy  for  the  central  pentamer  of  one  turn  (three  pentamers)  of  the  P-spiral  structure  of 
poly(GVGVP)  in  the  relaxed  state  in  A.  and  in  the  130%  extended  state  in  B.  See  text  for 
discussion.  Reproduced  with  permission  from  Urry  and  Venkatachalam,  1983. 

Figure  13:  Lambda  plots  of  the  torsion  angles  attending  the  peptide  moieties  of  the  suspended 
segment  of  poly(GVGVP)  and  an  analogue.  The  slopes  of  negative  one  demonstrate  the 
correlated  motions  of  the  P-spiral  structure  of  poly(GVGVP)  in  A  and  B  and  of  the  D-Ala  (A’) 
analog,  poly(GVA’VP).  Extension  to  130%  dramatically  decreases  the  number  of  states  of  the 
polymer  that  are  available  within  a  given  energy  cut-off  as  does  the  substitution  of  the  D-residue. 
The  decrease  in  number  of  states,  as  discussed  in  the  text,  represent  a  decease  in  entropy  of  the 
structure.  Reproduced  with  permission  from  Urry  et  al,  1982,  Urry  and  Venkatachalam,  1983, 
Venkatachalam  and  Urry  1986. 

Figure  14:  Representation  of  entropy  as  a  volume  in  a  A<t>j-A\|/j  configuration  space.  See  text  for 
further  discussion. 
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Figure  15:  Plot  of  the  above  equation  for  entropy  as  a  function  of  log(oscillator  frequency)  using 
the  partition  function  of  an  harmonic  oscillator.  On  the  left-hand  ordinate  is  plotted  the  entropy 
in  cal/mol-deg  (EU,  entropy  units),  and  on  the  right-hand  ordinate  is  plotted  the  Gibbs  free 
energy  in  kcal/mol  at  298°K.  This  illustrates  the  increasing  contribution  of  low  frequency 
oscillations  to  the  entropy  of  a  polypeptide  chain. 


Table  1  Listing  of  the  <(>  and  i|>  torsion  angles  for  the  central  six  pentamers  of  (VPGVG)n 
and  the  root  mean  square  (RMS)  fluctuations  of  those  torsion  angles,  i.e.,  the  Ac)),  and  , 
in  columns  for  the  relaxed  and  130%  extended  states  arising  from  45  picosecond 
equilibration  time  and  80  picosecond  of  molecular  dynamics  simulation  for  each  state. 
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Figure  1:  Plot  of  the  percent  water  by  weight  as  a  function  of  temperature  for  poly(GVGVP) 
showing  a  stable  intermediate  structured  state  of  63%  water  by  weight  between  30  and  60°C,  a 
state  that  slowly  denatures  above  60°C  to  form  a  state  containing  32%  water  by  weight.  Plotted 
from  the  data  of  Urry  et  al  1985b.  Amphiphilic  (LCST,  lower  critical  solution  temperature) 
polymers  such  as  PNIPAM,  poly(N-isopropyl  acrylamide),  exhibit  only  a  single,  sharper 
transition,  and  it  is  to  what  is  generally  considered  a  disordered  state  with  about  30%  water  by 
weight.  Plotted  from  the  data  of  Grinberg  et  al  1999. 
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Figure  2:  Schematic  diagram  of  an  atomic  force  microscope  apparatus  adapted  for  varying  the  z- 
direction  for  obtaining  force-extension  curves  of  single-chains. 
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Successive  traces  from  bottom  to  top,  shifted  by  250  pN.  Before  the  traces  with  the  *  where  taken, 
the  chain  was  at  least  for  30  s  in  the  relaxed  state. 


Figure  3:  Single-chain  force-extension  curves  for  Cys-(GVGVP)nx25]-Cys  at  a  temperature 
below  the  onset  temperature  for  hydrophobic  folding  and  assembly  for  this  composition.  The 
initial  trace  is  given  at  the  bottom  and  subsequent  traces  of  the  same  chain  without  intervening 
detachment  are  displaced  vertically  by  250  pN.  The  second  and  fifth  curves  exhibit  perfect 
reversibility.  The  curves  marked  by  an  asterisk  and  the  initial  curve  were  the  result  of  a  period 
of  at  least  30  s  in  the  relaxed  state,  where  one  possibility  is  that  some  hydrophobic  back-folding 
occurs.  Reproduced  with  permission  Urry  et  al,  2002 
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Figure  4:  Single-chain  force-extension  curves  for  Cys-(GVGIP)nx26o_Cys  at  a  temperature  above 
the  onset  temperature  for  hydrophobic  folding  and  assembly  for  this  composition.  In  each  case 
the  energy  expended  in  deformation  is  greater  than  that  recovered  on  relaxation.  In  addition 
there  appears  to  be  a  complex  relaxation  curve.  Reproduced  with  permission  Urry  et  al,  2002. 
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Figure  5:  A.  Block  diagram  of  the  instrumentation  used  to  determine  the  absorption  per  unit 
volume  of  elastic  samples  in  the  acoustic  range.  B.  Example  of  the  use  of  two  samples  that  differ 
in  length  by  a  factor  of  two  as  required  for  the  analysis  of  Sheiba  (1996)  to  be  utilized. 
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Figure  7:  A.  The  real  part  of  the  dielectric  permittivity  for  poly(GVGIP). 

B.  The  imaginary  (absorptive)  part  of  the  dielectric  permittivity  of  poly(GVGEP).  The  same 
phenomenon  of  increased  absorption  on  passing  through  the  temperature  range  of  the  inverse 
temperature  transition  is  observed  in  the  1  MHz  to  100  MHz  range  as  exhibited  by  the  same  composition 
in  the  acoustic  absorption  range  in  Figure  6.  Reproduced  with  permission  from  Buchet  et  al  1988.  Such 
low  frequency  motions  provide  a  source  of  entropic  elasticity  and  low  energy  for  the  structural  state  of 
the  elastic  protein-based  polymers.  Reproduced  with  permission  from  Buchet  et  al,  1988. 
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Permittivity  (Imaginary  Part,  e") 


Frequency  (MHz) 


Frequency  (MHz) 

Figure  8:  See  the  figure  legend  list. 


5. 


Temperature,  °C 


Figure  9:  Thermoelasticity  curves  for  A.  fibrous  elastin  and  B.  20  Mrad  cross-linked 
poly(GVGVP).  See  text  for  discussion.  Whether  or  not  a  near  zero  slope,  i.e.,  an  fE/f  «  1,  is 
obtained  depends  on  the  rate  at  which  the  temperature  is  raised.  This  is  due  to  a  slow 
denaturation  at  temperatures  above  60°C.  Reproduced  with  permission  from  Urry  1988a,  1988b. 
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B.  (3-tum  E.  (3-spiral 


h-1.8  nm-H 


F.  twisted  filament 


H — 5  nm — H 


Figure  10:  Description  of  the  proposed  molecular  structure  of  poly(GVGVP).  See  text  for 
discussion. 
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Figure  11:  Unrolled  perspective  of  one-half  turn  of  a  p-spiral  showing  the  suspended 
segment  between  P-tums.  Due  to  the  near  collinearity  of  the  bonds  on  each  side  of  a  peptide 
moiety,  the  peptide  moiety  can  undergo  rocking  motions.  In  particular,  the  indicated 
suspended  segment  there  are  two  peptide  moieties  with  the  pair  of  torsion  angles  and 
the  t|/5-<j>i  pair  of  torsion  angles.  These  pairs  of  torsion  angles  couple  in  each  case  to  result 
in  rocking  motions  of  the  respective  peptides,  referred  to  as  peptide  librations.  Such 
motions  provide  much  entropy  to  the  3-spiral  structure  and  to  all  protein  structures  that 
allow  for  these  coupled  motions  and  it  is  a  decrease  in  this  chain  dynamics  on  extension  that 
is  a  dominant  source  of  entropic  elasticity  exhibited  by  single  chains. 
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Figure  12:  Stereo  view  of  the  torsional  oscillations  possible  within  a  2  kcal/mol-residue  cut-off 
energy  for  the  central  pentamer  of  one  turn  (three  pentamers)  of  the  (3-spiral  structure  of 
poly(GVGVP)  in  the  relaxed  state  in  A.  and  in  the  130%  extended  state  in  B.  See  text  for 
discussion.  Reproduced  with  permission  from  Urry  and  Venkatachalam,  1983. 
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Figure  13:  Lambda  plots  of  the  torsion  angles  attending  the  peptide  moieties  of  the  suspended  segment  of 
poly(GVGVP)  and  an  analogue.  The  slopes  of  negative  one  demonstrate  the  correlated  motions  of  the  P- 
spiral  structure  of  poly(GVGVP)  in  A  and  B  and  of  the  D-Ala  (A’)  analog,  poly(GVA’VP).  Extension  to 
130%  dramatically  decreases  the  number  of  states  of  the  polymer  that  are  available  within  a  given  energy 
cut-off  as  does  the  substitution  of  the  D-residue.  The  decrease  in  number  of  states,  as  discussed  in  the 
text,  represent  a  decease  in  entropy  of  the  structure.  Reproduced  with  permission  from  Urry  et  al,  1982, 
Urry  and  Venkatachalam,  1983,  Venkatachalam  and  Urry  1986. 
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Figure  14:  Representation  of  entropy  as  a  volume  in  a  configuration 

space.  See  text  for  further  discussion. 
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Figure  15:  Plot  of  the  above  equation  for  entropy  as  a  function  of  log(oscillator  frequency)  using  the 
partition  function  of  an  harmonic  oscillator.  On  the  left-hand  ordinate  is  plotted  the  entropy  in 
cal/mol-deg  (EU,  entropy  units),  and  on  the  right-hand  ordinate  is  plotted  the  Gibbs  free  energy  in 
kcal/mol  at  298°K.  This  illustrates  the  increasing  contribution  of  low  frequency  oscillations  to  the 
entropy  of  a  polypeptide  chain. 
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T  Si,  kcal/mole  pentamer  (298°K) 
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Abstract 

Description  of  the  mechanics  of  elastin  requires  the  understanding  of  two  interlinked  but  distinct  physical  processes: 
the  development  of  entropic  elastic  force  and  the  occurrence  of  hydrophobic  association.  Elementary  statistical- 
mechanical  analysis  of  AFM  single-chain  force-extension  data  of  elastin  model  molecules  identifies  damping  of 
internal  chain  dynamics  on  extension  as  a  fundamental  source  of  entropic  elastic  force  and  eliminates  the 
requirement  of  random  chain  networks.  For  elastin  and  its  models,  this  simple  analysis  is  substantiated 
experimentally  by  the  observation  of  mechanical  resonances  in  the  dielectric  relaxation  and  acoustic  absorption 
spectra,  and  theoretically  by  the  dependence  of  entropy  on  frequency  of  torsion-angle  oscillations,  and  by  classical 
molecular-mechanics  and  dynamics  calculations  of  relaxed  and  extended  states  of  the  p-spiral  description  of  the 
elastin  repeat,  (GVGVP)„.  The  role  of  hydrophobic  hydration  in  the  mechanics  of  elastin  becomes  apparent  under 
conditions  of  isometric  contraction.  During  force  development  at  constant  length,  increase  in  entropic  elastic  force 
resulting  from  decrease  in  elastomer  entropy  occurs  under  conditions  of  increase  in  solvent  entropy.  This  eliminates 
the  solvent  entropy  change  as  the  entropy  change  that  gives  rise  to  entropic  elastic  force  and  couples  association  of 
hydrophobic  domains  to  the  process.  Therefore,  association  of  hydrophobic  domains  within  the  elastomer  at  fixed 
length  stretches  interconnecting  dynamic  chain  segments  and  causes  an  increase  in  the  entropic  elastic  force  due  to 
the  resulting  damping  of  internal  chain  dynamics.  Fundamental  to  the  mechanics  of  elastin  is  the  inverse 
temperature  transition  of  hydrophobic  association  that  occurs  with  development  of  mechanical  resonances  within 
fibrous  elastin  and  polymers  of  repeat  elastin  sequences,  which,  with  design  of  truly  minimal  changes  in  sequence, 
demonstrate  energy  conversions  extant  in  biology  and  demonstrate  the  special  capacity  of  bound  phosphates  to 
raise  the  free  energy  of  hydrophobic  association. 


Introduction 

As  a  brief  introduction  to  the  mechanics  of  elastin, 
noted  are  the  composition  and  sequence  of  bovine 
elastin,  the  inverse  phase  transition  behavior  of  the 
elastin  and  related  protein  components,  and  the  confor¬ 
mational  aspects  of  the  most  representative  repeating 
sequence  of  elastin. 

Composition  and  sequence  of  elastin 

The  single-residue  7\-based  hydrophobicity  plots  of 
Figure  1  (Urry,  1997)  provide  a  simple  way  to  gain 
perspective  of  the  composition  and  sequence  of  bovine 
elastin  and  to  make  comparisons  with  other  extracellu¬ 
lar  proteins,  for  example  human  fibronectin.  The  longest 
repeating  sequence  (GVGVP)n  is  labeled  W4  based  on 
the  nomenclature  of  the  Sandberg  group  (Sandberg 
et  al.,  1981,  1985).  Additional,  less  extensive  repeats  are 
also  apparent.  Of  the  nearly  40  lysine  (Lys,  K)  residues, 
all  but  two  or  three  are  used  in  forming  cross-links. 
Generally,  four  Lys  residues  come  together  to  form  a 
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tetrasubstituted  pyridinium.  By  contrast,  in  the  vulca¬ 
nization  of  rubber  every  monomer  is  a  potential  cross¬ 
link,  which  allows  for  cross-link  formation  whenever 
two  chains  are  in  contact.  In  elastin  only  5%  of  the 
residues  participate  in  cross-link  formation  and  these  are 
often  distributed  along  the  chain  in  pairs.  Such  a 
situation  requires  a  regular  structure  in  order  to  bring 
this  limited  number  of  residues  into  adequate  proximity 
for  covalent  cross-link  formation. 

There  are  only  three  negatively  charged  aspartic  acid 
(Asp,  D)  and  no  glutamic  acid  residues  in  bovine  elastin. 
The  more  hydrophobic  phenylalanine  (Phe,  F)  and 
tyrosine  (Tyr,  Y)  residues  are  also  apparent.  The  almost 
complete  absence  of  charged  groups  combined  with  a 
substantial  quantity  of  hydrophobic  residues  are  re¬ 
sponsible  for  the  hydrophobic  association  between 
chains  that  result  in  parallel-aligned  twisted  filaments 
as  seen  in  electron  micrographs  of  negatively  stained 
tropoelastin  and  a-elastin  (Cox  et  al. ,  1973,  1974)  and  in 
fibrous  elastin  (Gotte  et  al. ,  1974). 

By  comparison,  a  sequence  that  results  in  a  series  of 
globular  repeating  units  of  hydrophobically  folded 
p-barrels  is  seen  in  the  single-residue  Tt-based  hydro¬ 
phobicity  plots  for  fibronectin  in  the  lower  part  of 
Figure  1  (Komblihtt  et  al .,  1985).  In  this  case  there  are 
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Fig.  1.  Single  residue  7>based  hydrophobicity  plots  for  bovine  elastin  (A)  and  human  fibronectin  (B)  (see  text  for  discussion).  Part  A  reproduced 
with  permission  from  Urry  et  at.  (1995),  Part  B  reproduced  with  permission  from  Urry  and  Luan  (1995a). 


many  more  charged  Glu,  Asp  and  Lys  residues,  and  the 
periodicity  of  the  repeating  globular  units  becomes 
apparent  through  the  repeating  tryptophan  (Trp,  W) 
residues.  This  most  hydrophobic  residue,  W,  repeating 
approximately  every  90  residues  identifies  type  III 
domains  of  fibronectin,  which  also  appear  as  a  repeating 
unit  in  titin  (connectin).  Such  repeating  globular  units 
give  a  sawtooth  pattern  to  the  single-chain  force- 
extension  curves  of  titin  (Rief  et  al.,  1997;  Gaub  and 
Fernandez,  1998),  which  is  very  different  from  the 
smooth  single-chain  force-extension  curves  given  by 
elastin  models  (Urry  et  al. ,  2002a,  b). 

Inverse  temperature  transition  of  elastin  and  component 
sequences 

Increase  in  order  of  elastic  protein  upon  raising 
the  temperature 

Tropoelastin  (precursor  protein  of  fibrous  elastin),  a- 
elastin  (chemical  fragmentation  product  from  fibrous 


elastin)  and  high-molecular-weight  polymers  of  repeating 
sequences  of  elastin  are  water  soluble  at  low  temperatures, 
but  produce  phase  separation  upon  raising  the  tempera¬ 
ture.  The  initial  aggregates  of  the  phase  transition,  when 
placed  on  a  carbon-coated  EM  grid  and  negatively 
stained  with  uranyl  acetate/oxalic  acid,  are  observed  to 
form  parallel  aligned  filaments  comprised  of  twisted 
filaments  with  interesting  periodicities  in  the  optical 
diffraction  patterns  of  the  micrographs  (Volpin  et  al, 
1976a,  b).  Furthermore,  cyclic  analogues  containing 
pentapeptide  and  hexapeptide  repeats  crystallize  during 
temperature  increase  and  redissolve  upon  lowering  the 
temperature.  Clearly,  the  order  in  elastin  and  elastin- 
based  polymers  increases  with  increasing  the  temperature. 

This  ordering  of  protein  upon  raising  the  temperature 
is  consistent  with  the  second  law  of  thermodynamics:  a 
structured  hydration  that  surrounds  the  dissolved 
hydrophobic  groups  becomes  less  ordered  bulk  water 
during  the  phase  separation  of  hydrophobic  association 
(Urry  et  al 1997a,  b). 
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This  phase  transition  leading  to  a  fundamental 
increase  in  order  of  the  polymer  on  raising  the  temper¬ 
ature  is  called  an  inverse  temperature  transition. 

Phase  diagram 

Phase  diagrams  of  the  inverse  temperature-dependent 
behavior  of  elastin-based  polymers  is  seen  in  Figure  2 
(Sciortino  et  al ,  1990, 1993).  In  the  usual  polymer  phase 
diagram  the  polymers  are  insoluble  below  and  soluble 
above  the  binodal  (coexistence  or  coacervate)  line.  For 
elastin-based  polymers  solubility  is  inverted.  With  ela¬ 
stin-based  polymers  solubility  occurs  below  and  insolu¬ 
bility  occurs  above  the  coexistence  line.  The  shape  of  the 
coexistence  line  is  also  inverted:  while  for  the  usual 
petroleum-based  polymers  the  line  is  concave  towards 
the  concentration  (volume-fraction)  axis,  for  elastin- 
related  polymers  it  is  convex.  In  case  of  the  elastin-based 
polymers  the  coexistence  line  is  also  called  the  Tt-divide. 
Each  point  (called  Tt)  along  the  7>divide  corresponds, 
for  a  given  concentration,  to  the  onset  temperature  of 
aggregation. 

Elastin-based  polymers  exhibit  a  latent  heat  of  tran¬ 
sition.  Accordingly,  within  the  temperature  range  of  the 
transition  the  chemical  potential  of  the  polymers  in 
solution  and  in  the  phase-separated  state  are  the  same. 
This  enables  the  derivation  of  the  Gibbs  free  energy  of 
hydrophobic  association,  AGHa  (see  below).  The  more 
hydrophobic  elastin-based  polymers  that  associate  at 
lower  temperatures  have  lower  values  of  AGHa- 

Conformational  aspects  of  the  (GVGVP)„  of  elastin 

The  structure  of  (GVGVP)„,  shown  in  Figure  3,  was 
developed  by  using  methodologies  of  peptide  synthesis, 
NMR,  molecular-mechanics  and  dynamics  calculations 


volume  fraction 


Fig.  2.  Schematic  representation  of  phase  diagrams  for  several  model 
elastic  proteins  based  on  the  elastin  repeat,  (GVGVP),,  (see  text  for 
discussion).  Solid  curves  adapted  with  permission  from  Sciortino  et  al. 
(1990)  and  (1993). 


(Venkatachalam  and  Urry,  1981;  Chang  et  al,  1989), 
crystal  structure  of  cyclic  analogues  (Cook  et  al.,  1980) 
and  the  NMR-  and  computationally  characterized 
relationship  between  cyclic  and  linear  structures  (Urry 
et  al,  1981,  1989),  and  Raman,  absorption  and  circular 
dichroism  spectroscopies. 

The  structure  involves  a  series  of  p-turns,  ten-atom 
hydrogen-bonded  rings  from  the  Val1  C— O  to  the  Val4 
N— R  (Figure  3A  and  B)  which,  upon  raising  the  tem¬ 
perature,  wrap  up  into  a  helical  structure  called  the  p- 
spiral  (shown  schematically  in  Figure  3C  and  D  and  in 
detail  in  Figure  3E).  It  is  shown  that  the  P-turns 


B.  P-tum  E.  (3-spiral 


H.8  nnrH 


F.  twisted  Filament 
of  (3-spirals 


H — 5  nm — H 


Fig.  3.  Representations  of  the  proposed  molecular  structure  of  (GVGVP)„  (see  text  for  discussion). 
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function  as  hydrophobic  spacers  between  the  turns  of 
the  0-spiral  and  that  the  0-spirals  associate  hydropho- 
bically  to  form  twisted  filaments  of  the  dimensions 
found  in  the  transmission  electron  micrographs  of 
negatively  stained  incipient  aggregates  as  noted  above 
(Volpin  et  al .,  1976a,  b).  Hydrophobic  folding  of  0- 
spirals  and  the  hydrophobic  association  of  0-spirals  to 
form  twisted  filaments  occur  in  a  cooperative  manner. 


General  considerations  of  entropic  elastic  force 

Components  of  elastic  force  and  their  delineation 

Internal  energy  (fe)  and  entropy  (f$)  components 
of  elastic  force 

The  total  elastic  force  (/)  can  be  thermodynamically  de¬ 
scribed  as 

f  =  (dE/dL)ViTin~T(dS/dL)VT^  (1) 

where  E  is  the  internal  energy;  S  the  entropy;  T  the 
absolute  temperature  in  °K,  and  V,  T ,  and  n  indicate 
that  the  change  in  length  (6L)  occurs  at  constant 
volume,  temperature  and  composition.  Accordingly, 
the  total  elastic  force  comprises  two  components,  the 
internal  energy  component  (/*E)  and  the  entropic  com¬ 
ponent  (fs): 

/  =  /e+/s-  (2) 

Experimental  delineation  of  the  relative  magnitudes 
of  fa  andfs 

Following  Flory  et  al.  (1960)  the  total  force  can  also  be 
written  as 

f  =  (dE/dL)V'Ttn  +  m/  QT)yU  (3) 

which  is  equivalent  to 

A/f--n*W/r\/*T)w  (4) 

Equation  (4)  shows  that  the  ratio  of  the  internal  energy 
component  of  force  to  the  total  force  can  be  determined 
experimentally  by  plotting  ln[//7]  as  a  function  of 
temperature  while  maintaining  the  elastic  element  at 
constant  volume,  at  fixed  length  and  without  a  change  in 
composition.  Under  these  experimental  conditions,  the 
slope  of  the  plot  multiplied  by  (-7°K)  provides  th e/E// 
ratio. 

Statistical  mechanical  expression  for  entropy 
The  Boltzmann  relation 

The  Boltzmann  relation  provides  the  bridge  from  a 
statistical  mechanical  description  of  molecular  structure 
to  experimentally  determined  thermodynamic  quanti¬ 
ties.  It  is  an  elegant  yet  simple  statement  of  entropy  OS) 


in  terms  of  thermodynamic  probability  (W,  the  number 
of  a  priori  equally  probable  states  accessible  to  the 
system)  (Eyring  et  al .,  1964),  i.e., 

S  =  R\nW .  (5) 

R  (1.987  cal/deg  mol)  is  the  gas  constant.  R=Nk ,  where 
N  is  Avogadro’s  number  (6.02  x  1023/mol)  and  k  Boltz¬ 
mann’s  constant  (1.38  x  10-16  erg/deg  K).  IF  is  a  vol¬ 
ume  in  phase  space,  a  2AT-dimensional  space,  which 
fundamentally  provides  a  description  of  a  molecule  in 
terms  of  the  momentum  (pt)  and  coordinate  (q,)  of  each 
its  i  atoms. 

Description  of  entropy  using  partition  functions 
for  different  degrees  of  freedom 

In  practice,  the  description  of  the  molecule  is  achieved 
by  the  product  of  partition  functions,  which  group 
according  to  the  degrees  of  freedom  of  the  molecular 
system.  In  general,  the  IN  degrees  of  freedom  group  as 
three  translational  degrees  of  freedom,  three  rotational 
degrees  of  freedom  and  the  remainder  are  internal 
degrees  of  freedom  that  comprise  vibrations  and  tor¬ 
sional  oscillations  (rotations  about  bonds).  In  the 
measurement  of  elasticity,  the  single  chain  molecule  or 
the  cross-linked  matrix  is  fixed  at  both  ends.  In  this  case, 
there  are  holonomic  constraints  on  the  molecular  system 
in  that  there  are  neither  whole-molecule  translational 
nor  rotational  degrees  of  freedom.  Thus,  we  are  left  only 
with  internal  chain  dynamics. 

Frequency  dependence  of  entropy  for  the  harmonic 
oscillator  representation  of  internal  chain  dynamics 
The  harmonic  oscillator  partition  function  (fv)  provides 
a  relatively  simple  and  yet  informative  representation  of 
internal  chain  dynamics. 

fv  =  [1  -  exp (-hvi/kT)]~l.  (6) 

When  placed  in  Equation  (5)  and  expressed  in  terms  of 
the  frequency-dependence  of  entropy  (S)  we  obtain 
(Dauber  et  ai,  1981), 

Si  =^{ln[l  -  cxp(-hVi/kT)]~l 

+  (hVi/kT)[exV(hvi/kT)  -  I]" 1 }.  (7) 

Equation  (7)  is  plotted  in  Figure  4  with  S{  on  the  left- 
hand  ordinate  and  with  the  free-energy  contribution 
(TSt)  on  the  right-hand  ordinate. 

The  interesting  features  of  Figure  4  are  that  the 
vibrational  modes  occur  at  high  frequencies  where  the 
contribution  to  entropy  is  small,  and  the  torsional 
oscillations  occur  at  much  lower  frequencies  where  the 
contribution  to  entropy  is  much  larger.  As  discussed 
below,  the  mechanical  resonances  exhibited  by  elastin 
models  and  fibrous  elastin  occur  at  frequencies  near 
1  kHz  and  5  MHz,  i.e.,  near  3  and  6  on  the  abscissa  of 
Figure  4.  Accordingly,  it  becomes  reasonable  to  neglect 


S,  =  tf[1n(l  -  e^,tT )“*  +  ( hv ,  tkT^e^"  - 1)"1  ] 
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Fig.  4.  Plot,  based  on  the  harmonic  oscillator  partition  function,  of  the 
entropy  of  the  oscillator  as  a  function  of  log(oscillator  frequency).  On 
the  left-hand  ordinate  is  plotted  the  entropy  in  cal/mol  deg  (EU, 
entropy  units),  and  on  the  right-hand  ordinate  is  plotted  the  entropic 
contribution  to  the  Gibbs  free  energy  in  kcal/mol  pentamer  at  298°K. 
This  illustrates  the  increasing  contribution  of  low  frequency  oscilla¬ 
tions  to  the  entropy  and  free  energy  of  a  protein  chain  segment. 
Adapted  with  permission  from  Urry  et  al.  (1988). 


changes  in  vibrational  frequencies  on  extension  and  to 
consider  those  configurational  and  dynamic  changes  due 
to  changes  in  torsion  angles. 

Expression  for  the  change  in  entropy  on  extension 
The  change  in  entropy  upon  extension  is  written, 

AS  ={SC-  ST)  =  R\n(Wc/WT),  (8) 

where  e  and  r  stand  for  the  extended  and  relaxed  states. 
With  Equation  (8)  it  becomes  possible  to  use  both 
molecular  mechanics  and  dynamics  calculations  of  the 
molecular  structure  described  in  Figure  3  to  calculate 
the  contribution  of  torsional  freedom  of  the  structure  to 
entropy  in  both  the  relaxed  and  extended  states. 


The  physical  basis  for  entropic  elasticity  in  elastin 

In  an  ideal  or  perfect  elastomer  the  energy  repeatedly 
invested  in  extension  is  repeatedly  and  completely 
recovered  during  relaxation.  Ideality  increases  as  the 
elastic  force  results  from  a  decrease  in  entropy  upon 
extension,  because  this  occurs  without  stressing  bonds  to 
the  breaking  point.  Elastin  models  and  elastin  itself  in 
water  provide  examples  of  such  entropic  elastomers  with 
about  90%  of  the  elastic  force  being  entropic,  that  is,  the 
/E//ratio  of  Equation  (4)  is  about  0.1.  This  is  essential  to 
human  life  expectancy,  because  the  half-life  of  elastin  in 
the  mammalian  elastic  fiber  is  on  the  order  of  70  years. 
This  means  that  the  elastic  fibers  of  the  aortic  arch  and 
thoracic  aorta,  where  there  is  twice  as  much  elastin  as 
collagen,  will  have  survived  some  billion  demanding 
stretch-relaxation  cycles  by  the  start  of  the  seventh 
decade  of  life.  This  represents  an  ultimate  in  ideal 
elasticity. 
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Three  main  mechanisms  proposed  for  the  entropic 
elasticity  of  elastin 

The  classical  (random  chain  network)  theory  of  rubber 
elasticity 

In  1958  Hoeve  and  Flory  reported  studies  on  the  nature 
of  elasticity  of  the  bovine  fibrous  protein,  elastin,  as 
representative  of  the  mammalian  elastic  fiber.  Because 
they  were  able  to  determine  a  very  low  /E//  ratio,  they 
concluded  that  the  elastic  fiber  was  without  order,  that 
is,  it  was  a  random  chain  network.  They  re-affirmed  it  in 
the  spring  of  1974  in  a  Biopolymers  paper  (Hoeve  and 
Flory,  1974)  stating  that  ‘A  network  of  random  chains 
within  the  elastic  fibers,  like  that  in  a  typical  rubber,  is 
clearly  indicated.’  In  order  to  emphasize  this  perspective 
further,  the  following  statement  appeared  in  the  legend 
to  Figure  1  of  that  paper:  ‘Configurations  of  chains 
between  cross-linkages  are  much  more  tortuous  and 
irregular  than  shown.’  In  part  because  Paul  Flory 
received  the  Nobel  Prize  in  1974  for  his  work  on 
macromolecules,  this  perspective  became  firmly  set  in 
the  minds  of  the  interested  scientific  community. 

The  solvent  (bulk  water  <->  hydrophobic  hydration) 
entropy  mechanism 

It  appeared  that  the  sole  purpose  of  the  Hoeve  and 
Flory  1974  paper,  which  presented  neither  new  experi¬ 
mental  data  nor  theoretical  analysis,  was  to  refute  the 
publication  of  a  new  mechanism.  The  new  mechanism 
for  the  entropic  elasticity  of  elastin  was  published  by 
Weis-Fogh  and  Anderson  (1970)  with  the  title  ‘New 
Molecular  Model  for  the  Long-range  Elasticity  of 
Elastin.’  The  proposed  new  mechanism,  stated  in  the 
terms  of  present  adherents,  was  that  upon  extension 
hydrophobic  side-chains  become  exposed  to  the  water 
solvent.  The  result  is  formation  of  low-entropy  water  of 
hydrophobic  hydration.  Here  we  refer  to  this  as  the 
solvent  (bulk  water  <-»  hydrophobic  hydration)  mecha¬ 
nism  for  entropic  elasticity.  This  mechanism  has  been 
given  credibility  by  groups  using  computational  meth¬ 
ods  on  (GVGVP)„  in  water  (Wasserman  and  Salemme, 
1990;  Alonso  et  aly  2001).  One  of  the  many  arguments 
marshaled  by  Hoeve  and  Flory  (1974)  against  this 
mechanism  was  that  polymer  backbone,  not  solvent, 
must  bear  the  force. 

The  damping  of  internal  chain  dynamics  on  extension 
In  the  hydrophobically  folded  (3-spiral  structure  (Fig¬ 
ure  3),  the  presence  of  suspended  segments  between  the 
(3-tums  is  immediately  apparent.  With  water  surround¬ 
ing  the  suspended  segments  and  no  steric  hindrance  to 
limit  rotation  about  backbone  bonds  of  the  suspended 
segments,  the  peptide  moieties  would  be  expected  to 
rock  or  ‘libra te’.  These  suspended  segments  are  prob¬ 
ably  free  to  undergo  large  torsion-angle  oscillations, 
and  the  amplitude  of  these  torsional  oscillations  might 
be  damped  during  extension  (Urry,  1982).  We  devised 
experimental  (physical  and  chemical)  and  computation¬ 
al  tests  of  the  structure  with  a  focus  on  the  freedom  of 
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motion  of  the  suspended  segments.  The  chemical  tests 
(replacement  of  the  Gly  residues  with  D-  and  L- 
alanines)  resulted  in  the  predicted  limited  torsional 
oscillations  and  elasticity  (Urry  et  al. ,  1983a,  b,  1984, 
1991).  The  principal  physical  tests  involved  dielectric 
and  NMR  relaxation  methodologies  (Henze  and  Urry, 
1985;  Urry  et  al .,  1985a,  b,  c,  1986;  Buchet  et  al.,  1988). 
Remarkably,  in  the  dielectric  relaxation  characteriza¬ 
tion,  where  the  only  dipole  moments  were  those  of  the 
backbone  peptide  moieties,  intense  localized  relaxations 
developed  as  the  model  elastic  protein  underwent  the 
inverse  temperature  (phase)  transition.  Furthermore, 
molecular  mechanics  and  dynamics  calculations  dem¬ 
onstrated  the  decrease  in  amplitude  of  torsion-angle 
oscillations  within  the  suspended  segment  during  exten¬ 
sion  (Urry,  1982;  Chang  and  Urry,  1989).  Indeed,  by 
treating  the  dynamics  of  the  pentameric  unit  the 
oscillating  dipole  moment  reproduced  the  dielectric 
relaxation  results  (Venkatachalam  and  Urry,  1986). 
Thus,  the  concept  of  the  damping  of  internal  chain 
dynamics  by  extension  as  a  source  of  entropic  elastic 
force  became  established  (Chang  and  Urry,  1989).  Key 
experimental  and  computational  results  are  treated  in 
more  detail  below. 

Key  experimental  data  relevant  to  the  proposed 
mechanisms 

Atomic  force  microscopy  (AFM)  single-chain 
force-extension  results 

The  remarkable  capacity  to  obtain  single-chain  force- 
extension  curves  provides  new  insights  into  the  mecha¬ 
nism  of  entropic  elasticity.  Figure  5A  shows  a  series  of 
single-chain  force-extension  curves  for  a  single  molecule 
of  Cys-(GVGVP)„x25i-Cys  (Urry  et  al. ,  2002a).  Starting 
from  the  bottom,  a  series  of  extension-relaxation  cycles 
are  displayed.  The  second  trace,  during  which  the 
molecule  was  extended  from  an  end-to-end  distance  of 
200-700  nm  and  then  returned  to  the  starting  length, 
demonstrates  perfect  reversibility  at  the  rate  of  stretch 
and  within  the  resolution  of  the  technique.  Since  all  of 
the  energy  of  deformation  was  recovered  during  relax¬ 
ation,  we  observe  an  ideal  elastomer!  For  the  third, 
fourth  and  sixth  traces,  a  resting  time  of  30  s  near 
200  nm  was  allowed,  and  these  traces  demonstrated 
barely  detectable  imperfect  overlap  of  extension  and 
relaxation  curves.  The  fifth  trace  without  the  30  s  resting 
time  near  200  nm  again  exhibited  perfect  reversibility. 
Finally  on  the  seventh  extension  of  the  chain  the 
continuity  from  substrate  to  cantilever  tip  became 
severed  and  the  force  dropped  to  zero  just  above  700  nm. 

The  slight  hysteresis  of  traces  one,  three,  four  and  six, 
may  be  due  to  backfolding  of  the  chain  on  itself, 
interaction  of  a  second  chain  picked  up  at  low  extension 
or  possibly  even  non-specific  adsorption  of  the  single 
chain  to  surface  at  low  extension  and  at  a  position  along 
the  chain  of  within  100  nm  of  extended  length  from  the 
attachment  point.  Regardless  of  the  source  of  the  barely 
detectable  non-overlap  of  several  of  the  traces,  ideal 


Fig.  5 .  AFM  single-chain  force-extension  curves  of  the  model  elastic 
proteins  based  on  elastin,  Cys-(GVGVP)„X2si-Cys  (A)  and  Cys- 
(GVGIP)„x26o-Cys  (B),  showing  the  ideal  elasticity  due  to  perfect 
reversibility  of  traces  2  and  5  from  bottom  of  A  and  marked  hysteresis 
in  the  lower  trace  of  B.  See  text  for  discussion.  Reproduced  with 
permission  from  Urry  et  al.  (2002a). 


(entropic)  elasticity  has  been  observed  in  a  single  chain. 
No  longer  can  the  change  from  a  Gaussian  distribution 
of  end-to-end  chain  lengths  within  a  random  chain 
network  be  insisted  upon  as  the  structural  representa¬ 
tion  of  entropic  elasticity.  The  ideal  single-chain  force- 
extension  curves  require  high  dilution,  and  a  higher 
dilution  yet  is  required  for  more  hydrophobic  protein- 
based  polymers  where  the  tendency  for  aggregation  is 
greater.  This  and  additional  data  tell  us  that  these 
perfectly  reversible  force-extension  curves  are  due  to 
single-chains.  A  single  chain  does  not  constitute  a 
random  chain  network;  and  a  single  chain  with  ends 
fixed  in  space  does  not  comprise  a  Gaussian  distribution 
of  end-to-end  chain  lengths.  It  seems  that  this  eliminates 
the  basic  tenets  of  the  random  chain  network  theory  of 
entropic  elasticity  as  regards  these  elastic  protein-based 
polymers.  Furthermore,  random  chains  do  not  exhibit 
mechanical  resonances,  e.g.,  intense  localized  Debye- 
type  relaxations  as  observed  in  the  dielectric  relaxation 
spectra  near  5  MHz  and  3  kHz,  and  intense  localized 
absorption  maxima  as  observed  in  the  acoustic  absorp¬ 
tion  and  loss  permittivity  spectra.  These  points  consti¬ 
tute  only  the  more  apparent  elements  of  the  argument. 

In  addition,  on  the  basis  of  the  elementary  statistical 
mechanical  analysis  given  above,  it  would  seem  that  the 
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source  of  entropic  elasticity  must  come  from  a  decrease 
in  internal  chain  dynamics  upon  extension. 

Reduction  of  mean  solvent-entropy  change  increases 
rather  than  decreases  entropic  elastic  force 
Formation  of  hydrophobic  hydration  is  exothermic 
(Butler,  1937).  Accordingly,  when  the  temperature  of 
the  dissolved  protein  with  its  hydrophobic  hydration  is 
raised  from  below  to  above  the  inverse  temperature 
transition,  an  endothermic  transition  due  to  the  conver¬ 
sion  of  hydrophobic  hydration  to  bulk  water  occurs. 
The  transition  from  hydrophobic  hydration  to  bulk  water 
represents  a  positive  change  in  entropy.  By  finding  a 
suitable  solvent  that  allows  the  transition  but  reduces 
the  heat  of  the  transition  to  near  zero,  it  becomes 
possible  to  determine  whether  a  decrease  in  elastic  force 
occurs  due  to  the  loss  of  the  contribution  of  the  negative 
solvent-entropy  change  to  the  entropic  elastic  force.  In 
other  words,  if  a  decrease  in  solvent  entropy  occurs  as 
bulk  water  becomes  ordered  (in  the  form  of  hydropho¬ 
bic  hydration)  by  exposure  of  hydrophobic  groups 
during  extension,  then  the  extent  to  which  this  contri¬ 
butes  to  entropic  elastic  force  could  be  measured  as  a 
decrease  in  entropic  elastic  force. 

In  our  experimental  approach  (based  on  Hoeve  and 
Flory  (1958))  ethylene  glycol  was  added  to  suppress  the 
solvent  interaction  with  the  protein.  In  particular,  as 
shown  in  the  differential  calorimetry  data  of  Figure  6A 
increasing  the  amount  of  ethylene  glycol  decreases  the 
heat  and  lowers  the  temperature  of  the  transition.  By 
30%  ethylene  glycol  in  water  the  heat  of  the  transition  is 
minimal.  The  entropy  change  attending  the  transition  is 
calculated  by  dividing  an  increment  of  heat  released  over 
the  temperature  of  the  increment  and  summing  over  all 
increments  of  the  transition.  Accordingly,  the  decrease 
in  heat  of  the  transition  seen  in  Figure  6A  would  be 
expected  to  result  in  a  decrease  in  the  entropic  elastic 
force  due  to  a  decrease  in  the  magnitude  of  the  negative 
entropy  of  hydration  on  extension.  In  fact,  as  seen  in 
Figure  6B,  the  entropic  elastic  force  reached  is  greater  in 
30%  ethylene-70%  water. 

The  thermoelasticity  studies  on  y-irradiation-cross- 
linked  poly(GVGVP)  of  Figure  6B  are  carried  out  by 
stretching  the  elastomer  to  a  fixed  extension  of  50%  at 
37°C.  Then  held  at  that  fixed  length,  the  temperature  is 
re-equilibrated  below  that  of  the  onset  for  the  inverse 
temperature  transition,  and  the  temperature  is  very 
slowly  raised  to  55°C.  As  regards  the  contribution  of 
solvent  entropy  change  to  entropic  elastic  force,  it  is 
expected  that  the  force  reached  would  be  less  for  the 
ethylene  glycol  containing  solution  by  a  fraction  indi¬ 
cative  of  the  contribution  of  solvent  entropy  change  to 
the  elastic  force.  The  observation  that  the  force  actually 
increased,  rather  than  decreased  does  not  favor  the 
solvent  entropy  mechanism. 

As  shown  in  Equation  (4),  the  value  of /e//is  given  by 
-T(d\n[f/T\/dT)KL%n.  Using  the  temperature  range  from 
45  to  55°C,  the  /E// ratio  is  found  to  be  no  more  than 
0.1.  Based  on  this  analysis  using  Equation  (4),  the/s// 
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Fig.  6.  (A)  Differential  calorimetry  curves  as  a  function  of  ethylene 
glycol  (EG)  in  water.  Note  the  decrease  in  the  heat  and  lowering  of  the 
temperature  range  of  the  transition  as  the  amount  of  ethylene  glycol  is 
raised  to  30%.  B.  Thermoelasticity  curves,  plots  of  log(//7)  vs.  T  at 
fixed  length,  for  20  Mrad  y-irradiation  cross-linked  poly(GVGVP). 
The /E//  ratio,  determined  from  the  slope  above  45°C  for  both  solvent 
conditions,  is  less  than  or  equal  to  0.1,  that  is  the  entropic  component 
of  elastic  force  is  90%  or  greater.  Significantly,  addition  of  ethylene 
glycol  (EG)  results  in  an  increase  in  entropic  elastic  force  suggesting 
that  solvent  entropy  change  is  not  a  contribution  to  the  entropic  elastic 
force.  Reproduced  with  permission  from  Luan  et  at.  (1989). 
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value  would  be  0.9  or  greater,  and  yet  this  entropic 
elastic  force  would  appear  not  to  come  from  a  change  in 
solvent  entropy.  Even  more  telling  is  the  increase  in 
force  at  fixed  length,  (df/dT)L,  as  the  temperature  is 
raised  through  the  range  of  the  inverse  temperature 
transition.  As  will  be  discussed  below,  this  thermally 
driven  isometric  contraction  shows  an  increase  in  force 
under  conditions  of  an  increase  rather  than  the  required 
decrease  in  solvent  entropy. 


Presence  of  mechanical  resonances  in  elastin  models 
and  elastin  itself 

(i)  Mechanical  resonances  in  the  acoustic  frequency 
range :  The  elastic  models  of  elastin  based  on  the 
pentameric  repeat  (GVGVP)  and  analogues  thereof, 
such  as  (GVGIP),  exhibit  mechanical  resonances  in  the 
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acoustic  frequency  range  of  100  Hz  to  7  kHz.  This 
is  seen  in  Figure  7A  for  y-irradiation-cross-linked 
(GVGIP)26o>  which  exhibits  a  mechanical  resonance  of 
increasing  absorption  intensity  as  the  temperature  is 
raised  from  below  to  above  the  range  of  the  inverse 
temperature  transition  (Urry  et  al .,  2002a).  By  contrast, 
natural  rubber  and  polyurethane  (a  particularly  good 
elastomer  for  sound  absorption),  exhibit  only  broad 
low-intensity  absorption  as  might  be  expected  for  a 
random  chain  network.  For  a  random  chain  the 
frequency  for  rotation  about  each  polymer  backbone 
bond  will  be  different,  whereas  for  an  elastomer  with  a 


regular  repeating  structure  such  as  that  in  Figure  3 
mechanical  resonances  such  as  those  in  Figure  7A  can 
result.  For  our  purposes  the  observation  of  a  mechan¬ 
ical  resonance  is  evidence  for  a  non-random  structure, 
and  therefore  challenges  the  classical  theory  of  rubber 
elasticity. 

Low-frequency  dielectric  relaxation  studies  also  dem¬ 
onstrate  mechanical  resonances  with  absorption  in  the 
acoustic  frequency  range.  As  seen  in  Figure  7B,  the 
imaginary  part  (e")  of  the  dielectric  permittivity  of 
cross-linked  (GVGIP)32o  exhibits  a  relaxation  that  also 
increases  in  intensity  as  the  temperature  is  raised  from 
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Fig .  7.  Mechanical  resonances  seen  in  the  acoustic  absorption  frequency  range.  The  repeating  pentamers,  as  they  fold  into  the  regularly  repeating 
structure  of  Figure  3  on  raising  the  temperature,  develop  a  mechanical  resonance  wherein  all  pentamers  absorb  energy  over  the  same  frequency 
range.  Although  the  physical  means  of  exciting  the  mechanical  resonance  is  different  in  A  and  B,  an  acoustic  wave  in  A  and  an  oscillating  electric 
field  in  B,  the  maxima  of  these  low  frequency  mechanical  resonances  are  only  shifted  by  a  few  kHz.  A  second  higher  frequency  mechanical 
resonance  is  observed  for  this  same  elastic  model  protein  in  the  dielectric  relaxation  near  5  MHz,  1000  kHz  to  higher  frequency,  as  seen  in 
Figure  8  for  both  the  real  and  imaginary  parts  of  the  dielectric  permittivity.  Part  A,  acoustic  absorption/unit  volume  (loss  factor),  reproduced 
with  permission  from  Urry  et  al.  (2002a)  and  part  B,  imaginary  part  of  the  low  frequency  dielectric  relaxation  spectra,  reproduced  with 
permission  from  Urry  et  al.  (2002b). 


551 


below  to  above  the  temperature  of  the  inverse  temper¬ 
ature  transition  (Urry  et  al. ,  2002b).  As  the  excitation 
energy  is  not  that  of  a  compressional  acoustic  wave,  but 
rather  an  oscillating  electric  field,  it  is  interesting  to  see 
that  the  repeating  unit  of  five  residues  acts  as  a  single 
unit  with  an  oscillating  net  dipole  that  displays  mechan¬ 
ical  resonance  in  a  similar  frequency  range.  It  may  be 
further  noted  in  Figure  7B  that  the  absorption  intensity 
increases  by  more  than  50%  upon  going  from  20  to 
40°C.  Because  of  this,  the  acoustic  absorption  exhibited 
by  (GVGIP)„  at  20°C  in  Figure  7A  would  be  expected  to 
increase  by  an  additional  50%  upon  increasing  the 
temperature  further,  to  40°C.  This  would  make  the 
absorption  per  unit  volume  of  (GVG IP)„  much  larger 
than  that  exhibited  by  natural  rubber  or  urethane. 


(ii)  Mechanical  resonances  in  the  MHz  range :  Dielec¬ 
tric  relaxation  studies  demonstrate  mechanical  reso¬ 
nances  near  5  MHz  for  (GVGVP)„,  (GVGIP)„, 
(GVGLP)„  (Buchet  et  al.,  1988),  a-elastin  and  fibrous 
elastin  purified  from  bovine  ligamentum  nuchae.  The 
data  of  Figure  8A  give  the  real  part  of  the  dielectric 
permittivity  (s')  for  a-elastin,  a  70  kDa  chemical  frag¬ 
mentation  product  of  natural  fibrous  elastin  (Urry  et  al., 
1988a)  and  for  the  elastin  model  poly(GVGVP)  (inset) 
(Henze  and  Urry,  1985).  e"  for  poly(GVGVP)  and 
fibrous  elastin  is  shown  in  Figure  8B  (Luan  et  al,  1988). 

As  seen  in  Figure  1,  there  is  much  more  to  elastin  than 
poly(GVGVP),  of  the  W4  sequence,  although  it  is  the 
most  prominent  repeating  sequence  in  bovine  elastin. 
The  relaxations  for  100%  poly(GVGVP)  are  expected  to 


Fig.  8.  Dielectric  relaxation  spectra  in  the  1-1000  MHz  frequency  range  of  a-elastin  (part  A),  the  poly(GVGVP)  of  elastin  (inset  of  part  A  and 
part  B),  and  fibrous  elastin  (part  B).  In  all  cases  as  the  temperature  is  raised  there  develops  an  intense  mechanical  relaxation  centered  near  5  MHz 
as  repeating  elements  of  these  protein  systems  hydrophobically  fold  into  regular,  albeit  obviously  dynamic,  structures.  In  the  case  of 
poly(GVGVP)  each  pentamer  folds  into  the  same  conformation,  in  which  the  peptide  moieties  (the  only  entities  with  dipole  moments)  undergo 
coordinated  rocking  motions  with  a  resulting  oscillating  mean  dipole  moment.  The  rocking  of  the  mean  dipole  moments  of  the  pentamers 
resonate,  in  this  case  move,  at  the  same  frequency  in  response  to  an  alternating  electric  field,  that  is,  they  exhibit  a  mechanical  resonance  centered 
near  5  MHz.  A  similar  mechanical  resonance  is  observed  near  1  kHz  in  Figure  7.  Part  A  reproduced  with  permission  from  Urry  et  al.  (1988a), 
and  part  B  reproduced  with  permission  from  Luan  et  at.  (1988). 
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be  more  intense  than  that  of  oc-elastin  and  fibrous 
elastin.  It  is  indeed  surprising  that  the  frequency  overlap 
and  intensity  is  as  close  as  found.  It  is  apparent  that 
other  repeating  sequences,  or  quasi  repeats,  exhibit 
dielectric  relaxations  at  slightly  lower  frequencies. 

Importantly,  mechanical  resonances  near  5  MHz  and 
1  kHz  indicate  the  presence  of  regularly  repeating 
dynamic  structures.  Extension  of  such  structures  rea¬ 
sonably  gives  rise  to  an  entropic  elastic  force  upon 
extension  by  the  damping  of  internal  chain  dynamics 
represented  by  mechanical  resonance.  Molecular  me¬ 
chanics  and  dynamics  calculations  based  on  the  molec¬ 
ular  structure  shown  in  Figure  3  will  demonstrate  just 
how  effectively  the  structure  can  explain  the  experimen¬ 
tal  elasticity  and  relaxation  data. 

Entropy  calculations  based  on  the  elastin  model f 
(GVGVP)n 

Calculations  of  structures  and  the  entropies  of  the 
structures  have  used  conventional  molecular  mechanics 
and  dynamics  programs  such  as  CHARMM  (Karplus 
and  McCammon,  1983;  Karplus  and  Kushick,  1981), 
AMBER  (Weiner  and  Kollman,  1981),  and  Scheraga’s 
ECEPP  (Momany  et  al,y  1974,  1975).  The  satisfying 
feature  of  these  quite  distinct  computational  approaches 
is  that  they  give  essentially  identical  results  whether 
carried  out  within  our  group  by  different  researchers  or 
by  different  groups.  It  is  important,  however,  that  the 
constraints  and  boundary  conditions  should  properly 
reflect  the  conditions  considered.  Paramount  among  the 
constraints  is  adherence  to  the  elasticity  experiment  in 
which  the  ends  of  the  molecular  structure  are  fixed  at 
whatever  selected  degree  of  extension. 

Molecular  mechanics  calculations  using  Scheraga’s 
ECEPP  approach  for  energy  surfaces  in  <)>-\[/ 
configuration  space 

(i)  Calculation  of  entropy  change  during  extension  by  an 
enumeration  of  states  approach  for  relaxed  and  extended 
states'.  In  Scheraga’s  approach  the  internal  energy  of  a 
chosen  chain  segment  (in  our  case  the  pentamer  permu¬ 
tation,  V1P2G3V4G5)  is  calculated  as  a  function  of  a  pair 


of  adjacent  torsion  angles.  Normally,  the  4>  and  ^ 
torsion  angles  are  identical  to  those  in  the  Ramachan- 
dran  plot.  In  our  case  since  the  primary  sites  of  motion 
are  the  peptide  moieties  of  the  suspended  segment 
(V4G5V,),  such  that  the  two  /  plots,  now  called 
lambda  plots,  become  $5-^4  and  <t>i-\J/5. 

In  the  enumeration  of  states  approach,  a  5°  change  in  a 
single  torsion  angle  is  counted  as  a  new  state,  and  the 
number  of  states  is  counted  for  the  chosen  cut-off  energy 
for  both  the  relaxed  and  130%  extended  structures  (Urry 
et  al ,  1982;  Urry  and  Venkatachalam,  1983;  Urry  et  ai , 
1985c;  Urry,  1991).  A  single  pentamer  is  calculated 
within  the  relaxed  and  extended  (J-spiral  structures,  and 
0.6,  1.0  and  2.0  kcal/mol-pentamer  cut-off  energies  are 
used  as  well  as  an  energy  weighting  using  the  Boltzmann 
summation  over  states.  Table  1  gives  the  number  of 
states  where  the  entropy  is  simply  calculated  by  Equation 
(9)  using  the  example  of  1  kcal/mol-pentamer, 

AS  =  (S6  -  S')  =  R\n(Wc/WT)  =  R ln(58/762) 

=  -5.1  cal/mol-pentamer.  (9) 

This,  of  course,  is  -1.02  cal/mol-residue  as  given  in 
Table  1. 

(ii)  Calculation  of  dielectric  relaxation  data  from 
pentamer  molecular  mechanics :  With  the  aid  of  the 
Onsager  equation  for  polar  liquids  (Onsager,  1931; 
Bottcher,  1973)  it  becomes  possible  to  use  the  structures 
obtained  above  and  to  sum  the  dipole  moments  of  the 
states  of  Table  1 .  In  addition,  the  same  can  be  done  for 
V1P2G3V4A’5)  where  the  A’  stands  for  the  D-Ala  amino 
acid  residue  (Venkatachalam  and  Urry  1986  and  print¬ 
ers  erratum  as  corrected  below;  Urry,  1991).  Compar¬ 
ison  of  experimental  and  calculated  values  is  particularly 
satisfying.  ‘In  view  of  this  result,  it  seems  reasonable  to 
consider  a  cutoff  energy  of  1.5  kcal  moP1.  This  value  of 
energy  cutoff  leads  to  a  mean  dipole  moment  change  of 
3.8  Debye  per  pentamer  in  good  agreement  with  the 
value  obtained  from  dielectric  relaxation  studies.  Similar 
dielectric  studies  on  D.Ala5-polypentapeptide  indicate  a 
dipole  moment  of  about  one-third  of  that  found  for  the 
polypentapeptide  at  37°C. . .  This  is  in  very  good 
agreement  with  the  dipole  changes  obtained  from 


Table  /.  Perspective  of  entropy  of  the  poly(GVGVP)  p-spiral  by  enuneration  of  states 


Number  of  states 

Cutoff  energy  (kcal/mol) 

Relaxed  (r) 

Extend  (e) 

Entropy  change  per  residue 

2 

1853 

162 

-0.97 

1 

762 

58 

-1.02 

0.6 

342 

Using  Boltzmann  sum  over  states 

/= 2>~s/‘r 

i 

24 

-1.06 

2 

AS  =  R  ln(/* /f)  +  —y— 

-1.01 
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molecular  mechanics  calculations  presented  here’  (Venk- 
atachalam  and  Urry  1986). 

The  above  calculated  results  utilize  the  5  MHz  me¬ 
chanical  resonance.  They  further  demonstrate  how 
readily  calculations  based  on  the  dynamic  structure 
given  in  Figure  3  calculate  experimental  dielectric  relax¬ 
ation  data  with  inherent  relevance  to  entropic  elastic 
force.  The  values  for  entropy  change  calculated  using 
the  molecular  mechanics  approach,  when  used  in  the 
second  term  of  Equation  (1),  provide  excess  entropic 
elastic  force  as  expected  (see  Urry  et  al. ,  2002a,  b). 

Molecular  dynamics  calculations  using  Karplus ’ 
CHARMm  and  Kollman  s  AMBER  programs 
As  the  motions  of  the  torsional  oscillations  are  dynamic, 
it  is  appropriate  to  consider  molecular  dynamics  calcu¬ 
lations  using  the  Karplus  software  program  adapted  by 
Polygen,  Inc.  (now  Molecular  Simulations),  called 
CHARMm,  version  20.3.  For  these  Newtonian  classical 
mechanics  simulations  the  potential  energy  functions  and 
parameters  suggested  by  the  Karplus  group  (Brooks 
et  al. ,  1983)  were  used  in  a  software  program.  The  root- 
mean-square  (RMS)  fluctuations  of  the  torsion  angles, 
A<t>,-  and  Av| /,-  for  the  relaxed  and  extended  (130%)  of 
(GVGVP)i  i  as  given  in  Table  2  are  used  in  Equation  (10). 

AS  =  R  lntn/A^A^/n/A^A^]  (10) 

The  calculated  change  in  entropy  upon  extension  of 
-1.1  cal/mol  deg-residue  is  in  remarkable  agreement 
with  the  previously  described  molecular-mechanics  cal¬ 
culations. 


Furthermore,  Wasserman  and  Salemme  (1990)  using 
the  AMBER  software  program  of  the  Kollman  group 
(Weiner  and  Kollman,  1981),  but  also  including  a 
medium  of  water  molecules,  obtained  essentially  iden¬ 
tical  results.  It  must  be  concluded,  damping  of  internal 
chain  dynamics  is  an  abundant  source  of  decrease  in 
chain  entropy  that  is  sufficient  to  account  for  the 
entropic  elastic  force. 

As  Wasserman  and  Salemme  (1990)  included  water 
molecules  in  their  calculations,  and  an  ordering  of  water 
molecules  was  found  as  the  hydrophobic  side  chains 
became  exposed  on  chain  extension,  they  considered  this 
decrease  in  solvent  entropy  as  a  possible  source  of  the 
entropic  elastic  force.  As  will  be  shown  below,  solvent 
entropy  change  does  not  make  a  significant  contribution 
to  the  entropic  elastic  force  development  during  isomet¬ 
ric  contractions. 


Relationship  between  hydrophobic  association-dissocia¬ 
tion,  elastic  force  and  energy  conversion  in  elastin 
mechanics 

By  introducing  glutamate  (Glu,  E)  the  ideal  elastomer  in 
Figure  5A  becomes  a  chemo-mechanical  transducer 
capable  of  converting  the  chemical  energy  of  proton 
concentration  changes  into  the  mechanical  energy  of 
isotonic  and  isometric  contractions  (See  Figure  9A  and 
B).  An  analysis  of  this  data  in  combination  with  that  of 
Figure  6B  resolves  the  issue  of  the  role  of  changes  in 
solvent  entropy  in  elastic  force  development  and  energy 
conversion. 


Table  2.  RMS  fluctuations  of  torsion  angles  (cj>  and  \|/)  of  (VPGVG)n  (45  ps  of  equilibration  time  and  80  ps  of  molecular  dynamics  simulation). 
Reproduced  with  permission  from  Chang  and  Urry  (1989)  _ _ _ 
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Mcchano-chcmical  Coupling  in  Cross-linked  4%  Clu-Poly  penlapcptidc 
Due  to  Changes  in  pH  at  37X 


Fig.  9.  Mechano-chemical  coupling  by  20  Mrad  y-irradiation  cross-linked  of  poly[0.8(GVGVP),00.2(GEGVP)]  under  isotonic  (A)  and  isometric 
(B)  conditions.  In  part  B  under  isometric  contraction  conditions  there  occurs  an  increase  in  entropic  elastic  force  resulting  from  a  decrease  in 
entropy  of  the  elastomer  while  there  occurs  an  increase  in  the  entropy  of  the  solvent  as  hydrophobic  hydration  becomes  bulk  water.  Therefore  the 
increase  in  entropic  elastic  force  cannot  be  the  result  of  the  solvent  entropy  change.  Reproduced  with  permission  from  Urry  et  al.  (1988b). 


Chemically  driven  isotonic  (dLjd[iH)f  and  isometric 
(W^h)l  contractions 

The  elastin  model  protein  of  interest  is  poly  [xv- 
(GVGVP),xe(GEGVP)],  where  xv  and  xE  are  the  mole 
fractions  with  *v  +  xE—\  and  for  the  case  where  xE= 
0.2.  Upon  y-irradiation  cross-linking  an  elastic  matrix  is 
formed  which  can  be  studied  under  isotonic  and 
isometric  conditions.  Results  are  shown  in  Figure  9A 
and  B  (Urry  et  al .,  1988b).  The  isotonic  contraction  (0L/ 
0p//)y-  in  Figure  9A  is  the  change  in  length  at  constant 
load  (force)  resulting  from  an  increase  in  the  concen¬ 
tration  of  proton,  that  is,  the  increase  in  chemical 
potential  of  proton  (0p./f).  The  isometric  contraction  (0// 
0p H)L  in  Figure  9B  is  the  change  in  force  at  constant 
length  due  to  an  increase  in  proton. 

Both  of  these  contractions  result  from  hydrophobic 
association  due  to  protonation  of  the  carboxylate  group. 
The  underlying  physical  process  is  the  competition  for 
hydration  between  the  charged  carboxylate  and  the 
hydrophobic  side  chains  of  the  valine  (Val,  V)  residues 
(Urry,  1997).  This  competition,  described  as  an  apolar- 
polar  repulsive  free  energy  of  hydration,  results  in 


increasing  positive  cooperativity  of  the  acid-base  titra¬ 
tion  curves  that  correlate  with  increasingly  larger  ptfa 
shifts.  As  Val  residues  are  replaced  by  more  hydropho¬ 
bic  phenylalanine  (Phe,  F)  residues,  the  p K&  shifts 
increase  (up  to  6  pH  units  or  more)  and  positive 
cooperativity  increases  with  Hill  coefficients  of  up  to  8 
or  more  (Urry,  1997).  The  apolar-polar  repulsive  free 
energy  of  hydration  is  related  to  the  change  in  Gibbs 
free  energy  for  hydrophobic  hydration,  AGHa(x)»  °f  a 
phosphate  as  discussed  in  relation  to  Equations  (17)  and 
(19)  below,  but  first  the  source  of  entropic  elastic  force  is 
addressed. 

Isometric  contractions,  (df/dT)L  and  (df/d\xH)L>  confirm 
basis  for  entropic  elastic  force 

From  Equation  (1)  the  entropic  component  of  the  elastic 
force,  fs ,  is  proportional  to  -AS,  that  is,  an  increase  in 
entropic  elastic  force  results  from  a  decrease  in  entropy. 
Also  the  formation  of  hydrophobic  hydration  from  bulk 
water  constitutes  an  inherently  negative  AS,  and,  of 
course,  the  loss  of  hydrophobic  hydration  constitutes  a 
positive  change  in  entropy.  This  change  in  solvent 
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entropy  due  to  formation  of  hydrophobic  hydration 
has  been  credited  as  an  important  contribution  to 
entropic  elastic  force  (Weis-Fogh  and  Anderson,  1970; 
Wasserman  and  Salemme,  1990,  Alonso  et  al. ,  2001). 
Therefore,  the  immediate  question  becomes  whether 
experimental  results  support  or  eliminate  a  decrease  in 
solvent  entropy  as  a  source  of  entropic  elastic  force. 

Thermally  driven  isometric  contraction  ( Figure  6B) 
Figure  6B  shows  the  classic  thermoelasticity  study  for 
determining  the  /e// ratio  for  y-irradiation-cross-linked 
poly(GVGVP).  Because  extension  is  fixed,  a  thermally 
driven  isometric  contraction  also  occurs  over  the  tem¬ 
perature  range  of  the  inverse  temperature  transition  of 
hydrophobic  association.  During  hydrophobic  associa¬ 
tion,  hydrophobic  hydration  becomes  bulk  water. 
Therefore,  the  solvent  entropy  change  is  positive.  Any 
contribution  to  the  entropic  elastic  force  due  to  changes 
in  hydrophobic  hydration  during  this  thermally  driven 
isometric  contraction  would  necessarily  result  in  a 
decrease  in  force.  Yet  the  experimental  result  of  Fig¬ 
ure  6B  shows  the  development  of  an  elastic  force  that  is 
90%  entropic,  while  the  change  in  solvent  entropy  is  of 
the  opposite  sign.  Clearly,  the  change  in  solvent  entropy 
is  not  contributing  to  the  estimated  90%  entropic  elastic 
force.  Having  eliminated  solvent  entropy  change  as  the 
source  of  entropic  elastic  force  during  a  thermally  driven 
isometric  contraction,  we  now  analyze  the  solvent 
entropy  change  attending  a  chemically  driven  isometric 
contraction. 

Chemically  driven  isometric  contraction  ( Figure  9B) 

As  seen  in  Figure  9B,  at  fixed  extension  of  the  y- 
irradiation-cross-linked  elastic  matrix  comprised  of 
poly[0.8(GVGVP),0.2(GEGVP)],  protonation  of  four 
carboxylates  per  100  residues  results  in  development  of 
elastic  force.  A  thermoelasticity  characterization  of  this 
matrix  at  low  pH  gives  the  same  result  of  dominantly 


entropic  elasticity  as  found  in  curve  b  of  Figure  6B  for 
poly(GVGVP)  in  the  absence  of  carboxyl  moieties. 

The  process  of  protonation  allows  reconstitution  of 
hydrophobic  hydration  to  such  an  extent  that  the 
temperature  range  for  hydrophobic  association  drops 
below  that  of  the  operating  temperature  (Urry,  1993, 
1997).  The  result  is  a  contraction  due  to  hydrophobic 
association.  Again,  during  an  isometric  contraction  (this 
time  chemically  driven),  hydrophobic  hydration  be¬ 
comes  less  ordered  bulk  water.  The  solvent  entropy 
increases  during  the  development  of  entropic  elastic 
force  due  to  a  decrease  in  entropy. 

In  addition  to  eliminating  solvent  entropy  change  as  a 
source  of  entropic  elastic  force,  the  isometric  contrac¬ 
tion  results  provide  the  conceptual  bridge  between  the 
fundamental  source  of  entropic  elastic  force  and  con¬ 
traction  by  hydrophobic  association. 

Hydrophobic  association  effects  elastic  force  development 
by  extending  ( and  thereby  damping  internal  dynamics  of) 
interconnecting  chain  segments 

Clam-shaped  globular  proteins  that  open  and  close  due 
to  hydrophobic  association 

The  top  part  of  Figure  10  shows  a  series  of  clam-shaped 
globular  proteins  strung  together  near  the  open  end  by 
elastic  bands  and  maintained  at  fixed  extension.  In  this 
isometric  state,  there  is  depicted  an  equilibrium  between 
open  and  closed  states.  Obviously,  shifting  the  equilib¬ 
rium  toward  the  closed  state  would  increase  the  force 
measured  at  the  force  transducer,  whereas  shifting 
toward  the  open  state  would  lower  the  force  on  the 
interconnecting  elastic  bands.  Thus,  any  process  that 
increased  hydrophobic  association  would  increase  the 
elastic  force.  Now  in  the  lower  part  of  Figure  10,  we 
relate  this  perspective  to  the  p-spiral  representation  of 
Figure  3C. 


hydrophobic  hydrophobic  hydrophobic  association 

association  surface  between  turns 


delineating  stretched  chain  segments 


Fig.  JO.  Cartoon  of  the  relationship  between  hydrophobic  association  and  entropic  elastic  force  development.  Above:  A  series  of  clam-shaped 
globular  protein  strung  together  by  elastic  bands  with  an  equilibrium  between  open  and  hydrophobically  associated  closed  states.  Clearly,  as  the 
equilibrium  shifts  toward  more  closed  states,  the  force,  sustained  by  the  interconnecting  elastic  segments,  increases.  Below:  Representation  of  the 
p-spiral  structure  of  poly(GVGVP)  at  an  intermediate  state  of  hydrophobic  association  between  turns  of  the  p-spiral.  Clearly,  as  more 
hydrophobic  association  occurs  between  turns,  the  interconnecting  elastic  chain  segments  become  further  extended  with  an  increase  in  damping  of 
internal  chain  dynamics  giving  rise  to  greater  entropic  elastic  force.  Reproduced  with  permission  from  Urry  (1990). 
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Schematic  representation  of  an  equilibrium  between 
hydrophobically  associated  and  dissociated  turns 
of  a  spiral  held  at  fixed  extension 
The  bottom  half  of  Figure  10  uses  the  representation  of 
the  p-spiral  given  in  Figure  3C  to  depict  an  equilibrium 
between  associated  and  dissociated  turns  of  the  p-spiral. 
As  more  of  the  turns  are  involved  in  hydrophobic 
association,  those  not  involved  become  extended  to  a 
greater  extent.  Thus,  as  more  hydrophobic  association 
occurs,  the  interconnecting  segments  become  increas¬ 
ingly  extended.  By  the  entropic  mechanism  of  the 
damping  of  internal  chain  dynamics  on  extension,  the 
entropic  elastic  force  increases. 

Opening  and  closing  of  the  clam-shaped  globular 
proteins  of  the  upper  part  of  Figure  10  (Urry,  1990) 
would  be  detectable  under  conditions  of  force-clamp 
atomic  force  spectroscopy.  Such  traces  have  been 
observed  using  titin-based  elastic  constructs  comprised 
of  a  series  of  hydrophobically  folded  p-barrels  (Ober- 
hauser  et  al ,  2001).  Because  of  the  continuous  nature  of 
increases  in  hydrophobic  association  between  turns  of 
the  p-spiral  represented  in  the  lower  part  of  Figure  10, 
such  force-clamp  traces  for  poly(GVGVP)  would  not 
have  such  steps,  unless  sensitivity  were  sufficient  to 
detect  the  addition  of  a  repeating  unit  to  hydrophobic 
association. 

Thus,  the  changes  in  entropic  elastic  force  and  in 
hydrophobic  association  are  two  distinct  physical  pro¬ 
cesses  that  become  interlinked  when  occurring  within 
the  same  chain  system.  Hydrophobic  association  causes 
the  extension  of  interconnecting  dynamic  chain  seg¬ 
ments  with  the  result  of  the  damping  of  internal  chain 
dynamics  within  those  non-hydrophobically  associated 
interconnecting  chain  segments. 

Derivation  of  Gibbs  free  energy  for  hydrophobic 
association  A GHA 

The  role  of  hydrophobic  association  in  the  development 
of  entropic  elastic  force  has  been  delineated  above.  The 
next  step  is  to  derive  an  expression  for  the  Gibbs  free 
energy  of  hydrophobic  association  and  to  relate  it  to 
processes  that  direct  elastic-contractile  processes. 

(i)  At  the  temperature  of  the  inverse  temperature 
transition ,  pH a  =  Phd  and  AG  =0:  The  chemical  po¬ 
tential  (p)  is  the  Gibbs  free  energy  per  mole  (AG /An).  By 
the  Gibbs  phase  rule,  at  the  temperature  of  the  inverse 
temperature  transition  for  hydrophobic  association  the 
chemical  potential  of  the  hydrophobically  disassociated 
(dissolved)  state  (phd)  and  the  chemical  potential  of  the 
hydrophobically  associated  state  (pha)  are  equal.  There¬ 
fore,  at  the  temperature  of  the  inverse  temperature 
transition  (7\),  AGt  =  A Ht  -  TtASt  =  0,  where  A Ht  is 
the  heat  of  the  transition  and  ASt  is  the  entropy  change 
for  the  transition.  Accordingly  at  the  value  of  Tt  for  a 
given  elastic  model  protein  composition,  AHt=TtASt. 

By  way  of  example,  consider  the  two  different  compo¬ 
sitions  of  model  elastic  proteins,  (GVGVP)„  and 
(GVGIP)„.  These  two  polypentapeptides  differ  only  by 


a  single  CH2  moiety  per  pentamer.  In  particular,  the  R- 
group  for  the  V  residue,  — CH(CH3)2,  differs  from  the 
R-group  of  the  I  residue,  —CH(CH3)— CH2— CH3,  by 
the  addition  of  a  single  CH2  moiety.  At  the  temperature 
for  each  of  their  respective  phase  transitions  where 
Hha=Phd>  we  can  write, 

A//t(GVGVP)  =  7i  (G  V  GVP)  ASt  (GV  GVP)  (11) 

and 

A//t(GVGIP)  =  Tt(GVGIP)ASt(GVGIP).  (12) 

Butler  (1937)  showed  that  each  addition  of  a  CH2 
moiety  in  the  series  of  normal  alcohols  from  methanol  to 
/i-pentanol  on  dissolution  in  water  resulted  in  an  average 
exothermic  heat  (A/^)/CH2  of -1.4  kcal/mol-CH2  and  a 
change  of  the  term,  (~TAS)/CH2,  of  a  +1.7  kcal/mol- 
CH2.  While  there  are  many  important  points  to  make 
about  this  fundamental  finding,  we  use  it  here  simply  as 
a  justification  to  rewrite  Equation  (12)  as 

Atft(GVGIP)  -  A//t(GVGVP)  +  A//t(CH2),  (13) 

and 

7i  (GV GIP)  ASt  (GV GIP)  =  Tt  (GV GIP)  [A5t  (GV GVP) 

+  AiSt(CH2)].  (14) 

Substituting  Equations  (13)  and  (14)  into  Equation  (12) 
and  subtracting  Equation  (1 1)  gives 

[A//t(CH2)  -  Tt(GVGIP)A5t(CH2)] 

=  [rt(GVGIP)  -  Tt (GV GVP)]  ASt  (GV GVP) .  (15) 

The  left  hand  side  of  Equation  (15)  is  recognized  as  an 
expression  for  the  change  in  Gibbs  free  energy  due  to 
addition  of  the  CH2  moiety.  On  consideration  of  the 
phase  transition  being  analyzed,  the  left  hand  side  of 
Equation  (5)  becomes  the  change  in  Gibbs  free  energy  of 
hydrophobic  association  due  to  the  addition  of  a  single 
CH2  moiety,  AGHa(CH2),  that  causes  the  inverse 
temperature  transition  of  hydrophobic  association  to 
occur  at  a  lower  temperature  for  (GVGIP)„. 

AGHa(CH2)  =  [Tl(GVGIP)  -  7HGVGVP)] 

x  ASt(GVGVP).  (16) 

The  right  hand  side  of  Equation  (16)  is  negative  because 
rt(GVGIP)  <  rt(GVGVP)  and  ASt(GVGVP)  is  posi¬ 
tive  for  a  transition  whereby  hydrophobic  hydration 
becomes  less  ordered  bulk  water  during  hydrophobic 
association,  that  is,  as  insolubility  ensues. 

The  addition  of  a  CH2  moiety  per  pentamer  is  one  of 
very  many  ways  whereby  the  value  of  Tt  can  be  changed. 
A  rt-based  hydrophobicity  scale  has  been  developed  for 
substitution  of  each  of  the  naturally  occurring  amino 
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acid  residues  and  chemical  modifications  thereof  by 
systematically  increasing  the  mole  fraction  of  pentamers 
containing  the  change  X  and  extrapolating  to  the  value 
of  Tx  that  would  occur  for  poly(GXGVP)  (Urry,  1997). 
Therefore,  the  value  of  Tt  may  be  greater  than  100°C  as 
for  poly(GEGVP)  when  the  glutamic  acid  residue,  E,  is 
ionized  as  the  carboxylate,  where  7t= 250°C.  On  the 
other  hand,  the  value  of  Tt  may  be  less  than  0°C  as  for 
poly(GWGVP)  when  the  tryptophan  is  the  residue, 
where  rt=-90°C.  For  the  amino  acids  with  functional 
side  chains  the  revalue  changes  upon  changing  the  state 
of  the  side  chain  from  carboxyl  to  carboxylate,  or  the 
redox  state  of  an  attached  redox  couple.  The  most 
dramatic  change  in  the  revalue  occurs  upon  phospho¬ 
rylation  of  a  Ser,  Thr,  or  Tyr  residue  (Pattanaik  et  aly 
1991).  Furthermore,  adding  salt  to  the  solution  of  a 
neutral  polymer  such  as  poly(GVGVP)  itself  lowers  the 
value  of  Tx  (Urry,  1993),  but  the  A7\  is  10  times  greater 
when  the  salt  ion-pairs  with  a  charged  side  chain  (Urry 
and  Luan,  1995a;  Urry  et  a/.,  1997a). 

So  we  generalize  to  any  experimental  variable,  that 
alters  the  value  of  Tx  of  a  reference  model  elastic  protein 
and  write, 

AGha(x)  =  [?t(x)  -  rt(reference)]A5t (reference), 

(17) 

or 

AGha(x)  =  Ar,(x)ASt(reference),  (18) 

where  Tx  is  the  value  of  a  reference  polymer  and  state 
has  been  changed  to  Tx(x)  by  the  experimental  variable 
%.  All  of  these  energy  inputs  can  now  be  described  in 
terms  of  their  effect  on  the  free  energy  of  hydrophobic 
association,  AGHa(x)* 

Relevance  of  Gibbs  free  energy  for  hydrophobic 
association  to  energizing  by  phosphorylation 

Phosphorylation  of  poly[30(GVGIP),(RGYSLG)]  by 
the  cardiac  cyclic  AMP-dependent  protein  kinase  at 
the  serine  (Ser,  S)  residue  of  the  polymer  goes  to  an 
average  of  47%  completion  by  the  following  reaction 
(Pattanaik  et  aly  1991;  A.  Pattanaik,  personal  commu¬ 
nication). 

ATP  +  poly[30(GVGIP),  (RGYSLG)] 

=  ADP  +  poly[30(GVGIP), 

(RGYS(— OPOf-)LG)l.  (19) 

The  value  of  Tx  for  this  polymer  before  phosphorylation 
was  approximately  20°C  and  the  change  in  the  value  of 
Tx  extrapolated  as  defined  to  a  reference  state  of  one 
phosphate  per  GVGIP  becomes  860°C,  which 
is  represented  as  rt(— OPOj-).  If  we  now  use  the  experi¬ 
mentally  derived  value  of  ASt(GVGIP)  =  9.32  cal/mol 
(pentamer)  deg  as  the  ASt(reference)  (Luan  and  Urry, 


1999)  of  Equation  (17),  the  change  in  Gibbs  free  energy 
for  hydrophobic  association  due  to  phosphorylation, 
AGHa(“”OP03-),  is  +7.8  kcal/mol.  For  the  reaction  of 
Equation  (19)  with  an  equilibrium  constant  of  K  ~  47/ 
53  =  0.9  =  exp(-AG/R7),  the  change  in  Gibbs  free 
energy  for  the  reaction  is  0.06  kcal/mol.  Thus,  the 
calculated  value  for  the  free  energy  contributed  by  the 
terminal  phosphate  of  ATP  would  be  0.06-7.8  =  -7.7 
kcal/mol,  which  is  remarkably  close  to  the  free  energy  of 
hydrolysis  of  the  terminal  phosphate  of  ATP  of  7.3  kcal/ 
mol  (Voet  and  Voet,  1995). 

Accordingly,  phosphorylation  raises  the  free  energy  of 
the  hydrophobically  associated  state,  that  is,  it  results 
in  hydrophobic  disassociation,  and  phosphate  removal 
drives  hydrophobic  association,  e.g.,  contraction,  more 
effectively  than  any  other  chemical  change  measured  so 
far.  Expectations  are  that  the  binding  of  ATP,  and  most 
effectively  ADP  plus  phosphate,  would  dramatically 
raise  the  free  energy  of  the  hydrophobically  associated 
state  as  evidenced  by  an  increase  the  value  of  Tx .  This 
A7\(— OPO3-)  and  AGHa(“OP03_)  can  be  only  partially 
modulated  by  ion-pairing,  for  example,  with  magne¬ 
sium,  of  the  phosphates  at  the  binding  site. 

Relationship  between  entropic  elasticity  and  efficiency 
of  energy  conversion 

The  process  of  stretching  an  elastomer  constitutes  an 
energy  input  with  the  expended  mechanical  energy  yAL, 
where  /  is  force  and  A L  is  the  change  in  length. 
Accordingly,  the  mechanical  energy  is  obtained  from 
the  area  under  the  force  vs.  extension  curve.  For  an  ideal 
elastomer,  the  plot  of  force  vs.  extension  is  perfectly 
reversible,  as  shown  by  the  second  and  fifth  traces  of 
Figure  5A  for  the  single-chain  force-extension  curve  of 
the  entropic  elastomer,  (GVGVP)„  x  251-  Thus,  the  input 
energy  is  completely  recovered  on  relaxation  for  an  ideal 
(entropic)  elastomer.  Since  the  energy  recovered  on 
relaxation  can  be  viewed  in  the  sense  of  an  energy 
output,  an  ideal  elastomer  could  be  considered  a  perfect 
machine  for  storing  energy  of  deformation. 

Often  the  force  vs.  length  curve  obtained  on  relaxa¬ 
tion  falls  below  the  force  vs.  length  curve  obtained 
on  extension  (as  seen  in  the  lower  curve  of  Figure  5B); 
then  the  material  is  said  to  exhibit  a  hysteresis.  In 
this  case  the  energy  recovered  on  relaxation  is  less 
than  that  expended  on  extension.  The  input  deformation 
energy  has  been  dissipated  in  some  way  and  to  an  extent 
indicated  by  the  difference  in  the  areas  below  the 
extension  and  relaxation  curves. 

When  energy  conversion  occurs  by  means  of  an  ideal 
elastic  material,  it  is  possible  for  the  energy  conversion 
to  occur  at  high  efficiency,  whereas  when  energy 
conversion  occurs  using  an  elastic  material  that  exhibits 
hysteresis,  the  efficiency  of  energy  conversion  becomes 
limited  at  least  to  an  extent  determined  by  the 
magnitude  of  the  hysteresis.  The  elasticity  of  polymeric 
materials  becomes  inextricably  intertwined  with  the 
efficiency  of  energy  conversion.  Thus,  the  increase  in 
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elastic  force  during  the  isometric  contraction  and 
relaxation  of  Figure  9B  would  be  exactly  reversible  for 
an  ideal  elastomer  as  would  the  isotonic  contraction  of 
Figure  9A.  This  would  not  be  the  case  for  a  molecular 
machine  comprised  of  the  elastic  polymer  functioning  in 
Figure  5B. 

An  important  element,  therefore,  of  an  ideal  elasto¬ 
mer  is  to  have  the  energy  uptake  into  the  elastomer 
during  extension  reside  entirely  in  the  backbone  modes 
where  it  can  be  recovered  on  relaxation.  Should  energy  of 
deformation  find  its  way  into  side-chain  motional  modes 
and  into  chains  not  bearing  the  deformation  and 
irreversibly  into  solvent,  this  deformation  energy  be¬ 
comes  dissipated  and  unavailable  during  relaxation, 
resulting  in  hysteresis.  Comparison  of  the  single-chain 
force-extension  curves  of  Figure  5  for  poly(GVGVP) 
with  poly(GVGIP)  provides  an  example  with  proposed 
loss  of  energy  into  the  side  chain  motions  and  interac¬ 
tions  of  the  bulkier  isoleucine  (I)  residue  with  its  added 
CH2  moiety  into  adjacent  non-load-bearing  chains. 
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